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INTRODUCTION 

1. FOREWORD 

The elucidation of the mechanism of electroweak symmetry breaking is one of the main goals of the LHC 
physics program. In the Standard Model (SM), mass generation is triggered by the Higgs mechanism, 
which predicts the existence of one scalar state, the Higgs boson [1,2]. The Higgs boson couplings to 
fermions and gauge bosons are a prediction of the model and the only unknown parameter is the Higgs 
boson mass. 

The Minimal Supersymmetric extension of the Standard Model (MSSM) requires the introduction 
of two Higgs doublets, in order to preserve supersymmetry and give mass to the fermions, and after 
spontaneous symmetry breaking five Higgs particles remain in the spectrum: two CP-even {h,H), one 
CP-odd {A) and two charged {H^) Higgs bosons. At lowest order the MSSM Higgs sector can be 
described by two parameters, generally chosen to be niA, the mass of the pseudoscalar Higgs boson, and 
tan/? = V2/V1, the ratio of the two vacuum expectation values. The lowest order predictions receive 
large radiative corrections which must be included when calculating Higgs couplings or masses. At tree 
level, the lightest neutral Higgs boson has an upper bound of Mz, which is increased to nih ^ 0(130 — 
140) GeV when radiative corrections are included [3]. 

The search for the Higgs at collider experiments has now being on-going for two decades. The 
present direct lower limit of the Higgs mass in the SM is 114.4 GeV (at 95% CL) [4], while precision 
measurements point to a rather light Higgs, nih ^ 180 GeV [5, 6]. The Tevatron has a chance to find 
evidence for a Higgs boson if enough integrated luminosity can be accumulated. At present, the Tevatron 
is performing well, and it is approaching the sensitivity limit required to exclude the existence of the SM 
Higgs for ruh ~ 160 GeV [7]. 

If it exists, the Higgs boson will be seen at the LHC, which can provide a measurement of the 
Higgs mass at the per-mille level and of the Higgs boson couplings at the 5-20 % level [8]. These 
tasks, however, require accurate theoretical predictions for both signal and background cross sections 
and distributions, and this is true in particular for an accurate determination of the properties of the 
discovered particle, such as spin, CP, and couplings. 
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In the following we review the status of theoretical predictions for both signal and background at 
the LHC, with emphasis on recent developments for the Standard Model Higgs boson. 

1.1 Gluon-Gluon Fusion 

The gluon fusion mechanism, mediated by a (heavy)-quark loop, provides the dominant production 
mechanism of Higgs bosons at the LHC in the full mass range. 

QCD corrections to this process at next-to-leading order (NLO) have been known for some 
time [9-11] and their effect increases the leading order (LO) cross section by about 80-100%. This 
calculation is very well approximated by the large-mtop limit. When the exact Bom cross section (with 
full dependence on the masses of top and bottom quarks) is used to normalize the result, the difference 
between the exact and the approximated NLO cross sections ranges from 1 to 4% when m/j < 200 
GeV. In recent years, the next-to-next-to-leading order (NNLO) corrections have been computed in this 
limit [12-17], leading to an additional increase of the cross section of about 10 — 15%. Soft-gluon 
resummation leads to a further increase of about 6% [18]. The latter result is nicely confirmed by the 
more recent evaluation of the leading soft contributions at N'^LO [19-21]. Two loop EW effects are also 
known [22-24]. 

In the MSSM, for large tan f3, the contribution from bottom quark loops becomes important and 
the large-mtop limit is not applicable. The full SUSY-QCD corrections are known in the limit of heavy 
squark and gluino masses at NLO [25-27]. Recently, the exact contribution of squark loops has been 
evaluated at NLO [28] and is discussed in Sect. II l.l The massive virtual corrections to the squark loops 
are given in Ref. [29, 30]. 

The higher order calculations mentioned above are certainly important but they refer to total cross 
sections, i.e., the experimental cuts are largely ignored. The impact of higher order corrections on the rate 
and the shape of the corresponding distributions may be strongly dependent on the choice of cuts. In the 
case in which one [31] or two [32] jets are tagged at large pT the NLO corrections for Higgs production 
from gluon fusion are known and implemented in parton level Monte Carlo programs. These predictions 
are obtained in the large mtop limit, which is a good approximation for small transverse momentum of 
the accompanying jet. For Higgs plus one jet production, there is a rather flat dependence of the K factor 
on pt and rapidity for moderate pT and y. In the MSSM, the Higgs plus 1 jet rate is known at lowest 
order only [33, 34]. The Higgs plus 2 jet process from gluon fusion is a background for vector boson 
fusion, as discussed below. Interference effects in the Higgs plus 2 jet channel are discussed in Sect. |4] 

The NNLO inclusive cross section when a jet veto is applied [35] has been known for some time. 
The first NNLO calculation that fully takes into account experimental cuts was reported in Ref. [36], 
in the case of the decay mode /i ^ 77 which is implemented in the FEHIP Monte Carlo program. 
In Ref. [37] the calculation was extended to the decay mode h W^W" l^l^vV. Recently, an 
independent NNLO calculation has been performed [38, 39], including all the relevant decay modes of 
the Higgs boson: h —>■ 77, h —>■ W^W~ —>■ l^l^vV and h —>■ ZZ —>■ 4 leptons. Such a calculation is 
implemented in a Monte Carlo program and is documented in this report in Sect. |2] 

Among the possible differential distributions, an important role is played by the transverse mo- 
mentum (j)t) spectrum of the Higgs boson [40]. When px ~ m/j the standard fixed order expansion is 
applicable. When px ^ mh, large logarithmic contributions appear that may invalidate the customary 
fixed order expansion. The resummation of such contributions has been performed at different levels of 
theoretical accuracy [41-46]. In Refs. [44-46] the resummed result up to next-to-next-to-leading loga- 
rithmic accuracy is matched to the fixed order NLO result [31,47,48] valid at large transverse momenta. 
It is important to note that transverse-momentum resummation is approximately performed by standard 
Monte Carlo event generators. A comparison of results obtained with different tools was presented in 
Ref. [49]. 

For Higgs boson masses below about 140 GeV the dominant decay mode /i ^ 66 is swamped by 



the huge QCD background and the Higgs boson can be found by looking at the rare /i ^ 77 decay mode. 
The 77 background can be measured precisely from the data using sideband interpolation, but accurate 
theoretical predictions are useful to estimate the expected accuracies and to better understand detector 
performances. The /i — > 77 decay width is known including full two loop QCD and EW effects [50-52]. 
The NNLO QCD effects are known in the large-mtop limit [53]. The 77 irreducible background has 
been computed up to NLO including the fragmentation effects [54]. The gg —>■ 77 contribution, which is 
formally NNLO, is enhanced by the large gluon luminosity and is known up to N^LO (i.e. 0{a^)) [55]. 

For Higgs masses between 140 and 180 GeV the W^W~ l'^l~vT' decay mode is the most 
important. Despite the absence of a mass peak, there are strong angular correlations between the charged 
leptons [56]. To suppress the background, a jet veto has to be applied to cut events with high-p^ 6-jets 
from the decay of the top quark. The impact of higher-order corrections on the Higgs signal is strongly 
reduced by the selection cuts [37,39,57]. This channel appears to be one of the most promising for an 
early discovery [58], but at the same time it is the most challenging as far as the background is concerned. 
Because of the missing energy, the Higgs mass cannot be directly reconstructed, and a straightforward 
background extrapolation from sidebands is not possible. The background has to be extrapolated from 
regions where the signal is absent and this requires a precise knowledge of the background distributions. 
The W~^W~ irreducible background is known up to NLO [59, 60] including spin correlations, and the 
effects of multiple soft-gluon emissions has been included up to NLL [61]. Spin correlations in the 
W decay are crucial for a correct prediction of angular distributions and are now implemented in the 
MC@NLO event generator [62,63]. The potentially large gg W'^W" contribution, formally NNLO, 
has also been computed [64,65]. The tt background, including the effect of spin correlations [66], is 
known up to NLO and is also included in MC@NLO. A complete calculation including finite width 
effects (and thus W+W-bb, Wtb) is available at LO only [67]. 

When the Higgs mass is larger than about 180 GeV, the decay h ZZ 4 leptons becomes 
dominant. This channel is much easier to observe than the W~^W~ channel because the invariant mass 
of the leptons can be reconstructed and thus the background can be measured from the data. Accurate 
predictions become important when the nature of the Higgs particles is studied. The irreducible ZZ 
background is known up to NLO including spin correlations [59,60]. The impact of soft-gluon effects 
on signal and background has been studied recently [68]. The calculation of the gg — > ZZ contribution 
is accounted for in this report in Sects. [6] and |T] We finally note that the full QCD-i-EW corrections to 
the decay modes h W^W^{ZZ) 4 leptons have been recently computed [69,70]. 

1.2 Vector-Boson Fusion 

The vector boson fusion (VBF) process plays an important role for a wide range of Higgs masses. The 
VBF cross section is typically one order of magnitude smaller than the one from gluon fusion, and it 
becomes competitive with the latter for very large Higgs masses. 

VBF occurs through the scattering of two valence quarks that exchange a.W or a. Z boson. Since 
valence quark distributions in the proton are peaked at relatively large Bjorken x (x ~ 0.1 — 0.2), the 
scattered quarks emerge with very large longitudinal momentum and transverse momentum of the order 
of a fraction of the boson mass. As a consequence, the typical signature of a VBF event is given by two 
hard jets with a large rapidity interval between them, and since the exchanged boson is colourless, there 
is no hadronic activity between them. Although this channel has a smaller cross section with respect to 
gluon fusion, it is very attractive both for discovery and for the measurement of the Higgs couplings. 

The NLO QCD corrections to the total rate were computed some time ago and found to be of the 
order of 5 — 10% [7 1 ]. In recent years, these corrections have been implemented for distributions [72-74] . 
Recently, the full EW-i-QCD corrections to this process have been computed [75,76]. A comparison of 
the different calculations is presented in this report in Sect. [3] 

The h+2 jets final state can be produced also by gluon-gluon fusion. This signature, although 



part of the inclusive Higgs boson signal, represents a background when trying to isolate the hWW 
and hZZ couplings through VBF. The gluon fusion contribution is known at LO with full top mass 
dependence [77]. The kinematical distributions of the tagging jets show remarkable differences in the 
two production mechanisms. The A(p distribution of the tagging jets is rather flat for the VBF signal. 
By contrast, the loop induced hgg coupling leads to a pronounced dip at Acj) = 90". Another significant 
difference is found in the rapidity distribution of the third hardest jet with respect to the rapidity average 
of the other two. The VBF signal has a dip in the central region, where the gluon fusion background is 
peaked. As such, a cut on jets with px > p^^° in the central rapidity region (central jet veto) enhances 
the relevance of the VBF signal. Recently, NLO QCD corrections to the /i + 2 jets process in the 
large mtop-limit have been computed [32], and also parton shower effects on the relevant distributions 
have been evaluated [78]. These studies show that the discriminating power of previous LO results 
is not significantly changed. We note, however, that when the p™*° is much smaller than the Higgs 
boson mass the coefficients of the perturbative series are enhanced by large logarithmic contributions 
that may invalidate the fixed order expansion. The latter point deserves more detailed investigation. An 
experimental study of central jet veto efficiencies is presented in this report in Sect. [S] 

The most important decay channels of the Higgs boson in VBF are /i ^ r+r" and h 
W^W^ —>■ l^l^vV. The h —>■ t~^t~ decay mode provides an important discovery channel in the 
MSSM. The r+r" invariant mass can be reconstructed at the LHC with an accuracy of a few GeV This 
is possible because VBF typically produces Higgs bosons with large transverse momentum. As a conse- 
quence, a sideband analysis can in principle be used to measure the background from the data. The most 
important backgrounds are QCD Zjj and EW Zjj from VBF. Both are known up to NLO [79, 80]. 

The h —>■ W^W^ —>■ l^l^vV decay mode is the most challenging, because, as for gluon fusion, 
it does not allow a direct Higgs mass reconstruction. The irreducible W^W^ background is known up 
to NLO [81]. The other important background is ti + jets, and has the largest uncertainty. Recently, the 
NLO corrections to tt+ jet have been computed [82]. It will be essential to include the decay of the top 
quark with full spin correlations. In addition, finite width effects could be relevant. 

1.3 Associated Production With a bh Pair 

In the Standard Model, Higgs production in association with b quarks is never important, since this rate 
is suppressed by mb/v. This channel is important in MSSM scenarios at large tan/3, since the Higgs 
coupling to bottom quarks is enhanced in this regime. For tan /? ^ 7, Higgs production in association 
with a b quark is the dominant production mechanism at the LHC. The cross section for b- Higgs produc- 
tion can be computed using two different formalisms, which represent different orderings of perturbation 
theory. In the four-flavour scheme the cross section starts with gg — > bbh at LO. The cross section for 
the associated production of the Higgs boson with zero, one or two high-transverse momentum 6-jets 
is known up NLO [83-86]. In the five flavour scheme, the LO process is 66 ^ /i and bottom quark 
parton distributions are introduced to sum the potentially large logarithms, \og{mh/mb). The inclusive 
cross section has been computed up to NNLO [87], and the cross section for the associated production 
with one high-pj- 6 jet is known at NLO [88]. In recent years, a detailed comparison between the results 
of the two approaches has been performed with the conclusion that the two approaches lead to similar 
results. For a discussion see Ref. [89,90]. In addition, the electroweak corrections to 66 ^ /i [91], the 
dominant top quark contributions to the NNLO rate for the exclusive bbh process [92] and the SUSY 
QCD corrections to 66 h,bg ^ bh are known [93]. The effects of SUSY-QCD on 6-Higgs production 
is discussed in Sect. 112.1 

1.4 Associated Production With a tt Pair 

The htt channel offers the possibility of a clean measurement of the top quark Yukawa coupling. The 
NLO corrections to the htt signal were independently computed by two groups [94-97], and found to 
increase the signal cross section by 20 — 40%. The /itt channel was initially thought to be an important 
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Fig. 1: Total cross sections for Higgs production at the LHC. The gluon fusion result is NNLO QCD with soft gluon resumma- 
tion effects included at NNLL and uses MRST2002 PDFs with renormalization/factorization scales equal to mh- The vector 
boson fusion curve is shown at NLO QCD with CTEQ6M PDFs and renormalization/factorization scales equal to m/i. The Vh 
results (V = W, Z) include NNLO QCD corrections and NLO EW corrections and use MRST2002 PDFs with the renormal- 
ization /factorization scales equal to the — Mv invariant mass. The bb ^ h result is NNLO QCD, with MRST2002 PDFs, 
renormalization scale equal to rrih and factorization scale equal to mh/4. The results for tth production are NLO QCD, use 
CTEQ6M PDFs and set the renormalization/factorization scale to mt + mh/2 [100]. 

discovery channel in the low Higgs mass region, looking at the h ^ bb decay mode and triggering on 
the leptonic decay of one of the top. The main backgrounds are tibb and tibjj. Recently, more detailed 
investigations based on a more careful background evaluation and full detector simulation lead to a more 
pessimistic view on the possibility of observing the Higgs signal in this channel [98]. This channel could 
be important for measuring the top quark Yukawa coupling [8,99]. 

1.5 Associated production with aW or a Z boson 

This channel is essential for the Higgs search at the Tevatron for Higgs masses below 130 GeV. The 
leptonic decay of the vector boson provides the necessary background rejection to allow for looking at 
the h ^ bb decay mode. The signal cross section is known up to NNLO in QCD, the corrections being 
about +30% [101, 102]. These corrections are identical to those of Drell-Yan, but in the case of Zh an 
additional contribution from the gg initial state must be included [103]. Full EW corrections are known 
and decrease the cross section by 5 — 10% [104]. 

1.6 Conclusions 

The important Higgs production channels are known at NLO QCD and in a few cases to NNLO and 
progress is being made in implementing these results in Monte Carlo programs. A summary of the total 
rates for the most important Higgs production channels is shown in Fig. [T][100]. 



Part II 

STANDARD MODEL HIGGS BOSONS 

2. HNNLO: A MONTE CARLO PROGRAM FOR HIGGS BOSON PRODUCTION AT THE 
LHCI 

2.1 Introduction 

Gluon-gluon fusion is the main production channel of the Standard Model Higgs boson at the LHC. At 
leading order (LO) in QCD perturbation theory, the cross section is proportional to Og, as being the 
QCD coupling. The QCD radiative corrections to the total cross section are known at the next-to-leading 
order (NLO) [9-11] and at the next-to-next-to-leading order (NNLO) [12-17]. The effects of a jet veto 
on the total cross section has been studied up to NNLO [35]. We recall that all the results at NNLO have 
been obtained by using the large- Af^ approximation, Mt being the mass of the top quark. 

These NNLO calculations are certainly important but they refer to situations where the experi- 
mental cuts are either ignored (as in the case of the total cross section) or taken into account only in 
simplified cases (as in the case of the jet vetoed cross section). The impact of higher-order corrections 
may be strongly dependent on the details of the applied cuts and also the shape of the distributions is 
typically affected by these details. 

The first NNLO calculation that fully takes into account experimental cuts was reported in 
Ref. [36], in the case of the decay mode H 77. In Ref. [37] the calculation is extended to the 
decay mode H —>■ WW — > Ivlv. 

In Ref. [38] we have presented an independent NNLO calculation of the Higgs production cross 
section. The method is completely different from that used in Refs. [36, 37]. Our calculation is imple- 
mented in a fully-exclusive parton level event generator. This feature makes it particularly suitable for 
practical applications to the computation of distributions in the form of bin histograms. Our numerical 
program can be downloaded from [105]. The decay modes that are currently implemented are if ^ 77, 
H WW lulu and H ^ ZZ ^4 leptons [39]. 

In the following we present a brief selection of results that can be obtained by our program. We 
consider Higgs boson production at the LHC and use the MRST2004 parton distributions [106], with 
parton densities and as evaluated at each corresponding order (i.e., we use [n + l)-loop as at N"LO, 
with n = 0, 1, 2). The renormalization and factorization scales are fixed to the value fJ-R = = Mh, 
where Mh is the mass of the Higgs boson. 

2.2 Results For the Decay Mode H ^ -yy 

We consider the production of a Higgs boson of mass Mh = 125 GeV in the ii ^ 77 decay mode 
and follow Ref. [98] to apply cuts on the photons. For each event, we classify the photon transverse 
momenta according to their minimum and maximum value, pTmm and pTmax- The photons are required 
to be in the central rapidity region, \r]\ < 2.5, with pTmin > 35 GeV and PTmax > 40 GeV. We also 
require the photons to be isolated: the hadronic (partonic) transverse energy in a cone of radius R = 0.3 
along the photon direction has to be smaller than 6 GeV. By applying these cuts the impact of the NNLO 
corrections on the NLO total cross section is reduced from 19% to 11%. 

In Fig. |2]we plot the distributions in pTmm and pTmax of the signal process gg H 77. We 
note that the shape of these distributions sizeably differs when going from LO to NLO and to NNLO. 
The origin of these perturbative instabilities is well known [107]. Since the LO spectra are kinematically 
bounded by px < Mh/2, each higher-order perturbative contribution produces (integrable) logarithmic 
singularities in the vicinity of that boundary. More detailed studies are necessary to assess the theoretical 
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Fig. 2: Distributions in pTmin and pTmax for the diplioton signal at tiie LHC. The cross section is divided by the branching 
ratio in two photons. 
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Fig. 3: Normalized distribution in the variable cos 9* 



uncertainties of these fixed-order results and the relevance of all-order resummed calculations. 

In Fig. |3]we consider the (normalized) distribution in the variable cos 9*, where 9* is the polar 
angle of one of the photons in the rest frame of the Higgs boson H At small values of cos 9* the distribu- 
tion is quite stable with respect to higher order QCD corrections. We also note that the LO distribution 
vanishes beyond the value cos < 1- The upper bound cos due to the fact that the photons are 

required to have a minimum pT of 35 GeV. As in the case of Fig. |2l in the vicinity of this LO kinematical 
boundary there is an instability of the perturbative results beyond LO. 



2.3 Results for the Decay Mode H ^ Ivlv 

We now consider the production of a Higgs boson with mass Mh = 165 GeV in the decay mode 
H —>■ Ivlv. We apply a set of preselection cuts taken from the study of Ref. [58]. The charged leptons 
have pt larger than 20 GeV, and |r/| < 2. The missing pT is larger than 20 GeV and the invariant mass of 
the charged leptons is smaller than 80 GeV. Finally, the azimuthal separation of the charged leptons in the 
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Fig. 4: Normalized A(f> distribution at LO, NLO, NNLO. 



transverse plane (Acp) is smaller than 135°. By applying these cuts the impact of the NNLO corrections 
on the NLO result does not change and is of about 20%. 

In Figj4]we plot the A(p distribution at LO, NLO and NNLO. As is well known [56], the charged 
leptons from the Higgs boson signal tend to be close in angle, and thus the distribution is peaked at small 
A(j). We notice that the effect of the QCD corrections is to increase the steepness of the distribution, from 
LO to NLO and from NLO to NNLO. 

2.4 Conclusions 

We have illustrated a calculation of the Higgs boson production cross section at the LHC up to NNLO 
in QCD perturbation theory. The calculation is implemented in the numerical program HNNLO, which 
at present includes the decay modes — > 77 and H WW — > Ivlv and H ZZ — > 4 leptons. 
The program allows the user to apply arbitrary cuts on the momenta of the partons and leptons (photons) 
produced in the final state, and to obtain the required distributions in the form of bin histograms. We 
have presented a brief selection of numerical results that can be obtained by our program. More detailed 
results for the decay modes H — > WW and H ZZ can be found in Ref. [39]. The fortran code 
HNNLO can be downloaded from [105]. 

3. TUNED COMPARISON OF QCD CORRECTIONS TO SM HIGGS-BOSON PRODUCTION 
VIA VECTOR BOSON FUSION AT THE LHC0 

3.1 Introduction 

The electroweak (EW) production of a Standard Model Higgs boson in association with two hard jets in 
the forward and backward regions of the detector — frequently quoted as the "vector-boson fusion" (VBF) 
channel — is a cornerstone in the Higgs search both in the ATLAS [108] and CMS [109] experiments at 
the LHC. Higgs production in the VBF channel also plays an important role in the determination of 
Higgs couplings at the LHC (see e.g. Ref. [8]). Even bounds on non-standard couplings between Higgs 
and EW gauge bosons can be imposed from precision studies in this channel [110]. 

Higgs-i-2jets production in pp collisions proceeds through two different channels. The first channel 
corresponds to a pure EW process where the Higgs boson couples to a weak boson. It comprises the 
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scattering of two (anti-)quarks mediated by t- and u-channel W- or Z-boson exchange, with the Higgs 
boson radiated off the weak-boson propagator. It also involves Higgs-boson radiation off a W- or Z- 
boson produced in s-channel quark-antiquark annihilation (Higgs-strahlung process), with the weak 
boson decaying hadronically. The second channel proceeds mainly through strong interactions, the Higgs 
boson being radiated off a heavy-quark loop that couples to any parton of the incoming hadrons via 
gluons [32,77]. 

In the weak-boson-mediated processes, the two scattered quarks are usually visible as two hard 
forward jets, in contrast to other jet production mechanisms, offering a good background suppression 
(transverse-momentum and rapidity cuts on jets, jet rapidity gap, central-jet veto, etc.). Applying appro- 
priate event selection criteria (see e.g. Refs. [78, 1 1 1-1 14] and references in Refs. [2, 115]) it is possible 
to sufficiently suppress background and to enhance the VBF channel over the hadronic Higgs-i-2jets pro- 
duction mechanism. 

In order to match the required precision for theoretical predictions at the LHC, QCD and EW 
corrections are needed. When VBF cuts are imposed, the cross section can be approximated by the 
contribution of squared t- and u-channel diagrams only, which reduces the QCD corrections to vertex 
corrections to the weak-boson-quark coupling. Explicit next-to-leading-order (NLO) QCD calculations 
in this approximation exist since more than a decade [71, 115], while corrections to distributions have 
been calculated in the last few years [72-74]. Recently, the full NLO EW and QCD corrections to this 
process have become available [75,76]. This calculation includes, for the first time, the complete set of 
EW and QCD diagrams, namely the t-, u-, and s-channel contributions, as well as all interferences. 

In this short note we compare the NLO QCD corrected cross-section results obtained by three 
different calculations using a common set of input parameters and a uniformly tuned setup. We also 
present, in order to better understand the different approximations, the full NLO QCD and EW corrected 
results as obtained in Refs. [75,76]. 

In the next section, the different approaches that we compare are briefly summarized. The precise 
setup is described in Section 1331 and Section [3]4] contains the numerical results. 

3.2 Different Approaches and Codes 

The following collaborations have contributed to the tuned comparison of NLO QCD conected results 
for Higgs-boson production via weak-boson fusion at the LHC: 

• CDD: References [75,76] present a detailed description of the calculation of the complete NLO 
EW and QCD corrections to Higgs-boson production in the VBF channel at the LHC. The NLO 
0{as) corrections include the complete set of QCD diagrams, namely the t-, u-, and s-channel 
contributions, as well as all interferences. The integrated cross section (with and without dedicated 
VBF selection cuts) was calculated, as well as different Higgs-boson and tagging-jet observables. 
In the EW corrections, real corrections induced by photons in the initial state and QED corrections 
implicitly contained in the DGLAP evolution of PDFs were also taken into account. All EW 
contributions have been switched off for this comparison. 

• VBFNLO [116] is a NLO parton-level Monte Carlo program which implements one-loop QCD 
corrections for a collection of relevant VBF processes, of which Higgs-boson production, in the 
narrow resonance approximation, is the simplest example. Higgs-boson production in weak-boson 
fusion is implemented following the results of Ref. [72]. VBFNLO generates an isotropic Higgs- 
boson decay into two massless "leptons" (which represent r+r" or 77 or hb final states), and 
imposes a cut on the invariant mass of the Higgs boson. This feature has been disabled during 
this comparison, and only a non-decaying Higgs boson has been considered. We have employed 
VBFNLO- V .1.0, and included only four flavours of the external quarks. 

• W2H [117] calculates the production cross section of Higgs bosons via WW/ZZ h, H at 
hadron colliders at NLO QCD according to the formulae presented in Refs. [71, 1 15]. Interference 



effects between W- and Z-boson fusion are neglected. The program allows to calculate the total 
production cross section for the scalar Higgs bosons of the SM and MSSM. For the present study 
we employed the W2H version dated July 23, 2007, which was modified in order to switch off the 
contributions from b quarks in the final and/or initial states. 

3.3 Common Setup for the Calculation 

3.31 Input parameters and scheme definitions 

We choose the following set of input parameters [118], which have also been used in Refs. [75,76]: 
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CDD uses the complex-mass scheme [119]. This requires a fixed width in the W- and Z-boson 
propagators in contrast to the approach used at LEP to fit the W and Z resonances, where running widths 
are taken. Following Ref. [120] to convert the "on-shell" values of M^^^ and T\f'^ {V = W, Z) to the 
"pole values" denoted by My and Fy, leads to 

Mw = 80.397 ... GeV, Fw = 2.123 ... GeV, 

Mz = 91.1535 . . . GeV, Fz = 2.4943 . . . GeV. (2) 

In VV2H and VBFNLO the W- and Z-boson masses are fixed according to Eq. Q and the vector-bosons 
widths are set to zero. 

The masses of the light quarks are adjusted to reproduce the hadronic contribution to the photonic 
vacuum polarization of Ref. [121]. Since quark mixing effects are suppressed we neglect quark mixing 
and use a unit CKM matrix. All quark masses are set to zero in VBFNLO. We use the scheme, i.e. we 
derive the electromagnetic coupling constant from the Fermi constant according to 

ac, = V2G^M^{1 - M^/Mi)/7T. (3) 

CTEQ6 parton distributions [122] are used. Processes with external bottom quarks are not in- 
cluded in this comparison. As discussed in Section 3.4 of Ref. [76] these contribute at the level of a 
few per cent. We use Afw as factorization scale both for QCD and QED collinear contributions. For 
the calculation of the strong coupling constant we employ Mw as the default renormalization scale, in- 
clude 5 flavours in the two-loop running of as, and fix as{Mz) = 0.118, consistent with the CTEQ6M 
distribution. 



3.32 Phase-space cuts and event selection 

We employ the same jet definition parameters, phase-space and event selection cuts as described in 
Refs. [73,75,76]. Jet reconstruction from final-state partons is performed using the /cx-algorithm [123] 
as described in Ref. [124]. Jets are reconstructed from partons of pseudorapidity \r]\ < 5 using a jet 
resolution parameter D = 0.8. In the EW corrections, real photons are recombined with jets according 
to the same algorithm. Thus, in real photon radiation events, final states may consist of jets plus a real 
identifiable photon, or of jets only. 

We study total cross sections and cross sections for the set of experimental "VBF cuts". These 
cuts significantly suppress backgrounds to VBF processes, enhancing the signal-to-background ratio. 
We require at least two hard jets with 

PTj>20GeV, |yj|<4.5, (4) 
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Table 1: Total integrated cross section for pp H + 2jets + X in LO and NLO without any cuts, calculated by CDD, ctJ 
with VV2H, crJ^^/NLO' and with VBFNLO, ctlq/nIo' for the setup defined in the text. 



where pxj is the transverse momentum of the jet and yj its rapidity. Two tagging jets ji and j2 are defined 
as the two jets passing the cuts ^ with highest pT such that pTj^ > PTJa- Furthermore, we require that 
the tagging jets have a large rapidity separation and reside in opposite detector hemispheres: 

^yjj = |yji - yj2l > 4, yji • y^^ < o. (5) 



3.4 Numerical Results 

In this section we present, for a range of Higgs-boson masses, LO and NLO QCD corrected results 
obtained by CDD, o-™o/nlO' VV2H, o-j^^^/NLO' VBFNLO, o-lq/^lo- These results were 

calculated approximating the cross section by the contribution of squared t- and u-channel diagrams 
only, without any interferences. We also present the QCD corrected results, including all diagrams and 
interference contributions, o"l^'/nlo' together with the results including both QCD and EW corrections, 

'^LO/NLO^' obtained by CDD. 

Table[T]contains results for the total integrated cross section without any cuts. The small difference 
between the results obtained by VV2H and VBFNLO is due to the different treatment of vector-boson 
widths. We observe that the approximate LO cross sections agree within 5 x 10^^, and the NLO corrected 
results within 2 x 10^^, a difference which is of the order of the statistical error. The complete predictions 
^QCD+EW differ from the results of VV2H and VBFNLO by up to 30% for low Higgs-boson masses and 
by a few per cent for high Higgs-boson masses. The bulk of this big difference for small values of Mh 
is due to the s-channel contributions, which are only considered by CDD. 

Table |2] shows results for the integrated cross section after imposing VBF selection cuts. We 
observe that the approximate LO cross sections agree within 8 x 10~^, and the NLO corrected results 
within 1 X 10~^, a difference which is of the order of the statistical error. The difference between the 
complete predictions q-QCD+ew j-j^g results of VBFNLO is half a per mille or less in LO, and 6-8%, 
the size of the EW corrections, in NLO. This shows that, in this configuration, s-channel and interference 
contributions can be safely neglected, but EW corrections are as large as QCD corrections. 



3.5 Conclusions 

We have presented results for NLO cross sections of Standard Model Higgs-boson production via weak- 
boson fusion at the LHC. A tuned comparison of QCD corrected results obtained by three different 
calculations has been performed. Taking into account only t- and u-channel diagrams we found good 
agreement. We have also presented full NLO EW and QCD corrected results to gain some insight 
into the nature of this approximation. We found agreement between the approximate and full 0{as) 
results when VBF cuts are applied. On the other hand, for the total integrated cross section, there is a 
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Table 2: Integrated cross section for pp ^ H + 2jcts + X in LO and NLO, including VBF selection cuts, calculated by CDD, 
CTlo/nlO' ™'1 VBFNLO, o-Lo/NLO' f*"" *e setup defined in the text. 



sizeable difference between those results, which arises almost exclusively from s-channel contributions. 
Furthermore, EW corrections are in general as large as QCD corrections. 
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4. LOOP INDUCED INTERFERENCE EFFECTS IN HIGGS PLUS TWO JET PRODUCTION 
AT THE LHC 

4.1 Introduction 

Understanding the mechanism of electro-weak symmetry breaking is one of the primary goals at the 
CERN Large Hadron Collider. Central to this study is the measurement of the couplings of any 
observed Higgs scalar to the electro-weak bosons. A useful production process in this context is 
pp Hjj [125-127] through weak boson fusion (WBF) [128], as shown in Fig.|5la), with contributions 
from all identifiable decay channels. The Higgs plus two jet signature also receives contributions from 
Higgs boson production through gluon-fusion mediated through a top-loop, as illustrated in Fig. |5](b). 
However, the Higgs plus dijet-sample can be biased towards WBF by suppressing the gluon-fusion chan- 
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Fig. 5: (a) The WBF process for Higgs production in the Standard Model and (b) the equivalent gluon-fusion diagram mediated 
through a top-loop. 

nel through a combination of cuts. 

The next-to-leading order corrections to Higgs plus two jet production are considered to be well 
under control. For WBF, both the radiative corrections within QCD [10, 71, 72, 129] and the electro- 
weak sector [75, 76] have been calculated; for the gluon fusion process, the first radiative corrections 
have been calculated within QCD [32, 130] using the heavy top mass effective Lagrangian [9, 10, 131]. 
The radiative corrections to the WBF channel are small, 3% — 6%, and there is even partial numerical 
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(a) (b) (c) 

Fig. 6: Example of contributing one-loop interference terms: (a) A4gzMlg and (b) MggM*z- There are four contributing 
topologies for each gluon-fusion and Z-fusion process, (c) shows a real emission processes at matrix element squared level. 

cancellation between the QCD and electro-weak contributions. It would therefore seem that the Higgs 
coupling to electro-weak bosons can be very cleanly studied with a i/jj-sample. 

However, there is an irreducible (i.e. unaffected by the WBF cuts) contamination in the extraction 
of the ZZiJ-coupling from interference between the gluon fusion and WBF processes, which was ig- 
nored in the literature until recently. At tree level, such interference is only allowed in amplitudes where 
there is a t •s-^ n-channel crossing which leads to a high level of kinematic suppression [132]. These 
and other crossing-suppressed one-loop amplitudes were later calculated together with the electro-weak 
corrections [75,76]. 

Here we will report on the calculation of the processes allowed at the one-loop level which do 
not suffer from the kinematic suppression stemming from the requirement of a t ^ u-crossing [133]. 
At order 0{a^a^g), one finds an interference term between the gluon- and Z-induced amplitude which 
is not allowed at 0{a^a1) by colour conservation. The VF-induced amplitudes are crossing-suppressed 
and therefore not taken into account. The diagrams where the vector boson is in the s-channel can also 
be safely neglected because they are strongly suppressed by the WBF cuts. As discussed in Ref. [132], 
for identical quark flavours the loop amphtudes are the first order which does not require a kinematically 
disfavoured crossing. 

In the following section we will briefly sketch the calculation before discussing our results in 
section 1431 which are summarized in the conclusions. 

4.2 The Calculation 

Our calculation of the loop interference terms and the real emission contributions is based on helicity 
amplitudes. The leading order amplitudes, denoted hy M.z and M.g (Fig. |5la) and (b)), are proportional 
to a colour singlet and a colour octet term. The colour singlet is formally of order 0{a^) whereas 
the octet is of order 0{a1). The one-loop amplitudes, which we call M.gz and M.gg respectively, are 
mixtures of octet and singlet terms. For the interference term we need to consider only the octet part 
of Mgz and the singlet part of M.gg. One finds that only four one-loop five-point topologies for each 
amplitude survive this colour projection. Sample diagrams are shown in Fig.|6l 

The loop amplitudes require the evaluation of one-loop five-point tensor integrals with a mixture 
of massless and massive configurations in both propagators and external legs. We apply the reduction 
algorithm outlined in Ref. [134, 135] to express each Feynman diagram as a linear combination of 1-, 2-, 
and 3-point functions in D = 4 — 2e dimensions and 4-point functions in D=6. The same algorithm has 
been successfully applied to a number of one-loop computations [64, 136-138], where further details can 
be found. The algebraic expressions were checked by independent implementations, both amongst the 
authors and with another group [139]. 

After the algebraic reduction, all helicity amplitudes are obtained as linear combinations of scalar 
integrals. No one-point functions appear in the reduction, and also two-point functions are absent in the 
amplitudes of M.gz- Furthermore, coefficients of some of the integrals which arise in several topologies 
sum to zero. If the tree resulting from a certain cut of a master integral corresponds to a helicity forbidden 
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Table 3: The cuts used in the following analysis which bias the Higgs Boson plus dijet sample towards WBF. The indices a, b 
label the tagged jets. 

tree level process, one can immediately infer the vanishing of the corresponding coefficient. In our 
algebraic tensor reduction approach we verify and use such cancellations before the numerical evaluation 
of the cross section. 

As most of the required scalar integrals are not provided in the literature, we have evaluated rep- 
resentations in terms of analytic functions valid in all kinematic regions. These can be found in [133] for 
use in other calculations. 

We used dimensional regularisation to extract the IR singularities from the divergent integrals. The 
leading poles cancel, but there remains al/e pole which is cancelled when the virtual corrections 
are combined with the real emission part shown in Fig. Oc). As to be expected, the coUinear IR diver- 
gences from the three-parton final states integrate to zero, leaving only a soft divergence proportional to 
1/e, which we isolated using the phase space slicing method [140, 141]. The phase space integration and 
the numerical evaluation of integrals and coefficients is coded in a C++ program allowing for a flexible 
implementation of cuts and observables. 

4.3 Results 

As the aim of our study was to investigate a possible pollution of the clean extraction of the ZZH vertex 
structure by the interference terms, we apply the cuts summarised in Table |3] These are generally used 
to single out the WBF events from the gluon fusion "background" [77]. Our input parameters for the 
numerical studies are taken from [106] and [142]. In addition, we use a Higgs boson mass of 115 GeV 
and the NLO parton distribution set from Ref. [106]. We use 2-loop running for q^, in accordance with 
the chosen pdfs. 

We observe that in all the flavour and helicity channels, the contribution from the 3-parton final 
state is numerically negligible. The only role of this real emission is to cancel the divergences which 
arise from the one-loop diagrams. 

As the interference effect is proportional to 2Re{AiggAi*^ + M-gzJ^g), the result is not neces- 
sarily positive definite. In fact, the sign of the interference contribution depends on the azimuthal angle 
between the two tagging jets, As the event topology has two well separated jets, it becomes pos- 

sible to define an orientation of the azimuthal angle which allows observability in the whole range of 
[—IT, vr], as pioneered in Ref. [1 10, 128]. ^(pjj is then defined through 

\P+T I \P-T I cos Acpjj = p+j, • p„ J, , 
2\p+J\p^JsmA(j)jj = e^upab+p+b^^pZ., 

where 6+ (6_) are unit vectors in positive (negative) beam direction, and likewise for the jet momenta 
p±. The cuts ensure that the two tagging jets lie in opposite hemispheres. Defined in this way, the observ- 
able A.(j)jj becomes a powerful discriminator for different CP-structures of the Higgs Boson production 
vertex [1 10]. 

Figure |7] displays the contribution to the distribution in ^.(pjj from the interference terms for var- 
ious helicity and flavour configurations. There is an accidental cancellation of sea and valence quark 
contributions which leads to the fact that the sum over all flavour and heUcity assignments peaks at 
around 2 ab/rad only, with an integrated effect of 1.19 ± 0.07 ab, where the error is due to the numerical 
integration. 




Fig. 7: The A<^jj -distribution for various flavour and fielicity-configurations. The purple histogram labelled "Sum" indicates 
the sum over the four contributions shown. The sum over all flavour and helicity assignments including all sea flavours is shown 
in the black histogram. 



Due to the oscillatory behaviour, the total integrated cross section does not represent the impact 
on the A(f)jj distribution. The integral of the absolute value of the A(f>jj distribution, J^^ d,^4>jj dA^ ' 
is a more useful measure of the impact of the interference effect on the extraction of the ZZif-vertex. 
This integral evaluates to 9.1 ± 0.1 ab, an order of magnitude larger than the integral over the oscillating 
distribution. The total integral over the the absolute value of the fully differential cross section leads to 
29.59 ± 0.07 ab. 

Using the same cuts and value for the mass of the Higgs boson as in the present study, we have 
checked that the total contribution to the A(/)jj -distribution from the leading order WBF process (both 
Z and W^l~ included) is relatively flat at around 240 fb/rad. Therefore, the result of the interference 
effect reported here is unlikely to be measurable. 

As can be readily seen in Fig. |71 there is also a cancellation between the contribution from each 
flavour and helicity assignment, as has also been pointed out in [143]; this is because the sign of quark 
couplings to the Z-boson becomes relevant as it is not squared for the interference. The flavour- and 
helicity sum for each quark hne therefore leads to some cancellation, which amounts to roughly 20% in 
the most relevant regions of the pdfs [133]. 

The complex phases arising from the full one-loop calculation of the amplitudes also give rise to 
some suppression. We find that the relevant products and sums for the interference effect project out only 
about 20% of the full complex loop amplitudes. 

We chose the factorisation and renormalisation scales as in accordance with the natural scales in 
the relevant high energy limit (as explained in Ref. [78]), i.e. the factorisation scales are set equal to the 
transverse momenta of the relevant jet, and the renormalisation scale for the strong couplings are chosen 
correspondingly, i.e. one evaluated at each value of the transverse momentum of the jets, and one at 
the Higgs mass. However, we find that varying the choice of factorisation and renormalisation scales, 
the exact numerical values of the cuts or the parameters, or the choice of pdf set has no impact on the 
conclusions: the numerical importance of the interference is basically unchanged. 



4,4 Conclusions 

We have outlined the calculation of the loop-induced 0{o?o?g) interference effect between the gluon 
fusion and weak boson fusion processes in Higgs boson plus two jet production at the LHC. 

We find by explicit calculation that this contribution is too small to contaminate the extraction 
of the ZZi/-coupling from WBF processes. Interestingly the effect which survives comes dominantly 
from the virtual corrections. We have analysed in detail why this contribution is so small, and instead of 
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Table 4: The cuts used in the following analysis which bias the Higgs boson plus dijet sample towards WBF. The suffices c, d 
label the tagged jets, j any (possibly further) jet in the event. 

a single effect we rather find a conspiracy of several mechanisms: 

• accidental cancellations between the sea quark and valence quark contributions 

• compensations between different weak isospin components of the valence quarks due to their 
SU{2) X U{1) couplings, in combination with their weights from the (valence) quark content of 
the proton 

• cancellations due to destructive interference of the phases from the different contributions. 

The exact impact of these partly accidental effects, in particular the latter, could not be assessed without 
an explicit calculation. 

Finally we would like to point out that anomalous couplings which affect the phases could change 
the interference pattern substantially. However, the first two cancellation mechanisms still being present, 
we still expect the overall contribution to be experimentally insignificant. Please see Ref. [144] for more 
details. 
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5. HIGGS BOSON PRODUCTION IN ASSOCIATION WITH MULTIPLE HARD JETS E 
5.1 Introduction 

It is widely hoped that the LHC will discover the source of electro-weak symmetry breaking, mediated by 
the Higgs scalar within the context of the Standard Model. In order to determine whether any observed 
fundamental scalar is the Higgs Boson of the Standard Model, it is imperative to determine its couplings, 
especially to the weak gauge bosons. This is possible both by measuring the decay of the Higgs boson 
through the weak bosons, but also by isolating the Higgs Boson production through weak boson fusion 
(WBF). This process contributes to the signal for the production of a Higgs boson in association with 
two jets. This channel also receives a significant contribution from higher order corrections to Higgs 
boson production through gluon fusion. In fact, it has recently been suggested [145] that the increased 
significance of the signal over the background obtained by requiring at least two hard jets in association 
with a Higgs boson may decrease the necessary integrated luminosity required for a discovery of the 
Higgs boson through gluon fusion processes. However, in order to measure the Higgs boson couplings 
to the weak bosons, it is necessary to suppress the gluon fusion contribution to the production of a 
Higgs boson in association with two jets. This is achieved [77] by applying the so-called weak boson 
fusion-cuts: It is expected that the contribution from the gluon fusion process will be further suppressed 
relative to WBF by vetoing further jet activity [146]. The efficiency of such cuts can only be assessed 
by calculating the higher order corrections to the gluon fusion contribution to the /ijj-channel. The first 
radiative corrections have recently been calculated [32]. While this fixed order approach certainly is the 
best tested and understood approach for predicting the first few perturbative corrections, the calculational 
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complexity means that currently the production of hjj through gluon fusion has only been calculated at 
next-to-leading order. 

It is possible to estimate final state jet emission in this process [78] using parton shower algorithms. 
In this contribution we examine a different approach, and consider how to best estimate hard jet emission 
in Higgs production via gluon fusion. We take as a starting point a factorised form for the scattering 
amplitudes, which formally applies in a certain kinematic limit (that of multi-Regge-kinematics (MRK)). 
We extend the domain of applicability of the amplitudes from Asymptotia into the region of relevance 
for the LHC by using known all-order constraints of scattering amplitudes. We validate the approach 
by checking the approximations in a comparison with fixed order results, where these are available. 
Furthermore, the resulting estimate for the n-parton final state (which includes some virtual corrections) 
is then matched to the known tree level results for hjj and hjjj. Finally, we implement the description in 
a Monte Carlo event generator for Higgs -i- multiparton production, and present a sample of results. 



5.2 Estimating Multijet Rates 

5.21 The FKL Amplitude 

Our starting point is the FKL factorised {2 —>■ n + 2)-gluon amplitudes [147] adapted to include also a 
Higgs boson 

^-^ 'HE ~ '^I'O yi-ysj y^iai-Lo 



W (^^exp[a(g2)(y.„i - y.)] {ig,f^'^^'^^+^^ C^^{qi,qi+i) 
(^exp[a{qf){yj - yh)]CHiqj+i,qh)j U) 

n , 



i=j+i 



^ exp[a{ql+M - yb)] 



2 ^^fl'^\Hn+lJ\yn yoji yysj y/MtlJ-n+i 
Qn+l 

where gs is the strong coupling constant, and qi, qh are the 4-momentum of gluon propagators (e.g. qi = 
Pa — X]fc=o Vk for i < j), C^jj is the Lipatov effective vertex for gluon emission, and Ch is the effective 
vertex for the production of a Higgs boson, as calculated in Ref. [148]. Furthermore, a{qf) occurs 
from the Reggeisation of the gluon propagator, and encodes virtual corrections (see e.g. [149]). This 
approximation formally applies in the MRK limit, which can be expressed in terms of the rapidities {y.i} 
of the outgoing partons and their transverse momenta {pi±}: 

yo > yi > • • • > yn+i; Pi±^Pi+i±; qf-q]- (8) 

This limit is particularly well suited for studies within the WBF cuts of Table HJ since a large rapidity 
span of the event is then guaranteed. 

In the true limit of MRK, the squared 4-momenta qf —>■ —q\^, and the square of the Lipatov 
vertices fulfil —C^^C^^ — > 4ii^ii+i_ Applying these limits, the sum over n to infinity of the amplitudes 
in Eq. dTjl, integrated over the full phase space of emitted gluons can be obtained by solving two coupled 
BFKL equations. This result would then apply to the phase space of 

yo > 2/1 > • • • > yn+i; m-Pi+i±; ^L-^ij- (9) 

While both expressions are equally valid in the region of Asymptotia, we extend the applicability 
of the results obtained in the High Energy Limit to the region of interest for particle physics phenomenol- 
ogy by adhering to the following guidelines: 
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Fig. 8: The 2 and 3 parton cross-sections calculated using the known LO matrix elements (solid), and the estimate gained from 
the modified high energy limit (dashed). One sees that the estimate is well within the scale variation of the LO result, obtained 
by varying the common renormalisation and factorisation scales in the range 0.5 < /i/^o < 2, where fio is the default choice 
(indicated by the shaded regions). Also shown is the result obtained from the fully inclusive n-parton sample of Eq. l|7j. 



1. Do NOT INTRODUCE NEW DIVERGENCES: Using the expression in Eq. ^ corresponds to re- 
moving some divergences from the full scattering amplitude (the collinear divergences), but not 
moving any divergences. The expression in Eq. (O is divergent only for momentum configura- 
tions for which the full scattering amplitude is also divergent. This is different to the case where 
the MRK limit of invariants has been substituted (resulting from the use of the BFKL equation), 
which displaces divergences within the phase space region of interest for the LHC. 

2. Do NOT APPLY THE FORMALISM WHERE IT FAILS: We choose minimal interception by only 
removing the small region on phase space where the expression of Eq. dTJl results in unphysical 
(negative) differential cross sections. This happens when the effective Lipatov vertex is applied to 
momentum configurations very far from the MRK, where it is possible to obtain —C^fi^^ < 0. 
It is perhaps interesting to note that restricting the region of phase space where the formalism is 
applied is similar to the kinematic constraint of Ref. [150-152], although in fact the latter fails to 
exclude all of the region where the formalism underpinning the BFKL equation fails. 

In figure [8] we compare the prediction for the production of a Higgs boson in association with 
two and three partons (in a hard two-jet and three-jet configuration respectively) within the WBF cuts of 
tablelH obtained using both the full matrix element (extracted from MADEvent/MADGraph [153]) and 
the relevant expression of Eq. ([7]l for two and three parton production, with the virtual corrections set to 
zero (d = 0). We choose renormalisation and factorisation scale in accordance with the study of Ref. 
[78]. One notes two things. Firstly, the approximation to the jet rates is well within the scale uncertainty 
of the known tree level results. We have therefore explicitly shown that the terms taken into account 
in this approach indeed dominate . Secondly, the cross section for the production of a Higgs boson in 
association with 3 jets is similar to the one for the production of a Higgs boson in association with two 
jets. The large size of the three-jet rate was already reported in Ref. [130], and clearly demonstrates the 
necessity of considering hard multi-parton final states in order to describe correctly the expected event 
topology and to answer questions on e.g. the effectiveness of a central jet veto in suppressing the gluon 
fusion channel. 




LO 2-jet 
f^-parton, = 2-jet 

LO 3-jet 
n-parton, > 2-jet 



0.504 ±0.0013 
0.267 ± 0.0034 
0.228 ±0.0018 
0.161 ±0.0087 



Table 5: The angular decorrelation parameter given by equation l llOt . subject to the cuts of table|4] Note that the 2 and 3-jet 
values are obtained from matrix elements matched to the known tree level results. 

5.3 All Order Results and Matching 

The divergence in Eq. ^ obtained when any ^ is regulated by the divergence of the virtual 
corrections encoded in a. By implementing the regularisation through phase space slicing it becomes 
possible to obtain the fully inclusive any-parton sample by summing Eq. Q over all j, n. This is very 
efficiently implemented by following the method for phase space generation outlined in Ref. [154]. 
Furthermore, since we can trivially expand the expressions to any order in Ug, it is possible to check the 
performance of the formalism against the available tree-level results, and to implement matching to these. 
We choose to implement In /^-matching at the amplitude-level for channels which have a contribution 
in the high-energy limit (e.g. ug hug and gg hggg), and i?-matching for those which do not 
(e.g. gg huu and uu hggg). 

It is now possible to cluster each event in the inclusive sample of a Higgs boson plus n partons 
into jets according to a given algorithm. As an example, we choose KtJet [155]. We use the parton 
distribution functions of Ref. [106]. The distribution of final state jets subject to the cuts of table |4] is 
shown with the dashed histogram in Figure [8] One sees a significant number of events with more than 3 
hard {p± > 40GeV) jets. More importantly though, the method outlined in this paper allows for an esti- 
mate of the emissions of partons not quite hard enough to be classified as jets, but still causing sufficient 
decorrelation. The azimuthal angular correlation between the tagging jets has been suggested previously 
as a good observable for differentiating between the GGF and WBF production modes. Furthermore, the 
nature of the distribution of the azimuthal angle between the two tagging jets can potentially be used to 
determine the nature of the Higgs coupling to fermions [145]. However, the usefulness of this observable 
is threatened by hard jet emission which acts to decorrelate the tagging jets. As suggested in Ref. [156] 
the structure of the distribution da /d(pj^j^^ can be distilled into a single number given by: 



The results using our approach are collected in Table[5] Of particular interest is the difference between the 
first two numbers. The first {A^ = 0.504 ± 0.0013) describes the result obtained in the two-jet tree-level 
calculation. The second (A^ = 0.267 ± 0.0034) is the result obtained for events classified as containing 
only two hard jets, but completely inclusive in the number of final state partons. The difference is mostly 
due to the decorrelation caused by the additional radiation not classified as hard jets. As expected, the 
further hard emissions have a stronger effect than estimated using a parton shower approach [78]. 

5.4 Conclusions 

We have outlined a new technique for estimating multiple hard parton emission, and demonstrated its 
application to Higgs boson production (via GGF) in association with two jets. Our starting point is the 
FKL factorised form of Higgs-i-multijet amplitudes, which formally applies in multi-Regge kinematics 
(MRK). We extend the region of applicability of the formalism by adhering to two simple rules. We 
compare the results obtained order by order to those obtained in a fixed order approach and find very good 
agreement. The approximations are well within the uncertainty estimated by varying the renormalisation 
and factorisation scale by a factor of two in the tree level results. 




^i^jgjt < ^/4) - <^(7r/4 < < 37r/4) + ajc^j^j^ > Stt/A) 
'^{Kib < ^/4) ± <^(vr/4 < < 37r/4) ± a{^j^j^ > Stt/A) 



(10) 



We have presented example results for the distribution of final state jets, and for the azimuthal 
decorrelation parameter A^. We find significant decorrelation arising from additional hard final state 
radiation not captured by present NLO calculations; significantly more than previously estimated using 
parton shower algorithms. 

The technique outlined here can be extended to e.g. VF+jet emission, as well as pure multijet final 
states. It would be very interesting to interface the final states found here with parton shower algorithms, 
thus resumming in principle both the number of jets (hard partons) and the structure of each (soft coUinear 
radiation). Furthermore, the results presented here are based upon effective vertices correct to leading 
logarithmic order. Work is in progress towards extending the accuracy to next-to-leading logarithmic 
order. 
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6. GLUON- INDUCED Z-BOSON PAIR PRODUCTION AT THE LHC: PARTON LEVEL RE- 
SULTS |!| 

6.1 Introduction 

The hadronic production of Z boson pairs provides an important background for Higgs boson searches in 
the H ZZ channel at the LHC. It has been studied extensively in the literature including higher order 
corrections [59, 60, 157, 158]. Production of Z boson pairs through gluon fusion contributes at 0{a'^) 
relative to qq annihilation, but its importance is enhanced by the large gluon iiux at the LHC. It was 
analyzed in Refs. [159, 160]. Leptonic Z decays were subsequently studied for on-shell [161] and off- 
shell [162] vector bosons. In this note we present the first complete calculation of the gluon-induced loop 
process gg —>■ ^(7)^(7) —>■ £££'i', allowing for arbitrary invariant masses of the Z bosons and including 
the 7 contributions. Our calculation employs the same methods as Refs. [64, 136]. The tensor reduction 
scheme of Refs. [134, 135] has been applied to obtain the amplitude representation implemented in our 
program. We compared it numerically with an amplitude representation based on FeynArts/FormCalc 
[163, 164] and found agreement. Note that single resonant diagrams (in the case of massless leptons) 
and the corresponding photon exchange diagrams give a vanishing contribution. A combination of the 
multi-channel [165] and phase-space-decomposition [67, 166] Monte Carlo integration techniques was 
used with appropriate mappings to compensate peaks in the amplitude. A more detailed description of 
our calculation can be found in a forthcoming article. 

6.2 Results 

In this section we present numerical results for the process pp ^(7)^(7) ill' I' at the LHC, i.e. for 
the production of two charged lepton pairs with different flavor. Note that no flavor summation is applied. 
First, we give the cross section when standard LHC cuts for Z boson production [60] are applied. More 
precisely, we require 75 GeV < < 105 GeV for the invariant masses of II and I'l' , which suppresses 
the photon contribution to less than 1%. Motivated by the finite acceptance and resolution of the detectors 
we further require pTi > 20 GeV and \ri(\ < 2.5 for all produced leptons. To obtain numerical results 
we use the following set of input parameters: A% = 80.419 GeV, Mz = 91.188 GeV, Gf, = 1.16639 x 
10"^ GeV "^, Tz = 2.44 GeV. The weak mixing angle is given by c„ = Mw/Mz, si = I - c^. 
The electromagnetic coupling is defined in the scheme as = \/2G^M^s^/7r. The masses of 
external fermions are neglected. The values of the heavy quark masses in the intermediate loop are set to 
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a{pp Z*{j*)Z*{j*) m'l') [fb] 


99 




0"NLO 


""NLO+ss 


LO NLO 


0"L0 


0"NLO 


O'std 


1.492(2) 


7.343(1) 10.953(2) 


1.49 


1.14 



Table 6: Cross sections for the gluon and quark scattering contributions to pp Z* {'y*)Z* (j*) — » £££'£' at the LHC 
= 14 TeV) where standard LHC cuts (75 GeV < Mu < 105 GeV, pre > 20 GeV, |r?i| < 2.5) are applied. The 
integration error is given in brackets. We also show the ratio of the NLO to LO cross sections and the ratio of the combined 
NLO+gij contribution to the NLO cross section. 

Mt = 170.9 GeV and M;, = 4.7 GeV. The pp cross sections are calculated at ^/s = 14 TeV employing 
the CTEQ6L1 and CTEQ6M [122] parton distribution functions at tree- and loop-level, corresponding 
to Ag^ = 165 MeV and A^^ = 226 MeV with one- and two-loop running for as(/x), respectively. The 
renormalization and factorization scales are set to Mz- 

We compare results for £££'£' production in gluon scattering with LO and NLO results for the quark 
scattering processes. Since we are interested in ^(7)^(7) production as a background, the gg ^ H ^ 
ZZ signal amplitude is not included. The LO and NLO quark scattering processes are computed with 
MCFM [60], which implements helicity amplitudes with full spin correlations [167] and includes finite- 
width effects and single-resonant corrections. Table[6]shows gluon and quark scattering cross sections for 
the LHC. We find a NLO A'-factor for qq ZZ of 1.5. Enhanced by the large gluon flux at the LHC, the 
gg process yields a 14% correction to the total ZZ cross section calculated from quark scattering at NLO 
QCD. This is substantially higher than the corresponding 6% increase for WW production [136], where 
no right-handed V f f coupling contributes. Relative to the LO qq — > ZZ prediction the gg contribution is 
about 20% in agreement with the finding in Ref. [162]. Without top and bottom quark contribution the gg 
cross section is 0.885(1) fb. The massive bottom and top loops increase the result based on intermediate 
light quarks by about 70%. This is much more than the corresponding 15% for gg — > WW [136], 
where all quark loops are suppressed by at least one top propagator. In the gg —>■ ZZ case on the other 
hand a pure b quark loop occurs and gives rise to a contribution that is similar to the massless first or 
second generation down quark loop. We estimate the remaining theoretical uncertainty introduced by 
the QCD scale by varying the renormalization and factorization scales independently between Mz/2 
and 2Mz- For the gg —>■ ZZ process we find a renormalization and factorization scale uncertainty of 
approximately 20%. The scale uncertainty of the qq —>■ ZZ process at NLO is approximately 4%. The 
scale uncertainties are thus similar for gg ZZ and gg — > WW. 

Selected differential distributions for pp —>■ ^(7)^(7) —>■ £££'£' at the LHC are shown in Fig. |9l 
where the standard set of cuts defined above is applied. Fig. |9^) shows the distribution in the invariant 
mass M41 of the four produced leptons. We compare the gluon-gluon induced contribution with the quark 
scattering processes in LO and NLO. We observe that the invariant mass distribution of the gluon-gluon 
induced process is similar in shape to the quark scattering contributions and suppressed by about one 
order of magnitude in normalization. Z boson pairs produced in quark-antiquark scattering at the LHC 
are in general strongly boosted along the beam axis. Gluon-induced processes on the other hand result 
in ZZ events at more central rapidities. This feature is born out by the distribution in the pseudorapidity 
of the negatively charged lepton shown in Fig. |9}5). In order to distinguish the shapes of the various 
contributions we have chosen a linear vertical scale and plot the gluon-gluon contribution multiplied 
by a factor 10. Compared to LO quark-antiquark scattering, the lepton distribution of the gluon-gluon 
process shows a more pronounced peak at central rapidities. We also observe an enhancement of the NLO 
corrections at central rapidities which is due to the substantial contribution of gluon-quark processes at 
NLO. To demonstrate the impact of the photon contribution, we show in Fig. [TO] distributions for the 
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Fig. 9: Distributions in tlie Ut! I' invariant mass M/a (a) and the pseudorapidity r\g- of the negatively charged lepton (b) 
for the gluon scattering process (sohd) and the quark scattering processes in LO (dotted) and NLO QCD (dashed) of pp 
Z* {'y*)Z* {-y*) £££'£' at the LHC. Input parameters as defined in the main text. Standard LHC cuts are applied (see main 
text and Table[6ll. The gg distribution of rj^- is displayed after multiplication with a factor 10. 

gluon-gluon induced process that include only ZZ, only 77 and all contributions. Here, a minimal set of 
cuts is applied, i.e. only Mu > 5 GeV in order to exclude the photon singularityJl With this minimal set 
of cuts the LHC cross section for gg Z{j)Z{-f) lll'l' increases to 7.874(5) fb. In Fig.[lOji) the Z 
and 7 contributions to the distribution in the invariant mass M^i are displayed. We observe that for Higgs 
masses below the Z-pair threshold, where one Z boson is produced off-shell, the photon contribution to 
the background is important. In Fig.fTOb) the contributions are shown for the distribution of the transverse 
momentum pj-^- of the negatively charged lepton. For this observable, the photon contribution becomes 
non-negligible for values below 70 GeV. 

6.3 Conclusions 

We have calculated the loop-induced gluon-fusion process gg — > Z*(7*)Z*(7*) lll'I' , which pro- 
vides an important background for Higgs boson searches in the H —>■ ZZ channel at the LHC. Our 
calculation demonstrates that the gluon-fusion contribution to the ZZ background yields a correction of 
about 15% to the qq prediction at NLO QCD and that the photon contribution is important for Higgs 
masses below the Z-pair threshold. We conclude that the gluon-gluon induced background process 
should be taken into account for an accurate determination of the discovery potential of Higgs boson 
searches in the pp ^ H ^ ZZ —>■ leptons channel. Our public program, named GG2ZZ, includes 
the ZZ, Z7 and 77 contributions with full spin and polarization correlations, off-shell and interference 
effects, as well as finite top and bottom quark mass effects. It is available on the Web [168] and can be 
used as Monte Carlo integrator or to generate unweighted parton-level events in Les Houches standard 
format [169, 170]. ATLAS and CMS are currently using our program to study the hadron-level impact 
of the gg — > ZZ background on H ^ ZZ searches at the LHC. 
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'^As in Ref. [162], a technical cut prz > 2 GeV is employed with standard cuts to exclude critical configurations. With 
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Fig. 10: Distributions in tlie Ut! I' invariant mass Mu (a) and the transverse momentum Pti- of ths negatively charged lepton 
(b) for the gluon scattering process gg Z* {'y*)Z* {j*) III' I' at the LHC with ZZ contributions only (dashed), 77 
contributions only (dotted) and all contributions (solid). Input parameters and minimal set of cuts as defined in the main text. 
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7. gluon- induced z*z* background simulation for higgs boson 
searchEI 

The contribution of the 55 ^ Z*(7*)Z*(7*) ^ ^1"//' process to the total ^ Z*(-i*)Z*{^*) Ul'i' 
production cross section has been evaluated after the application of the selection cuts adopted for the 
Higgs boson search through the H ZZ 2e2^ decay channel in the CMS experiment [171]. The 
minimal set of kinematical cuts needed to maximize the discovery significance has been used: upper 
thresholds for the transverse momenta {pTi) of the four produced leptons; upper threshold on the highest 
reconstructed M^^; lower threshold on the lowest reconstructed M^^; upper and lower thresholds on the 
M^n. The values of the selection cuts are mass dependent, optimized for different Higgs boson mass 
scenarios, from 120 GeV to 500 GeV. The selection procedure and the cut values are described in details 
in Ref. [171]. A sample of 9k gg Z*(7*)Z*(7*) Ul't' events has been generated at parton level 
using the generator program GG2 ZZ here presented. For the simulation of the shower evolution we have 
interfaced the generated parton-level events to the PYTHIA Monte Carlo generator [172]. In order to 
increase the event statistics in the kinematical region of interest the following set of pre-selection cuts 
has been used: pTt > 5 GeV, \ri(\ < 2.5, M^^ > 5 GeV. The cross section for the selected events is 2.8 fb. 
We compare the gluon induced contribution with 70k qq Z*(7*)Z*(7*) Ui'i' events generated 
with the MadGraph LO Monte Carlo generator [173]. The LO cross section of this process is 27.67 fb, 
where the same set of pre-selection cuts has been applied. In Fig [TT] we compare the distribution of 
the invariant mass (M41) of the four leptons produced in the gluon-gluon and in the quark scattering 
processes respectively. The peak at M41 ~ Mz in the qq distribution is due to the s-channel production 
process, that gives the main contribution to the cross section in the M41 mass region below and near Mz- 
Since the Higgs mass region below 114.4 GeV has been excluded by the LEP data [4], almost all the 
events produced by the s-channel process are removed by the selection cuts adopted in the Higgs boson 
search analyses. The effect of the mass dependent selection cuts on the M41 distribution is shown in 
Fig [12] The different curves correspond to the gg —>■ Z*(7*)Z*(7*) —>■ £££'£' events selected after the 
pre-selection cuts (solid curve) and for a few Higgs boson mass scenarios (dashed curves), when only 
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the cuts on the four leptons transverse momenta and on the di-lepton invariant masses (M^^, M^/^-/ ) have 
been appUed. The photon contribution is strongly suppressed for the Higgs boson search above 2Mz- 
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Fig. 1 1 : Distribution in tiie IWI' invariant mass, M^^i , for 
the gluon scattering process (solid) and ttie quark scatter- 
ing process (dashed) of pp ZZ ^ ill' I' at the LHC, 
after applying pre-selection cuts. 



Fig. 12: Selected distributions in the £U'£' invari- 
ant mass, M4e, for the gluon scattering process at 
the LHC, obteined applying the pre-selection cuts 
(solid) and the set of selection cuts optimized for 
the Higgs boson search in different mass scenarios 
from 120 GeV to 500 GeV. 



The contribution of the gluon scattering to the ZZ cross section is reported in Fig [131 in terms of 
the ratio of the selected gg Z*{Y)Z*{j*) Ui'i' events respect to the LO qq Z*{-f*)Z*{'y*) 
£££'£' selected events (solid square markers). The correction increases approximately linearly from 3% 
to 24% in the M41 region between 120 GeV and 200 GeV, and it is quite uniform, around ~ 24%, in 
the M41 region above 200 GeV. Superimposing to the graphic the ratio of the distributions shown in 
Fig [TT] (dashed curve), where only the pre-selection cuts were applied, we observe that the gluon induced 
contibution is enhanced by the selection cuts, especially in the mass region above 200 GeV. The empty 
round markers in Fig [13] show the gg contribution respect to the qq process calculated at the NLO. The 
mass dependent NLO k-f actor evaluated in Ref [174] has been used to rescale the quark scattering cross 
section at the NLO. The gluon-gluon contribution is reduced to 18% in the AI41 region above 200 GeV, 
a value compatible with our previous evaluation. 
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Fig. 13: Contribution of the gluon-induced background gg Z* {'y*)Z* {y*) £££'£' respect to the qg 

^*(7*)Z*(7*) —> £££'£' process, as a function of the four leptons invariant mass M41, after application of the pre-selection 
cuts only (dashed) and all selection cuts (other two curves). The reference cross section of the qg scattering is evaluated at the 
LO (solid markers) and at the NLO (empty markers). 



8. THE METHODS FOR THE CENTRAL RAPIDITY GAP SELECTION IN THE VECTOR 
BOSON FUSION SEARCHES IN CMS0 

8.1 Introduction 

In the VBF Higgs boson searches at LHC a selection of the events with the central rapidity gap between 
the two tagging jets is aimed to reduce the QCD Z+jets and other backgrounds like W+jets and ti while 
keep a high efficiency for the Higgs boson signal from the VBF production, VV H. The central jet 
veto was proposed and used in the first VBF Higgs boson analyses [113, 175] (see also references in it) 
and exploited in the recent, published ATLAS and CMS analyses [108, 176]. The central calorimeter 
jet veto technique is suffering from the pile-up and the electronic noise in the calorimeters which could 
create the fake jets. The method of the reduction of the fake calorimeter jets using the information 
from the event vertex and the tracks was proposed in [177] and successfully used in the CMS analyses 
[176,178]. 

We consider three methods to perform the hadron activity veto in the central rapidity region: the 
(traditional) central calorimeter jet veto (CJV), the track counting veto (TCV) and the veto on jets made 
from the tracks only (TJV). The idea of the track counting veto is inspired by the paper [179] where it 
was proposed to distinguish between the gluon and vector boson fusion processes for the Higgs boson 
production. The performance of methods is compared in terms of the signal efficiency and the QCD 
Z+jets background rejection. 

8.2 Studies at generator level 

The QCD Z+jets events were generated using the ALPGEN [180] generator with the MLM prescrip- 
tion for jet-parton matching [181, 182] in the PYTHIA6.4 shower generation [183]. The details on the 
ALPGEN generation and soft VBF preselections at the generator level can be found in [176]. 

The final VBF selections similar to the ones used in the full simulation analysis [176] were applied 
to the PYTHIA particle level jets. An event must have at least two leading Et jets reconstructed with a 
cone algorithm (cone size 0.5) that satisfy the following requirements: 

1. E^j, > 20GeV 

2. |r/J'| < 4.5 

3. Mjij2 > 1000 GeV 

4. lAr/Ji-^'^l > 4.2 

5. 1]^^ X 1]^'^ < 0. 

where jl and j2 are two leading £'t jets ordered in Et- 

The performance of the two methods, CJV and TCV was compared. The CJV requires to reject 
events with a third jet that satisfies 

• E^^ > 20 GeV 

, min _j_ ^ ^ max _ q 

where rf and tj^ are the minimum and maximum t] of the two leading jets (jl and j2). The 
TCV requires to reject events with a certain number of "tracks" (charged particles) within the tracker 
acceptance region, |?7| < 2.4 that satisfies 

• > p^^ GeV/c 

• T^i ™™ + 0.5 < 7?*™^'' < rj^ - 0.5, 

The effect of multiple parton interactions generated with Tune DWT [184] on the track counting veto 
was studied. 

Fig. O shows the number of charged particles within the tracker acceptance and between the two 
tagging jets {r]^ + 0.5 < r/*™^'^' < r]^ """"^ - 0.5) with pr > GeV/c (left plot), > 1 GeV/c (middle 
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plot) and > 2 GeV/c (right plot) for the signal (solid line) and the QCD Z+jets background (dashed 
line). The multiple parton interactions were switched off in PYTHIA. One can see a clear difference 
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Fig. 14: The number of charged particles within the tracker acceptance and between two tagging jets (t]^ + 0.5 < rf''^^^^ < 
^3 max _ ^-^j^ pT > GcV/c (left plot), > 1 GcV/c (middle plot) and > 2 GeV/c (right plot) for the signal (solid line) 
and the QCD Z+jets background (dashed line). The multiple parton interactions are switched off in PYTHIA. 



between the signal and the QCD Z+jets background distributions. This difference, however is spoiled 
when the multiple parton interactions are switched on. Fig.[T5]shows the same distributions as in Fig. [141 
but with the multiple parton interactions included in the generation. With no cut of the charged particle 
Pt applied, it is not possible to distinguish between the signal and the background. The cut on the 
"track" PT removes charged particles from the underlying event, thus giving the selection power for the 
TCV method. With the cut pf^=2 GeV/c the efficiency for the signal VV ^ H {Mh=120 GeV) is 
^ 0.8 and for the QCD Z+jets background is 0.54. For the same ^ SO % signal efficiency, the central 
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Fig. 15: The number of charged particles within the tracker acceptance and between two tagging jets (77-' + 0.5 < ff'''°-'^^ < 
max _ Q 5^ ^jjjj PT > GeV/c (left plot), > 1 GeV/c (middle plot) and > 2 GeV/c (right plot) for the signal (solid line) 
and the QCD Z+jets background (dashed line). The multiple parton interactions generated with Tune DWT are switched on in 
PYTHIA. 



jet veto efficiency for the background is smaller, 0.44, thus leading to the conclusion that at the particle 
level simulation the central jet veto provides the better performance than the track counting veto. The 
final conclusion, however should be resulting from the full detector simulation including the detector and 
reconstruction effects, like fake jet contribution from the pile up and the electronic noise, the track and 
jet reconstruction inefficiency. 



8.3 Studies with tlie full detector simulation 

The fully simulated datasets from the VBF Higgs boson analysis {H tt ^ £ + jet) [176] at an 
instantaneous luminosity L = 2 x 10'^^cm~^s^^ are used. The pile-up events (4.3 events per crossing) 
were included in the simulation. At the reconstruction level the same VBF selections 2-4 on tagging jets 
as described in the section 18.21 were used and the tagging jets were required to have >40 GeV. The 
CJV requires to reject events with a third jet that satisfies 

• > 10 GeV, where Et is a raw, non calibrated energy. 

• fake jet rejection parameter a-'^ = ptf"''^'' / E^ > 0.1 (see [176] for details) 

The TCV requires, on top of the selections mentioned in the previous section, the quality selections on 
the tracks: > 8 hits, AZ(track, vertex) <2 mm, Ai? (track, jet) >0.5. The lepton and tracks from r 
jet are not counted. 
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Fig. 16: Track multiplicity between the two forward tagging jets (ri^ + 0.5 < 7;*'''"='= < 77^ - 0.5) with pr > 1 GeV/c 
(left plot), > 2 GeV/c (middle plot) and > 3 GeV/c (right plot) for the signal (full circles) and the QCD Z+jets background 
(open squares). 



Fig. [16] shows the number of reconstructed tracks between the two forward tagging jets with pT > 
1 GeV/c (left plot), > 2 GeV/c (middle plot) and > 3 GeV/c (right plot). Both the Higgs boson signal 
(circles) and the QCD Z-i-jets background (squares) can be clearly separated when applying a cut on the 
track multiplicity and for different track pT thresholds. The left plot of Fig. [17] shows the performance 
of the algorithm, i.e the efficiency of selecting the signal versus the background. Starting from the 
bin on the left bottom comer, the points correspond to an increasing cut on the track multiplicity 
and pt up to the right top corner where 100% of events are selected. The black star indicates the 
performance of the central jet veto (CA'^) based on calorimeter jets. One can notice that this latter 
achieves a good performance: SO.Oit 3.3% efficiency for the signal and 39.7ib 5.% efficiency for the 
backgrouind. The TCV algorithm can reach this discrimination power rejecting events with more that 
one track of > 3 GeV/c. The right plot of Fig. [17] shows the ratio of the signal and the background 
selection efficiencies as a function of the signal efficiency. It shows that the better ratio can be achieved 
with the TCV at the price of losing a bit of signal. This would obviously depend on the overall tuning 
of the analysis. The third algorithm, the track-jet counting veto (TJV) is very similar to the TCV. Tracks 
are first clustered along the beam axis (Z) starting from the track with the highest px following the 
condition: AZ(cluster, track) <2 mm. Once a z-cluster is formed, the same procedure applies again 
with the remaining tracks. In a second step a traditional cone jet finding method is applied with AR = 0.5 
around seed tracks (with highest pt)- These jets are thus formed solely of tracks originating from the 
different z-clusters. They are finally associated with the signal vertex if their z-impact parameter is 
within 2 mm from the lepton z-impact parameter. This method allows to refine the description of the 
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Fig. 17: The track counting veto (TCV) performances (see text) for the different and track multiplicity thresholds. 



hadronization process usually producing several collimated particles with respect to the more exclusive 
approach of the TCV algorithm. The discrimination variable are then the multiplicity and the minimum 
Pt of the track-jets and its constituents lying in between the two forward tagging jets. The performance 
of this algorithm has been found to be very close to the TCV, reaching 80% for the signal efficiency and 
40% for the background when requiring no track-jet with AZ(track — jet, lepton) < 2mm, p^f* > 3 
GeV/c and with at least one track of px >0.9 GeV/c. 

8.4 The efficiency measurement of the central rapidity gap selection for Z ^ rr background. 

The efficiency of the central rapidity gap selection (CRGS) for the Z-f-jets, Z —>■ tt background in the 
VBF H TT search can be measured with the Z-f-jets, Z ^fi events passed the "signal like" VBF jet 
selections. We estimated the expected number of such events and the statistical accuracy of the CRGS 
for 100 pb^^ of integrated luminosity. Only QCD Z-i-jets events were used. The events from the EWK 
Z-i-2jets production still have to be added. The fully simulated events with no pile-up were required to 
pass the di-muon trigger. In the off-line analysis the events with two muons pT > 10 GeV/c, |r/| <2.4 
isolated in the tracker were selected within the di-muon mass window 70< M^^ <110 GeV/c^. The 
following VBF cuts relaxed for the early analysis with the first 100 pb~^ of the data are used. An event 
must have at least two leading Et jets that satisfy the following requirements: Ej. > 40 GeV, \r]^ \ <4.5, 
^ii X rj^^ < and: 

• soft VBF selections: Mjij2 > 400 GeV/c^, \Ari^'^^'^\ > 2.5 

• hard VBF selections: Mjij2 > SOOGeV/c^, \ Ar]^^^'^\ > 3.5 

The CJV used in this section requires to reject events with a third calorimeter jet that satisfies 

• E^ > 30 GeV, where Et is the calibrated jet energy 

• T]^ ™" + 0.5 < 7]^^ < r]^ ™<'^ - 0.5, 

Table |7] shows the expected number of events after selections for 100 pb~^ and the efficiency and 
the statistical accuracy of the CJV. Fig. [18] shows the distribution of rj^ and tj^ (left-upper plot), 
the 1]^^ (left-bottom plot) and the variable i]z=ri^^-0.5{rj^ + t]^ "*"^) (right plot) for 100 pb^^ of the 
"data" for one random experiment. All selections except the CJV were applied. 

8.5 Conclusions 

With the full detector simulation is was shown that both the central jet veto and the track counting 
algorithms achieve very similar performance. The robustness and the stability of the methods under a 
variation of the run and detector conditions will be tested with the real data using Z-f-jets, Z ^fi 
events. It is believed that the track counting algorithms relying on a single sub-detector, the tracking 
system, would perform with a higher reliability. 



Table 7: The expected number of events after selections for 100 pb ^ and the efficiency and the statistical accuracy of the CJV. 
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number of events with 100 pb ^ 


CJV efficiency 


"soft" VBF 


121 




"soft" VBF + CJY 


61 


0.50 ± 0.06 


"hard" VBF 


31 




"hard" VBF + CJV 


11 


0.35 ±0.11 




Fig. 18: The distribution of andr;-' (left-upper plot), the r;-'^ (left-bottom plot) and the variable ?7z=?7^^-0.5(r;-' + 

max^ (right plot). All selections except the CJV were applied. 
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9. PRODUCTION OF A HIGGS BOSON AND A PHOTON IN VECTOR BOSON FUSION AT 
THE LHC 

9.1 Introduction 

Higgs boson search is one of the main tasks of present and future collider experiments [2, 185]. At 
the CERN Large Hadron Collider (LHC), the Higgs boson is expected to be produced with high rate 
via gluon or vector-boson fusion (VBF) mechanisms and associate W{Z)H production. Apart from 
observing the Higgs signal, it would be crucial to make at the LHC also a measurement of the Hbb 
coupling [8]. To this aim, Higgs production via VBF, with the Higgs boson decaying into a bb pair, plays 
an important role [186]. In this contribution, we consider a further process that could help in determining 
the Hbb coupling, that is the Higgs boson production in association with a large transverse-momentum 
photon (with ^ 20 GeV) and two /orH'ar(/ jets [187] 



pp H'yjj bb-fjj + X , 
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(11) 



ruH (GeV) 


110 


120 


130 


140 




67.4 


64.0 


60.4 


56.1 


Bn{H bb) 


0.770 


0.678 


0.525 


0.341 



Table 8: Cross sections for the H j jj signal at LHC, for > 20 GeV, AR^j > 0.4, and a cut rrijj > 100 GeV on the invari- 
ant mass of the final quark pair. Also shown are the Higgs boson branching ratios to bb (computed through HDECAY [188]), 
that are not included in the cross sections shown. 

with H decaying to bb, where at the parton level the final QCD parton is identified with the coiresponding 
jet. In our study, we will not include diagrams where the photon is emitted from one of the two b-quarks, 
since the requirement of a large photon would shift in that case the bb invariant mass outside the 
experimental bb mass resolution window around the Higgs mass. 

There are a number of advantages in considering this QED higher-order variant of the VBF Higgs 
production process pp ^ H{-^ bb) jj. The fact that the production rate is penalized by the electromag- 
netic coupling is compensated by a few peculiarities of the channel in Eq. ([iTI ). First of all, the presence 
of an additional high px photon can improve the triggering efficiencies for multi-jet final states, such as 
those needed to select pp — > H{-^ bb) jj events. Second, there is a large gluonic component entering 
the partonic processes giving rise to the QCD backgrounds to the bb jj final state. As a consequence, 
the QCD backgrounds are in general much less active in radiating a large px photon with respect to the 
VBF signal. In addition there are dynamical effects that dramatically suppress the radiation of a central 
photon in the irreducible QCD background to bb 7 jj with respect to the VBF channel. When the photon 
is forced to be emitted in the central rapidity region, a destructive quantum interference arises between 
the photon emission off the initial quark exchanging a gluon (or any other neutral vector boson) in the 
t channel, and the photon emission off the corresponding final quark. For the signal case of the H jj 
production, the above mechanism of destructive interferences affects only the diagrams involving the ZZ 
fusion. On the other hand, in the diagrams involving WW fusion (that are responsible for the dominant 
part of the basic VBF H jj cross section) the charged currents in the qq'W vertices change the electric 
charges of the in-out partons, and consequently the interference is now additive rather than destructive. 
Therefore, the cross section for H 7 jj is expected to follow the usual pattern of QED corrections as far 
as its WW fusion component is concerned. The relative contribution of the ZZ fusion will be instead 
remarkably smaller than in the case of the basic VBF H jj process. 

To summarize, a measurement of the bb 7 jj rate could lead to a combined determination of the 
Higgs boson couplings to b quarks and W vector bosons, with less contamination from the HZZ cou- 
pling uncertainties. 

In Section 19. 2[ we go through the main kinematical and dynamical characteristics of the process 
in Eq. (fTTT) . We also discuss the features of the main QCD irreducible background. In Section 1931 the 
signal rates are computed at parton level for a set of kinematical cuts that optimizes the signal/background 
ratio, restricting the analysis to the case of the irreducible background. In Section l9!4l the main reducible 
background channels are included in the analysis. Finally, in Section [931 we draw our conclusions. 

9.2 Signal and Irreducible Background 

Cross sections for the H 7 jj production at = 14 TeV are shown in Table[8] In order to present results 
as inclusive as possible only a minimal set of kinematical cuts is applied (AR^j > 0.4, p^ > 20 GeV, 
and nijj > 100 GeV). The Higgs boson branching ratios to bb, which are not included in the cross 
section results, (computed through HDECAY [188]), are also shown. The full tree-level matrix elements 
for the electroweak process pp —>■ H j jj have been computed independently with ALPGEN [180], and 
MadEvent [173]. Details on the values of input parameters, such PDF's and scales are given in Section 



Before discussing the process in Eq. (fTTl) . it is useful to recall here the main kinematical properties 
of a typical VBF event, that is pp —>■ H jj, and the corresponding backgrounds. For the Higgs boson 
decaying to a bb pair, the main background to the basic VBF process comes from the QCD production 
of the final state bbjj, whenever the bbjj kinematical characteristics approach the typical VBF configu- 
ration. By imposing a large invariant mass cut for the two-forward-jet system [i.e., rrijj > 0(1) TeV], 
a minimal of a few tens GeV's, and requiring the bb invariant mass to be around niH within the m^,^ 
experimental resolution, one can obtain a signal significance (S/^/B) of the order of S/^/B ~ 3 — 5, 
assuming an integrated luminosity of 600 fb~^ [186]. 

Let us now consider the VBF Higgs production when a further central photon is emitted, namely 
pp — > H-yjj. According to the usual pattern of QED corrections, one might expect the request of a 
further hard photon to keep the relative weight of signal and background quite stable. Were this the case, 
the rates for pp — > H j jj and its background would be related to a 0{a) rescaling of the rates for the 
H jj signal and its background, respectively, where a is the fine electromagnetic structure constant. On 
this basis, one would conclude that there is no advantage in considering the H 7 jj variant of the H jj 
process, apart from the fact that the presence of a hard photon in the final state can improve the triggering 
efficiency of the detectors. However, as we explained in the introduction, this pattern does not hold in 
general when restricted regions of phase space are considered. 

In the next section we will study this effect on a quantitative level, showing that the requirement 
of a further central photon gives rise to a dramatic increase (by more than one order of magnitude) in the 
S/B ratio, while the signal cross section roughly follows the naive QED rescaling. 

9.3 Cross Sections for the Signal versus the Irreducible Background 

The numerical results presented in this section have been independently obtained by the Monte Carlo 
event generators ALPGEN [180], and MadEvent [173]. The signal is calculated in the narrow width 
approximation, i.e. we computed the exact lowest-order matrix element for the process pp jj, 
and then let the Higgs boson decay into a bb pair according to its branching ratio and isotropic phase 
space. After the decay, cuts on the 6— quark jets are implemented. For the irreducible pp bfry jj 
background, we computed all the matrix elements at 0(af a), neglecting 0{ala^), O(a^a^), ©(a^a^) 
and 0{a^) contributions and their interference with the ©(a^a) contribution. We checked that this has 
no numerical impact on the results. The present study is limited at the parton level. A more complete 
simulation, that takes into account showering, hadronization and detector simulation, even if crucial for 
the assessment of the potential of this channel, is beyond the scope of the present contribution. As PDF's, 
we use the parametric form of CTEQ5L [189], and the factorization/renormalization scales are fixed at 
= /^R = S ^iid /^F = Mr = + for the backgrounds and signal, respectively {Et is the 

transverse energy of any QCD parton). The three Higgs-mass cases 120, 130 and 140 GeV are analysed. 

We start by the definition of two basic event selections (sets 1 and 2) that differ only by the 
threshold on the photon transverse momentum p^: 

j4 > 30 GeV, > 30 GeV, AR,k > 0.7, 
^7! < 2.5, \vb\ < 2.5, < 5, 

rrijj > 400 GeV, mnil - 10%) < m^^ < "^^(1 + 10%), 

1) pI > 20 GeV, 

2) pI > 30 GeV, (12) 

where ik is any pair of partons in the final state, including the photon, and ARik = \/ A?rjik + A?(f)ik, 
with rj the pseudorapidity and (f) the azimuthal angle. The cross sections for the above basic event 
selections are reported in Table |9l Before comparing the signal and the background for the H 7 jj 
process, we tried to optimize our event selection in Eq. (fT2l) . Indeed, the signal detectability can be 
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Table 9: Cross sections for the signal and the irreducible background for the basic event selections in Eq. il2l . Higgs production 
cross sections include the Higgs branching ratios to bb. The signal and irreducible background production rates for the VBF 
process without photon are also shown. 
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Table 10: Cross sections for the signal and the irreducible background for the optimized event selections of Eq. l ll3t . added 
to the basic selection in Eq. l ll2b . Higgs production cross sections include the Higgs branching ratios to bb. The signal and 
irreducible background production rates for the basic VBF process are also shown. 

further improved by imposing optimized cuts, that can be deduced by looking at the following kinematical 
distributions: 

da da da da da 

drrijj' dpi^'' dp^' dm^yn' 

where jl and 61 denote the leading px hght jet and b— jet, respectively, and m^n is the invariant mass 
of the 'jbb system. By studying the variation of the significance S/ \fB as a function of the cuts on the 
distributions (for more details see [187]), we found an optimized event selection where, in addition to the 
basic cuts, we impose the following cuts 

> 800 GeV, p'^ > 60 GeV, > 60 GeV, 
\Ar]jj\ > 4, m^H > 160 GeV, AR^hf-yj > 1-2. (13) 

With the above additional requirements, we find the cross sections reported in Table [TOl One see that the 
requirement of the extra central photon with p^ > 20 GeV in the final state involves a reduction factor of 
order 100 for the signal rate with respect to the final state without photon, according to the expectations 
of the 0{a) QED naive scaling. On the other hand, the radiative background is suppressed by a factor 
of about 3000 with respect to the case of no photon radiation. Finally, a summary of the statistical 
significances, including only the irreducible background, with an integrated luminosity of 100 fb^^ is 
given in Table [TTl 

9.4 Reducible Backgrounds 

A complete analysis of the reducible backgrounds to the H 7 jj signal is beyond the scope of our study. 
However, in order to have a sensible estimate of the achievable S/B ratio and statistical significance at 
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Table 11: Statistical significances with the event selection of Eq. l ll2b and i\3l . with an integrated luminosity of 100 fb^^. 
The value et = 60% for the 6— tagging efficiency and a Higgs boson event reduction by e^;, ~ 70%, due to the finite (±10%) 
bb mass resolution, have been assumed. Jet- tagging efficiency and photon-identification efficiency are set to 100%. Only the 
irreducible background is included in this analysis. 



parton level, we computed with ALPGEN the cross sections, assuming mn = 120 GeV and with the 
optimized event selection of Eq. (fT2l ) and (fT3] ). for three main potentially dangerous processes 

• pp —>■ ^ + 4 jets, where two among the light jets are fake tagged as 6— jets; 

• — > 65 + 3 jets, where one of the light jets is misidentified as a photon; 

• — > 5 jets, where one of the light jets is misidentified as a photon, and two light jets are fake 
tagged as 6— jets. 

By including also the reducible backgrounds, the statistical significance decreases by about 14(12)% 
for p^^""^ = 20(30) GeV) with respect to Table [TT] where only the irreducible background has been 
considered. The most dangerous contribution to reducible backgrounds comes from pp bb + 3j. 

9.5 Conclusions 

In this contribution, we studied the detectability of the Higgs boson production signal, when the Higgs 
boson is accompanied by a high— px central photon and two forward jets at the LHC. The Higgs boson 
decay into a bb pair is considered. We analyzed the signal, the irreducible QCD background, and main 
reducible backgrounds at the parton level. The presence of a photon in the final state can improve the 
triggering efficiencies with respect to the basic VBF Higgs production without a photon. Moreover, 
we find that the requirement of a central photon in addition to the typical VBF final-state topology 
significantly suppresses the irreducible QCD background. In particular, due to dynamical effects, the 
latter has rates that are lower than the expectations of the 0{a) QED naive scaling by more than an order 
of magnitude. As a consequence, after optimizing kinematical cuts, we obtain a statistical significance 
S/^/B for the > bb)"fjj channel that goes from around 3, if rriH ~ 120 GeV, down to about 1.5, 
if niH — 140 GeV, for an integrated luminosity of 100 fb~^. These significances are not far from the 
corresponding values for the basic H{—>^ bb) jj process without a photon. The latter estimates are based 
on the irreducible QCD background. The impact of including a few main reducible backgrounds has 
found to be moderate. The same dynamical effects that are responsible for the irreducible background 
suppression also remarkably curb the relative contribution of the ZZ H boson fusion diagrams with 
respect to the WW — > H ones in the process pp —>■ H{-^ bb)j jj. As a consequence, we think that the 
study of the -ff(— > bb)'y jj signal at the LHC could have a role in the determination of both the Hbb and 
HWW coupUngs. 

10. THE Z PLUS MULTI-JET BACKGROUND FROM THE DOUBLE PARTON INTERAC- 
TIONS IN THE VECTOR BOSON FUSION H t+t- SEARCH^ 

The Z-f-jets production is the dominant background obtained in the VBF H —>■ tt searches at the LHC 
[108, 113, 175, 176]. We estimated an additional Z-i-jets background originated from double parton inter- 
actions (DPI) in a proton-proton collision when the Z boson is produced in one parton-parton interaction 
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and the QCD di-jets are produced in the second parton-parton interaction. In that case the two choices of 
the tagging jets are possible: (a) one tagging jet is taken from the QCD di-jet production and the second 
one is taken from the Drell-Yan production and (b) two tagging jets are both selected from the QCD 
di-jet production. 

The contribution from the double-parton interaction was estimated with PYTHIA6.4 [183] at the 
particle level At the first step the Drell-Yan and the QCD di-jet events were generated separately 
in PYTHIA. The Drell-Yan production was generated with the full underlying event (UE) using Tune 
DWT [184], while in the QCD di-jet production the UE was switched off (MSTP(81)=0). The Drell- 
Yan events were generated with the di-lepton mass mu > 70 GeV/c^ and the QCD di-jet events were 
generated with pr >20 GeV/c. The NLO cross section 2x10^ fb for the Drell-Yan production and the 
PYTHIA cross section 8.2 x 10^^ fb for QCD di-jet production were used in the estimates presented. At 
the second step two events (Drell-Yan and QCD di-jets) were mixed together and analyzed as one event. 
Jets were found at the particle level by the simple cone algorithm (cone size 0.5) implemented in the 
PYTHIA PYCELL routine. 

The cross section for the double parton interactions was evaluated with the factorization formula 

m aAXCTB 

crA,B = — , (14) 

where m=l, for indistinguishable parton processes and m=2 for distinguishable parton processes (in 
our case we use m=2). In the experimental study of double parton collisions CDF quotes (j^ff=l4.5 
mb [190]. For LHC energy we use currently the value cre//=20 mb IH It gives the o-^ £)=8.2x 10^ 
fb (j4=Drell-Yan, i?=QCD di-jets). More pessimistic value of 12 mb o will double our estimates of 
the Z-i-jets background from the double parton interactions. The longitudinal correlations in the double- 
parton structure functions neglected in the above formula can have a sizable effect at the LHC [191, 192]. 

We compare the Z-i-jets background from the double parton collisions with the "normal" QCD 
Z-f-jets background from one parton-parton collision. It was generated using the ALPGEN [180] genera- 
tor with the MLM prescription for jet-parton matching [181, 182] at the PYTHIA6.4 shower simulation. 
We generated ^i'-i-2jets exclusive and ^£-1-3 jets inclusive samples with the ALPGEN settings: mu > 
70 GeV/c^ p{, >20 GeV/c, |r/j | < 5, /^Rjj > 0.5. The user "soft" VBF pre-selections in ALPGEN 
generation were: Ar]jij2 >4, r]ji x r]j2 <0, Mjij2 >600 GeV/c^, where jl and j2 are two leading pT 
partons ordered in pT- The parameters for MLM jet-parton matching were: £'|?"*=20 GeV, i?'^'"'*=0.5 
and r]"^ ^'^^ =5.0. 

The VBF selections similar to that were used in a full simulation analysis [176] (except cut on 
Et of the tagging jets, which is lower here) were applied to the PYTHIA particle level jets. An event 
must have at least two leading Et jets reconstructed with a cone algorithm (cone size 0.5) that satisfy 
the following requirements: 

• E^j.>20 GeV 

• 1]^ < 5.0 

• Mjij2 > 1000 GeV 

• lArz-^'i-''^! > 4.2 

• rj^^ X T]-?^ < 0. 

where j 1 and j2 are two leading Et jets ordered in Et- The double parton scattering events where the 
two leading jets were both originated from the Drell-Yan production (the fraction of such events is ~ 
20%) were excluded from the consideration to avoid the double counting with "normal" QCD Z-f-jets 
background. 
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Table [12] shows the initial cross sections (in fb) for the Z+jets background from one and two 
parton-parton interactions and cross sections after the VBF cuts. After selections the contribution from 

Table 12: The initial cross sections (in fb) for the Z+jets background from one and two parton-parton interactions and cross 
sections after the VBF cuts. 



interaction 


one parton-parton 


two parton-parton 


process 


exclusive ££+2] 


inclusive £1+3] 


Drell-Yan 


QCD di-jets 


no cuts 


1.0x10^ 


2.0x10^ 


8.2x10^ 


> 2 jets, E^j, > 20 GeV 






4.0x10^ 


^Vji,j2 >4.2, Tjji X r]j2 <0 


2.4x102 


5.3x102 


5.0x10^ 


Mjij2 >1000GeV/c2 






3.2x10^ 



the double parton interactions is ~ 40% (320 fb) of the "normal" Z-f-jets background from the one 
parton-parton interactions (770 fb). Fig. [19] (left plot) shows an angle in the transverse plane between 
two tagging jets (Aipjij2) for the Z-i-jets backgrounds and the signal VV —>■ H. All distributions are 
normalized to unit. We have obtained that the fraction of the DPI events when the one tagging jet is 

DY and one from QCD di-jets 
from QCD di-jets 



<0 
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Fig. 19: Left plot: the angle in the transverse plane between two tagging jet for Z+jets backgrounds and the signal VV — > H. 
Right plot: the angle in the transverse plane between two tagging jet for Z+jets background from the DPI for the cases (a), (b) 
and total (see the text). All distributions are normalized to unit. 

selected from the Drell-Yan and another from the QCD di-jet production (case (a)) is ~ 70%; in the 
rest 30% of the DPI events the both two tagging jets are selected from the QCD di-jet production (case 
(b)). Fig. [19] (right plot) shows the Aipjij2 distributions for the cases (a) and (b) separately as well as 
their sum (the same curve as in Fig. [19] (left plot)). One can see that in the case (a), as expected there 
is no any correlations between two tagging jets, while in the case (b) they are forming the back-to-back 
configuration. 

Fig.[20](left plot) shows the transverse energy of the tagging jets from the Z-f-jets backgrounds and 
the signal VV — > H. One can see that the Z-f-jets background from the DPI can be largely suppressed 
with the cut on the tagging jet energy Ej, >40 GeV. This cut was used in the full simulation analysis 
[176]. After applying this selection the cross section of the Z-i-jets from the DPI is ~ 100 fb and the 
cross section of the "normal" Z-f-jets background is ~700 fb, thus the relative contribution from the 
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DPI is reduced to ~ 15 %. The further reduction of the relative Z+jets contribution from the DPI is 
expected when the cuts on the momentum of the lepton (from r — > Ivv decay) and the hadronic r jet 
(from T hadrons v decay) will be applied. It is due to the momentum of the Z boson from the DPI 
is softer that the one from the "normal" Z+jets production. It is shown in Fig. |20] (right plot) where the 
distributions of from DPI and the "normal" Z+jets production are normalized on the expected cross 
sections after the VBF cuts. 
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Fig. 20: Left plot: the transverse energy of the tagging jets from the Z+jets backgrounds and the signal VV H\ distributions 
are normalized to unit. Right plot: the distributions of the Z boson transverse momentum from the DPI and the "normal" 
Z+jets production are normalized on the expected cross sections after the VBF cuts. 



It is important to control and measure the Z+jets background from the double parton interactions. 
The possibility of the usage of the Z+2jets, Z— > ///x events with the VBF jet selections and looking at 
the unbalance in ~pT between the Z boson and the jets is under investigation. 

10.1 Conclusion 

The Z+jets background from the double parton interaction was estimated at the particle level to be 
less than 15% of the "normal" QCD Z+jets background in the VBF H —>■ tt searches at LHC after 
the experimental like event selections and assuming the a^j /=20 mb in the factorization formula. The 
fraction of the DPI events when the one tagging jet is selected from the Drell-Yan and another from the 
QCD di-jet production is ~ 70% while in the rest 30% of the DPI events the both two tagging jets are 
selected from the QCD di-jet production. 
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Part III 

MSSM HIGGS BOSONS 



11. SUSY-QCD CORRECTIONS TO SQUARK LOOPS IN GLUON FUSION TO HIGGS 
BOSONS0 

11.1 Introduction 

In the MSSM 2 Higgs doublets are introduced to generate masses of up and down type quarks. After 
electroweak symmetry breaking this leads to five physical Higgs particles, two light CP-even, /i, H, 
one CP-odd, A, and two charged H^. At tree level the MSSM Higgs sector can be described by 2 
independent parameters, usually chosen as the pseudoscalar mass Ma, and the ratio of the 2 vacuum 
expectation values tg/3 = V2/V1. The MSSM Higgs couplings to quarks are modified such that the 
couplings to down(up)-type quarks rise(decrease) with tg/3. The main neutral Higgs production at the 
Tevatron and LHC proceeds via gg fusion. The next-to-leading order (NLO) QCD corrections to this 
process have been known for a long time [11, 193] including the full quark mass dependence. They 
turn out to be important, increasing the cross section by up to 100%. Next-to-next-to leading order 
(NNLO) corrections, calculated in the large quark mass limit only [15-17, 194, 195], add another 20- 
30% and next-to-next-to-next-to leading order (NNNLO) corrections have been estimated [19, 196, 197], 
indicating improved perturbative convergence. NLO corrections to squark loops have been known so 
far only in the heavy squark mass limit [198], and the full SUSY-QCD corrections have been obtained 
for heavy SUSY masses [25-27, 199]. As a first step to a full SUSY-QCD result we present the QCD 
corrections to squark loops including the full squark and Higgs mass dependences [28]. 



11.2 NLO QCD Corrections 

For our calculation of the pure QCD corrections to squark loops we need a modified MSSM Lagrangian 
which separates the gluon and gluino contributions in a renormalizeable way. In this work we do not 
take into account the self-interaction among squarks, and the required Lagrangian is then given by 

C = -}-G''^^'^G1,--F^^''F^, + \[{d,,nf-Ml{H^] (15) 



Q ^ Q 



2 

Tn ~ 

\DA''-^%\Q?-9^^\Q?n 



with the covariant derivative = + igsG'^^T'^ + ieA^Q. Here G^j^ denotes the gluon field strength 
tensor and GJJ the gluon field accompanied by the color SU (3) generators T*^ (a = 1, . . . , 8), while F^^ 
is the photon field strength tensor and Af^ the photon field associated by the electric charge operator Q. 
The Higgs field H represents generically either the light scalar h or the heavy scalar H Higgs boson of 
the MSSJvT^ Since we do not take into account gluino exchange contributions, the coefficients gQ and 
g^ are not renormalized, thus leading to a renormalizeable model with strongly interacting scalars Q. 
Gluino corrections are expected to be small [25-27, 199]. 

For our numerical results we choose the gluophobic Higgs scenario [200] which maximizes the 
destructive interference effects between top and stop loops in the light Higgs coupling to gluons. It is 
defined by the MSSM parameters [mt = 174.3 GeV] MsusY = 350 GeV, i_i = M2 = 300 GeV, 
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"since there are no squark loop contributions to the pseudoscalar Higgs boson couplings to photons and gluons at leading 
order (LO), in this paper we will only deal with the scalar Higgs bosons h, H. 



Xt = At — /u/tg/3 = —770 GeV, Ah = At and nig = 500 GeV. The squark masses are given by 



tgP = 3: mi^ = 156 GeV tg/3 = 30 : m^^ = 155 GeV 

mi^ = 517 GeV m^-^ = 516 GeV 

= 346 GeV m^^ = 314 GeV ^ ' 

m-' = 358 GeV m-' = 388 GeV . 

02 02 

The results of this work look similar in other scenarios, whenever the squark masses are of the order of 
the top mass, or the Higgs mass reaches values beyond the corresponding squark-antisquark threshold. 



11.21 Scalar Higgs couplings to photons 

The leading order photonic Higgs couplings are mediated by top, bottom and W boson loops, with 
significant contributions from squark loops for stop and sbottom masses below ~ 400 GeV [2, 11, 115, 
193, 201-203]. The reverse processes -yj —>■ h, H play an important role for the MSSM Higgs boson 
searches at a photon collider [204-210]. The two-loop diagrams of the QCD corrections to squark loops 
lead to 5-dimensional Feynman parameter integrals. We have reduced these integrals in one calculation 
to 1 -dimensional integrals which have been integrated numerically. A second calculation has solved the 
integrals purely numerically. The two calculations agree within integration errors. In order to improve 
the perturbative behaviour of the squark loops they have been expressed in terms of the running squark 
masses mQ(M^/2), which are related to the pole masses Mq via mQ{fi) = MQ{as{fi) / ctsiMQ))^/ ^° 
where /3o = 33 — 2Nf with Np = 5 light flavors. Their scale is identified with /i = M>^/2 within the 
photonic decay mode thus insuring a proper definition of the QQ thresholds iV% = 2Mg. The LO scale 
dependence of the squark masses due to light particle contributions has been taken into account. 
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Fig. 21: Relative QCD corrections to the scalar MSSM Higgs decay widths to two photons for tg/9 = 3 and 30. The full curves 
include all loop contributions while in the dashed lines the SUSY contributions are omitted. The kinks and spikes correspond 
to the WW, titi , tt, bibi , TiTi , r2r2 and 6262 thresholds in consecutive order with rising Higgs mass. 

Figj2T] shows the relative QCD corrections to the photonic Higgs decay widths for the two cases, 
in which SUSY particles have been taken into account or not. The spikes which appear at the QQ 
thresholds are due to singularities originating from Coulomb singularities at the threshold since QQ 
pairs can form 0^^ states. This behaviour can be derived quantitatively from the Sommerfeld rescattering 
corrections, and we checked explicitly that this agrees with our numerical results. As can be inferred from 
Figl2T]the QCD corrections reach a size of 10-20% for moderate and large Higgs masses apart from the 
threshold regions, where the perturbative results are unreliable due to the Coulomb singularities. At a 77 
collider the photon fusion cross section can be measured with few per cent accuracy, and therefore these 
corrections have to be taken into account properly. The size of the QCD corrections with and without 
SUSY particle loops is of the same order of magnitude, but they can be of opposite sign. 

Figj22j in which the ratio of the fully massive photonic decay width at NLO T{h/H 77) and 
the NLO width with the relative QCD corrections in the heavy squark mass limit Too is plotted, quantifies 
the size of the squark mass effects beyond the heavy squark mass limit in the relative QCD corrections. 
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Fig. 22: Ratio of the QCD corrected partial decay widths of the scalar MSSM Higgs bosons to two photons including the full 
squark mass dependence and those obtained by taking the relative QCD corrections to the squark loops in the heavy mass limit 
as functions of the corresponding Higgs masses for tg/3 = 3 and 30. 



(The full squark mass dependence in the LO width has been kept in both expressions.) With a size 
of up to ~ 30% the squark mass effects are larger than the expected experimental uncertainty in the 
measurement of the Higgs production in 77 fusion and hence have to be taken into account in realistic 
analyses. 

11.22 Gluon Fusion 

The gluon fusion processes gg h, H axe mediated by quark and squark triangle loops with the latter 
contributing significantly for squark masses below ~ 400 GeV. The NLO QCD corrections consist of 
virtual two-loop corrections and the real corrections from the radiation processes, gg gh/H, gq 
qh/H and qq — > gh/H. The strong coupling constant has been renormalized in the MS scheme, 
with the top quark and squark contributions decoupled from the scale dependence, and the quark and 
squark masses in the on-shell scheme. The parton densities are defined in the MS scheme with five 
active flavors, i.e. the top quark and the squarks are not included in the factorization scale dependence. 
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Fig. 23: Production cross sections of the scalar MSSM Higgs bosons via gluon fusion as functions of the corresponding Higgs 
masses for tg/3 = 3 and 30. The full curves include the QCD corrections, while the dashed lines correspond to the LO 
predictions. The kinks and spikes correspond to the titi, bibi and 6262 thresholds in consecutive order with rising Higgs mass. 

Fig. |23] shows the LO and NLO cross sections. The QCD corrections increase the gluon fusion 
cross sections by 10-100% and are significantly larger in regions of large destructive interferences be- 
tween quark and squark loops. The corrections are of very similar size for the quark and squark loops 
individually. In spite of the large corrections the scale dependence is reduced from about 50% at LO to 
~ 20% at NLO thus indicating a significant stabilization of the theoretical predictions. Based on this and 
the approximate NNLO and NNNLO results in the limit of heavy squarks and top quarks the residual 
theoretical uncertainties of our NLO results can be estimated to less than about 20%. The spikes at the 
QQ thresholds are Coulomb singularities due to the formation of O"*""*" states. 

Fig. [24l which shows the ratios of the NLO cross sections including the full mass dependence 
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Fig. 24: Ratio of the QCD corrected production cross sections of the scalar MSSM Higgs bosons via gluon fusion 
full squark mass dependence and those obtained by taking the relative QCD corrections to the squark loops in th 
limit as functions of the corresponding Higgs masses for tg/3 = 3 and 30. 
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and of the NLO cross sections in the heavy squark Umits, exempUfies the squark mass effects on the 
K factors. In addition to the LO squark mass dependence of the cross section, the K factors develop a 
squark mass dependence of up to about 20% and hence support the relevance of our results compared 
to the previous results of Ref. [198]. The squark mass effects on the K factors are larger than the 
corresponding quark mass effects [211]. And they are larger than the residual theoretical uncertainties 
so that they cannot be neglected in realistic analyses. Since the gluino contributions are expected to 
be much smaller, the squark mass dependence will be the dominant part of the differences between the 
heavy mass limits and a full MSSM calculation at NLO. 

11.3 Conclusions 

We have discussed the NLO QCD corrections to the squark loop contributions to neutral MSSM Higgs 
boson production in gg fusion at the LHC and their decay modes into photons, including the full mass 
dependences. The corrections are sizeable and stabilize the theorectical predictions compared to the LO 
results. Squark mass effects on the relative QCD corrections are significant and larger than the mass 
effects from quark loops. They are always relevant for Higgs masses beyond the corresponding virtual 
squark-antisquark threshold or for squark masses of the order of the top mass. Since they are larger than 
the experimental uncertainties and the approximative results beyond NLO indicate sufficient perturbative 
convergence, the results of this work have to be taken into account for realistic analyses. 

12. HIGGS BOSON PRODUCTION IN ASSOCIATION WITH b QUARKS: SUSY QCD CON- 
TRIBUTIONS 

12.1 Introduction 

In the Standard Model, Higg production in association with b quarks is never important. However, 
in the minimal supersymmetric model (MSSM), the couplings of the Higgs bosons to b quarks can 
be significantly enhanced for large values of tan /3 and for a large range of parameter space, Higgs 
production in association with b quarks is the most likely discovery channel [83, 86, 89, 212, 213]. 

The production of a Higgs boson in association with a b quark has been extensively studied. In a 
4- flavor number scheme, the lowest order processes for producing a Higgs boson and a b quark are gg — > 
bb(j) and qq bbcf) [83,85,214] (The neutral Higgs bosons of the MSSM are generically ct> = h^, H^,A^). 
In a 5- flavor number scheme, the b quark appears as a parton and potentially large logarithms of the form 
In(^) are absorbed into b quark parton distribution functions. Although the 4- and 5- flavor number 
schemes represent different orderings of perturbation theory, the two schemes have been shown to yield 
equivalent numerical results. In the 5- flavor number scheme, the lowest order process for producing a 
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Higgs boson in association with b quarks is 66 ^ when no 6 quarks are tagged in the final state and 
bg b(j) when a single 6 quark is tagged [83, 85, 86, 89,214]. 

In this note, we consider the production process, bg — > bcp, for which the NLO QCD corrections 
are well understood, [86, 89,215]. Here we present the ©(a^) SUSY QCD (SQCD) corrections from 
gluino-squark loops to the 6- Higgs production cross section [93, 216]. We compare the results from 
an effective Lagrangian approach with those obtained from an exact one-loop calculation. Finally, we 
consider whether the process bg b(p +jet provides a useful signature and compare this channel with 
the irreducible background from bg —>■ bZ +jet. 



12.2 Effective Lagrangian 

The MSSM contains two Higgs doublets, Hu and H^. At tree level, the 6 quark couples to only one of 
the Higgs doublets {Hd) and there is no ipibuHu coupling (where ipL = {tb, 6l)). The coupling of the 
6 quark to the "wrong" Higgs doublet at one-loop leads to the effective interaction [217,218], 

L,ff = -AfcV^i {Hd + ^^«) + h.c. . (17) 

This effective interaction shifts the 6 quark mass from its tree level value, 

m, = ^(l + Am,), (18) 

and also implies that the Yukawa couplings of the Higgs bosons to the 6 quark are shifted from the 
tree level predictions. The shift of the Yukawa couplings [92] can be accounted for using an effective 
Lagrangian [217,219], 

L^ff = -^f-^V-^Vl- ^."^^ \bbh' 



vsMx^ + ^fnbJx COS P J \ tan /? tan a 
nih f 1 \ / cos a\ / Amh tan a 



1 + bbH' 



rO 



vsM V 1 + Am^ J \ cos P J \ tan /3 



vsM V 1 + Am^ J \ J \ tan^ /3 ^ 

^ 9,lJ>bh^ + 9,lHbbH'' + %lA^^l^bA\ (19) 

where vsm = 246 GeV, tan P = Vu/vd, and a is the mixing angle which diagonalizes the neutral Higgs 
boson mass matrix. The Lagrangian of Eq. [19] has been shown to sum all terms of 0(a" tan" P) for 
large tan/3 [217]. 

The expression for Am^ is found in the limit << M^,Mz « m-^^,m-^^,mg, (where 
"^fei ' "^62 ' ^9 sbottom and gluino masses) . The contribution to Am;, from sbottom/gluino loops 

is [218,220] 

where the function /(a, 6, c) is. 



// is the bilinear Higgs mixing parameter and OsifJ-R) should be evaluated at a typical squark or gluino 
mass. Note that Eq. |20]is valid for arbitrary values of tan p. 



Eq. |20]is a non-decoupling effect: If the masses of the squarks and gluino, along with the mixing 
parameter /x, become large for fixed Ma, Am^ does not vanish, 



In the large Ma limit, 



Ami, ^ —sign{fi) ^^anP + cot a] . (22) 



2M^ / j\^4 ■ 

tan/3 + cota^ — tan/3cos 2/3 + O ) , (23) 



and the decoupling limit of the MSSM is recovered [221]. 

The effective Lagrangian can be used to approximate the squark and gluino contributions to the 
rate for bg — > b(p [93]. We define an Improved Bom Approximation in which the Bom amplitude is 
normalized by the Yukawa couplings, g^-j^,, of Eq. [I9l 



da IB A duB orn 
dt ^ dt 



(24) 



The Improved Bom Approximation incorporates the effective Lagrangian approximation to the SQCD 
effects on the bb4> Yukawa couplings at low energy, but does not include the full SQCD calculation. 
In particular, the "Improved Born Approximation" does not include contributions from box diagrams 
including internal squarks and gluinos or the full momentum dependence of the SQCD contributions. 



12.3 Results 

In Figs. |25]and|26]we compare the results for bg —>■ hhP and bg bH^ at the LHC [93,222]. The curves 
labelled "LO" use CTEQ6L PDFs, as{fJ,R) evaluated at 1-loop, and use the tree level Yukawa couplings. 
The NLO results use CTEQ6M PDFs with the 2-loop evolution of Osi^iR) and af^°{Mz) = 0.118. 
The Yukawa couplings of both the IBA and NLO results are evaluated using the effective Lagrangian of 
Eq. \T9\ The outgoing b quark is required to have prib) > 20 GeV and | rjh |< 2.5. The renormalization 
and factorization scales are set to be M<^/4. The "Improved Born Approximation" (IBA) curves use NLO 
PDFS and the 2— loop evolution of as(/i/j). The b quark mass in the Yukawa couplings is the running 
MS mass evaluated at 2— loops for the NLO and IBA results and at 1— loop for the LO results. Finally, 
the MSSM parameters are evaluated using FeynHiggs to generate an effective Higgs mixing angle and 
radiatively corrected Higgs masses. 

From Fig. [25] we see that for relatively Ught squark and gluino masses, it is important to include 
the exact SQCD contributions and that the Improved Born Approximation is a poor approximation to the 
complete result. In this case, the SQCD contributions significantly reduce the rate. On the other hand, for 
squark and gluino masses on the TeV scale. Fig. |26] demonstrates that the effective Lagrangian approach 
to including the SQCD corrections is extremely accurate. Both Figs. |25] and [26] assume tan /3 = 40. For 
small values of tan (3, the SQCD corrections are insignificant. 

In Fig. [27] we compare the tree level rate for pp —>■ bh^ + jet, with the irreducible background 
from pp ^ bZ + jet at the LHC for tan/3 = 40. We require pxib) and pr(jet) > 20 GeV and 
I %) Vjet \< 2.5. Unfortunately, the rate is quite small. 

13. CHARGED HIGGS BOSONS IN THE MSSM AT CMS: DISCOVERY POTENTIAL 
13.1 Introduction 

Identifying the mechanism of electroweak symmetry breaking will be one of the main goals of the LHC. 
The most popular models are the Higgs mechanism within the Standard Model (SM) and within the 
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Fig. 25: Total cross section for pp ^ b(f) (tf) = h°,H°) at theLHC withpT(&) > 20 GeV and \ rjb \< 2.5. The curve labelled 
"Complete NLO" includes the full set of ©(a^) QCD and SQCD contributions, while the curve labelled "NLO (gluon only)" 
omits the SQCD contributions. The MSSM parameters are rrig = 250 GeV, mi_^ = 250 GeV, and m^^ = 350 GeV. 



Minimal Supersymmetric Standard Model (MSSM) [223-225]. Contrary to the case of the SM, in the 
MSSM two Higgs doublets are required. This results in five physical Higgs bosons instead of the single 
Higgs boson of the SM. These are the light and heavy CP-even Higgs bosons, h and H, the CP-odd 
Higgs boson, A, and the charged Higgs boson, H^. The Higgs sector of the MSSM can be specified 
at lowest order in terms of the gauge couplings, the ratio of the two Higgs vacuum expectation values, 
tan P = V2/V1, and the mass of the CP-odd Higgs boson, Ma (or Mjj± , the mass of the charged 
Higgs boson). Consequently, the masses of the CP'-even neutral Higgs bosons (and the charged Higgs 
boson) are dependent quantities that can be predicted in terms of the Higgs-sector parameters. The same 
applies to the production and decay properties of the MSSM Higgs boson^^ Higgs-phenomenology 
in the MSSM is strongly affected by higher-order corrections, in particular from the sector of the third 
generation quarks and squarks, so that the dependencies on various other MSSM parameters can be 
important. 

The charged Higgs bosons of the MSSM (or a more general Two Higgs Doublet Model (THDM)) 
have been searched at LEP [226], yielding a bound of Mfj± > 80 GeV [227, 228]. The Tevatron 
placed new bounds on the MSSM parameter space from charged Higgs-boson searches at large tan f3 
and low Ma [229]. At the LHC the charged Higgs bosons will be accessible best at large tan /3 up to 
Ma < 800 GeV [98,230,231]. At the ILC, for Mh± < ^/2 a high-precision determination of the 
charged Higgs boson properties will be possible [232-236]. 

The prospective sensitivities at the LHC are usually displayed in terms of the parameters Ma and 
tan /3 (or Mfj± and tan f3) that characterize the MSSM Higgs sector at lowest order. The other MSSM 
parameters are conventionally fixed according to certain benchmark scenarios [200,231]. We focus 
here [237] on the 5 a discovery contours for the charged MSSM Higgs boson for the two cases Mjj± < 
mt and Mh± > rnt, within the m™^^ scenario and the no-mixing scenario. For the interpretation of the 
exclusion bounds and prospective discovery contours in the benchmark scenarios it is important to assess 
how sensitively the results depend on those parameters that have been fixed according to the benchmark 
prescriptions. Consequently, we investigate how the 5 a discovery regions in the Mj|^±-tan /3 plane for 
the charged MSSM Higgs boson obtainable with the CMS experiment at the LHC depend on the other 
MSSM parameters, most prominently the Higgs mixing parameter /i. 



'if the production or decay involves SUSY particles at tree-level, also other MSSM parameters enter the prediction. 
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Fig. 26: Total cross section for pp ^ b(j> ((j> = h°,H°) at theLHC withpT(&) > 20 GeV and \ rjb \< 2.5. The curve labelled 
"Complete NLO" includes the full set of ©(a^) QCD and SQCD contributions, while the curve labelled "NLO (gluon only)" 
omits the SQCD contributions. The MSSM parameters are rrig ~ rrij;^ ~ mj^ ~ 1 TeV. 

13.2 Experimental Analysis 

The main production channels at the LHC are 

pp^tt^ H~b t or t H~^b (25) 

and 

gb H't or gb H^i . (26) 
The decay used in the analysis to detect the charged Higgs boson is 

TUr . (27) 

The analysis described below correspond to CMS experimental sensitivities based on full simulation 
studies, assuming an integrated luminosity of 30 fb~^. 

13.21 The light charged Higgs Boson 

The "light charged Higgs boson" is characterized by Mij± < rrit. The main production channel is given 
in eq. (|25] ). Close to threshold also eq. (l26l ) contributes. The relevant (i.e. detectable) decay channel 
is given by eq. ( [271 ). The experimental analysis, based on 30 fb^^ collected with CMS, is presented in 
Ref. [238]. 

A total number of events leading to final states with the signal characteristics is evaluated, in- 
cluding their respective experimental efficiencies. The various channels and the corresponding effi- 
ciencies can be found in Tab. [13] A 5 cr discovery can be achieved if a parameter point results in 
more than 5260 events (with 30 fb~^). We furthermore used BR(VF=^ Ivi) = 0.217 (/ = /i,e), 
BR{W^ TVr) = 0.1085, BR{W^ jets) = 0.67, BR(r hadrons) = 0.65. The next-to-leading 
order LHC cross section for top quark pairs is taken to be 840 pb. For the W+3 jets background the 
leading order cross section for the process pp + 3 jets, — > ('^u {£ = e, ix) as given by 

MadGraph Ref. [173,239] generator (840 pb) was used. 
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Fig. 27: Comparison of the transverse momentum distributions for the bottom quark jet in the signal process pp ^ b(j) + j and 
background process pp ^ bZ + j at the LHC. 



channel 


exp. efficiency 


pp^tt^ H+b t ^ (r+z^^) (W+b); r hadrons, W Ivi 


0.0052 


pp ^ ti ^ (rur) (li'i)', T hadrons 


0.00217 


pp^ti^ W+ W- {Ivi) {Ivi) 


0.000859 


pp^tt^ W+ W- Get jet) {li^i) 


0.000134 


pp^W + 3 jets, W ^ 


0.000013 



Table 13: Relevant channels for the light charged Higgs bosons and their respective experimental efficiencies. The charge 
conjugated processes ought to be included. The efficiency for the charged Higgs production is given for A/^+ = 160 GeV. / 
denotes e or fi. 

13.22 The heavy charged Higgs Boson 

The "heavy charged Higgs boson" is characterized by Mfj± > nit. Here eq. ( |26l ) gives the largest 
contribution, and very close to threshold eq. (1251 ) can contribute somewhat. The relevant decay channel 
is again given in eq. ( [27] ). The experimental analysis, based on 30 fb^^ collected with CMS, is presented 
in Ref. [240]. 

The number of signal events is evaluated as 

iVev = X a{pp + X)x BR{H^ rvr) x BR(r hadrons) x exp. eff. , (28) 

where C denotes the luminosity and the experimental efficiency is given in Tab. [I4]as a function of Mfj± . 
Aba discovery corresponds to a number of signal events larger than 14.1. 



Mh± [GeV] 


171.6 


180.4 


201.0 


300.9 


400.7 


600.8 


exp. eff. [IQ-^] 


3.5 


4.0 


5.0 


23 


32 


42 



Table 14: Experimental efficiencies for the heavy charged Higgs boson detection. 



The charged Higgs boson production with the mass close to the top quark mass (first column in 




Tab-O was generated with the PYTHIA generator processes 401 {gg tbH^) and 402 {qq — > tbH^) 
implemented as described in Ref. [241]. 



13.3 Calculation of Cross Section and Branching Ratios 

For the calculation of cross sections and branching ratios we use a combination of up-to-date theory 
evaluations. The Lagrangian for the interaction of the charged Higgs boson with the t/b doublet is given 
by [217] 

g mb 



2Mw 1 + Aft 



V2Vtb tan (3 H+tLbR 



+ h.c. (29) 



Here denotes the running bottom quark mass including SM QCD corrections. The prefactor 1/(1 + 
Afe) in eq. (l29l ) arises from the resummation of the leading corrections to all orders. The explicit form of 
Aft in the limit of heavy SUSY masses and tan /3 S> 1 reads [218,220,242] 

Ab = nig fitanP x /(m^^ , m^^ , rrig) + — At tan (3 x /(m^-^ , m^^ , /x) . (30) 

Here m^^, m^^, m^^, m^_^ denote the i and b masses, rrig is the gluino mass. Large negative fi can lead to 
a strong enhancement of the H^tb coupling, while a large positive value leads to a strong suppression. 
Concerning the m™'^^ and the no-mixing benchmark scenarios, as discussed in Refs. [231,243] the A^ 
effects are much more pronounced in the m™^^ scenario, where the two terms in eq. (|30l ) are of similar 
size. In the no-mixing scenario the first term in eq. (l30l ) dominates, and the total effect is smaller. 

For the production cross section in eq. (l25l) we use the SM cross section a{pp — > it) = 840 pb 
times the BR(i/^ — > tb) including the A;, corrections described above. The production cross section 
in eq. (l26l ) is evaluated as given in Refs. [244, 245] . In addition also the Af, corrections of eq. (l29l ) are 
applied. Finally the BR(ff^ — > rv-r) is evaluated taking into account all decay channels, among whom 
the most relevant are — > tb, cs, W^*^h. For the decay to tb again the Ab corrections are included. 
All the numerical evaluations are performed with the program FeynHiggs [246-249]. 



13.4 Numerical Analysis 

The numerical analysis has been performed in the m^*^^ and the no-mixing scenario [200, 231] for 
^ = -1000, -200, +200, +1000 GeV. In Fig.[28]we show the results combined for the 5 cr discovery 
contours for the light and the heavy charged Higgs boson, corresponding to the experimental analyses 
in Sects. [TT2T] and [T3.22[ respectively. As described above, the analyses were performed for the CMS 
detector and 30 fb~^. The top quark mass is set to mj = 175 GeV. 

Within the m™^^ scenario, shown in the left plot of Fig.|28]the search for the light charged Higgs 
boson covers the area of large tan f3 and Mf^± < 150 GeV. The variation with /i induces a strong shift 
in the 5 a discovery contours with A tan /? = 15 for Mh± = 100 GeV, rising up to A tan /3 = 40 for 
larger Mfj± values. The discovery region is largest (smallest) for = — (+)1000 GeV, corresponding 
to the largest (smallest) production cross section. 

The 5 a discovery regions for the search for heavy charged Higgs bosons show a similar behavior. 
For Mh± = 170 GeV the accessible parameter space starts at tan jS = 20(58) for fi = -(+)1000 GeV, 
i.e. the variation of jj, again induces a very strong shift in the 5cr discovery contours. For Mfj± = 
300 GeV the 5 a regions vary from tan /3 = 38 to tan /3 = 54. For = —1000 GeV and larger tan /? 
values the bottom Yukawa coupling becomes so large that a perturbative treatment would no longer be 
reliable in this region. 

The no-mixing scenario is shown in the right plot of Fig. |28] The qualitative behavior is the same 
as for the m™**^ scenario. However, as discussed above, the effects from a variation of /i are much less 
pronounced. The induced shifts stay below A tan /3 = 20(10) in the search for the light (heavy) charged 
Higgs boson. The 5 a discovery areas are slightly larger than in the m™'^^ scenario. 




Fig. 28: Discovery reach for tlie cliarged Higgs boson of CMS witli 30 fb in tiie Afjif±-tan /3 plane for tiie m™^" scenario 
(left) and the no-mixing scenario (right). 



13.5 Conclusions 

We have presented the 5 a discovery contours for the search for the charged MSSM Higgs boson. We 
combine the latest results for the CMS experimental sensitivities based on full simulation studies with 
state-of-the-art theoretical predictions of MSSM Higgs-boson properties. The experimental analyses are 
done assuming an integrated luminosity of 30 fb^^ for the two cases, Mu± < rrit and Mfj± > rrif. 

The numerical analysis has been performed in the m™^^ and the no-mixing scenario for fj, = 
±200, ±1000 GeV. The search for the light charged Higgs boson covers the are area of large tan /3 and 
Mh± < 160 GeV. The search for the heavy charged Higgs boson reaches up to Mh± ^ 400 GeV for 
large tb. The variation of /i induces a very strong shift in the 5 a discovery contours of up to A tan /? = 
40. The effect enters via the variation of A;,, affecting the charged Higgs production cross section and 
branching ratios. Large negative fi values give the largest reach, while large positive values yield the 
smallest 5 a discovery areas. 

Acknowledgements 

The work of S.H. was partially supported by CICYT (grant FPA 2007-66387). Work supported in part 
by the European Community's Marie-Curie Research Training Network under contract MRTN-CT-2006- 
035505 'Tools and Precision Calculations for Physics Discoveries at Colliders'. 

14. STUDIES OF SPIN EFFECTS IN CHARGED HIGGS BOSON PRODUCTION WITH AN 
ITERATIVE DISCRIMINANT ANALYSIS AT THE TEVATRON AND LHcEll 

14.1 Introduction 

The importance of charged Higgs boson searches has in the recent years been emphasized, including 
in the '2005 Les Houches' proceedings [90]. This work extends the charged Higgs boson '2005 Les 
Houches' studies. It is the purpose of this note to outline the possible improvements that can be achieved 
at the Tevatron and LHC in the search for charged Higgs bosons focussing on the spin effects and 
the —>■ Ti>r decay. In order to quantify the spin effect an Iterative Discriminant Analysis (IDA) 
method [250] has been applied, which is a powerful tool to separate signal and background, even in 
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cases such as the one presently under study when several selection variables with limited discriminant 
power are present. 

14.2 Tevatron Energy 

We start by studying charged Higgs production qq, gg — > thH^ with subsequent decays t — > bW , — > 
Ti/r at the FNAL Tevatron with ^/s = 1.96 TeV. In the following we analyze hadronic decays of the 
boson and r lepton {W^ —>■ qq', r —>■ hadrons+i/T-), which results in the signature 2b+2j +Tjct+pT^^^ (2 
6 jets, 2 light jets, 1 r jet and missing transverse momentum). The most important irreducible background 
process is qq,gg — > tt with the subsequent decays t bW^ and t — > bW~, one boson decaying 
hadronically (W^ —>■ qq') and one leptonically (W^ —>■ Ti>r), which results in the same final state 
particles as for the expected signal. 

14.21 Simulation and Detector Response 

The signal process qq,gg tbH^ is simulated with PYTHIA [251]. The subsequent decays t — > bW^ 
(or its charge conjugate), — > qq' and — > are also carried out within PYTHIA, whereas the 
r leptons are decayed externally with the program TAUOLA [252, 253], which includes the complete 
spin structure of the r decay. The background process qq, gg tt is also simulated with PYTHIA with 
the built-in subroutines for ti production. The decays of the top quarks and bosons are performed 
within PYTHIA and that of the r lepton within TAUOLA. 

The momenta of the final b and light quarks from the PYTHIA event record are taken as the mo- 
menta of the corresponding jet, whereas for the r jet the sum of all non-lep tonic final state particles as 
given by TAUOLA is used. The energy resolution of the detector and parton shower and hadronization 
effects are emulated through a Gaussian smearing {A{pt)/pt)'^ = of the transverse mo- 

mentum for all jets in the final state, including the r jet [254]. As typical for fast simulation studies, 
no effects of underlying events, are simulated. Events are removed which contain jets with less than 20 
GeV transverse momentun J^. corresponding to about |ry| > 3. The transverse momentum of the lead- 
ing charged pion in the r jet is assumed to be measured in the tracker independently of the transverse 
momentum of the r jet. The identification and momentum measurement of the pion is important to fully 
exploit the r spin information. In order to take into account the tracker performance we apply Gaussian 
smearing on 1/pJ with a{l/pJ)[TeV''^] = y/0.522 + 222/(pnGeV])2 sin6'^, where 9^ is the polar an- 
gle of the IT. The missing transverse momentum pf^^^^ is constructed from the transverse momenta of all 
visible jets (including visible r decay products) after taking the modelling of the detector into account. 
The generic detector description is a good approximation for both Tevatron experiments, CDF and DO. 

1 4. 22 Expected Rates 

For completeness we present a brief discussion of the expected cross section of the charged Higgs boson 
signature under investigation. The signal cross section has been calculated for tan /3 = 30 and mij± = 
80, 100, 130 and 150 GeV with PYTHIA, version 6.325, using the implementation described in [241], in 
order to take the effects in the transition region into account. Furthermore, it has been shown in [255] that 
the signal cross section for tbH^ agrees with the one from the top-decay approximation tt — > tbH^ for 
charged Higgs boson masses up to about 160 GeV if the same factorization and renormaUzation scales 
are used. Thus, we have used everywhere in this study the factorization scale {mt + mu± ) /4 and the 
renormalization scale mjj± for both signal and background (i.e., those recommended in [241] as most 
appropriate for the tbH^ signaljM since the primary purpose of our study is to single out variables that 
show a difference between our and data samples and that this can unambiguously be ascribed to 

^"In order to be largely independent of the specific detector performance, no requirement on the jet resolution is applied. 
^ Clearly, for a proper experimental study, factorization and renormalization scales for our background process qq, gg 
ti — > tbW^ ought to be chosen appropriately, i.e., unrelated to the charged Higgs boson mass. 



the different nature of the two kinds of bosons (chiefly, their different mass and spin state). In addition, 
the running b quark mass entering in the Yukawa coupUng of the signal has been evaluated at mfj± . This 
procedure eventually results in a dependence of our background calculations on tan (3 and, especially, 
that is more marked than the one that would more naturally arise as only due to indirect effects 
through the top decay width. Hence, the cross sections have been rescaled with a common factor such 
that the total tt cross section is (J^^°'^ = 5.2 pb [256]. To be more specific, we have first calculated 
the total cross section (jP^^^'PYTHIA •) ^^^q built-in routine for tt production in PYTHIA for 
all rriHi = 80, 100, 130 and 150 GeV and then calculated from this the respective rescaling factors 
c{mu±) = 5.2 pb/(T^j°'^'^^^^^^(mj:^±) for each ■mu±. Then we have calculated the background cross 
section for m}j± = 80 GeV into the final state with the signature 26 + 2j + rjet + p™'^^ by enforcing 
the respective decay channels in PYTHIA using the built-in routine for ti production and multiplied it 
with c(80 GeV). In the same manner we have calculated the signal cross sections with the PYTHIA 
routines for tbH"^ production by enforcing the respective decay channels in PYTHIA and multiplying 
with the rescaling factors c{mu±) for mu± = 80, 100, 130, 150 GeV. The resulting cross sections are 
given in Table [T5] before (cr*^) and after (a) applying the basic cuts j^^^^ > 20 GeV and the hard cut 
ptnms ^ j^QQ Pqj- tijg fQuj- signal masses, the tbH^ and ti — > tbH^ cross section calculations agree 
numerically. 



Table 15: Tevatron cross sections of background qq, gg tt and signal qq, gg —> tbH for tan = 30 and m^± = 
80, 100, 130 and 150 GeV into the final state 26 + 2j + rjct + vT'"" before (cr"') and after (cr) the basic cuts (pt > 20 GeV for 
all jets) and the hard cut (pr'"" > 100 GeV). 





gg, gg tt 


g 


q^gg - 


tbH 


± 


mH± (GeV) 


80 


80 


100 


130 


150 


cj*'^ (fb) 


350 


535 


415 


213 


85 


a (fb) for pJ/*' > 20 GeV 


125 


244 


202 


105 


32 


a (fb) for {j^t^^pT^^) > (20, 100) GeV 


21 


30 


25 


18 


7 



14.23 Event Preselection and Discussion of Discriminant Variables 

The expected cross sections of the 2b + 2j + rjet + p^^^^ signature are of the same order of magnitude 
for the signal and background reactions, as shown in Table [151 Thus, the same number of signal and 
background events is assumed for the analysis of different kinematic selection variables. For the signal 
5-10^ events have been simulated with PYTHIA for each charged Higgs mass at the Tevatron energy 
of 1 .96 TeV using the built-in tt routine in the ti — > tbH^ approximation, while for the ti background 
also 5 • 10^ events have been simulated using the built-in routine. Then the basic cuts p'^'^^ > 20 GeV 
are applied. An additional hard cut on the missing transverse momentum p^^^^ > 100 GeV is used to 
suppress the QCD background, as for example demonstrated in Ref. [257]. After the additional anti-QCD 
cut about 28000 to 42000 signal events, depending on the simulated charged Higgs bosons mass, and 
about 30000 ti background events remain. Other background reactions, for example W-i-jet production, 
are expected to be negligible because they have either a much lower production cross section or are 
strongly suppressed compared to ti background, as quantified for example in Ref. [257]. In addition to 
the previous study (based on 5000 x BR(r hadrons) events each) [90], the present one applies an IDA 
method [250] to explore efficiencies and purities. As already mentioned, particular attention is devoted 
to the study of spin sensitive variables in the exploitation of polarization effects for the separation of 
signal and background events. 

Examples of the signal and background distributions of some of the kinematic variables used in 
the IDA method and the respective difference between signal and background distributions are given in 
Ref. [258], namely: 

• the transverse momentum of the r jet, p]-^"^ , 



• the transverse momentum of the leading vr in the r jet, , 

• the ratio p^^ /pl'"^ , 

• the transverse momentum of the second (least energetic) b quark jet, , 

• the transverse masS in the rjet + pf^^^^ system, mt = 2p['°^ pf-^^^\l — cos(A(/))], where is 
the azimuthal angle between p^"^ and p™'^^, 

• the invariant mass distribution of the two light quark jets and the second h quark jet, rnjjh^, 

• the spatial distance between the rjet and the second b quark jet, AR{t, 62) = \/(A^p~+~(A?7p^, 
where Acp is the azimuthal angle between the r and 6 jet, and 

• the sum of the (scalar) transverse momenta of all the quark jets, i/jets = Pt^ + pi^ + Pt^ + Pt^- 

The distributions of signal and background events are normalized to the same number of 10^ events, in 
order to make small differences better visible. 

The signal and background distributions for the variables shown in Ref. [258] are as expected 
rather similar for mfj± = my/± and are hence mostly important to discriminate between signal and 
background in the IDA for > m^yi. Especially the transverse mass, shows a large variation 

with the charged Higgs boson mass. However, the different spin of the charged Higgs boson and the 
boson has a large effect on the r jet variables and p^ resulting in significantly different 
distributions of signal and background even for mu± = m^r± . Moreover, the spin effects in the p['°^ and 
p^ distributions are correlated where the distributions of the ratio pj^ /pl'"^ [259-261] show even larger 
differences [258]. This highlights the importance of the additional variable p^ (and hence p^ /Pt"^), 
compared to a previous study [90]. The large separation power of this variable is indeed due to the 
different r polarizations in signal and background [258]. There the signal and background distributions 
for pp"', Pt and pj are shown for reference samples where the r decay has been performed 

without the inclusion of spin effects with the built-in routines of PYTHIA and hence the differences 
between signal and background nearly vanish. 

14. 24 Iterative Discriminant Analysis ( IDA ) 

The IDA method is a modified Fisher Discriminant Analysis [250] and is characterized by the use of a 
quadratic, instead of a linear, discriminant function and also involves iterations in order to enhance the 
separation between signal and background. 

In order to analyze our events with the IDA method, signal and background have been split in two 
samples of equal size. With the first set of samples the IDA training has been performed and then the 
second set of samples has been analyzed. We have used the following 20 variables in the IDA study: 
the transverse momenta pi , pf^^'^, Pi, pj-^, pj-^, pj--'; the transverse mass mt, the invariant 
masses rrijj, rajji,^, mjjiy^, and s = mjjbbr', the spatial distances Ai?(r, 61), Ai?(r, 62), AR{T,ji), 
Ai?(r, ^2); the total transverse momenta of all quark jets ^/jets and of all jets ifaii = -f^jcts + In 
the analysis of real data, b-quark tagging probabilities and the reconstruction of t and W masses could 
be used to improve the jet pairing, and replace the allocation of least and most energetic 6-jet by a 
probabilistic analysis. 

The results of the IDA study are obtained for the event samples with spin effect in the r decays for 
m}j± = 80, 100, 130, 150 GeV and for the reference samples without the spin effect for mj:^± = 80 GeV 
in order to illustrate the spin effect. In all plots of the IDA output variable the number of background 
events has been normalized to the number of signal events. Two IDA steps have been performed. After 
the first step, 90% of the signal is retained when a cut at zero is applied on the IDA output variable. 
The signal and background events after this cut are then passed to the second IDA step. A cut on IDA 
output variable distributions after the second step leads to the efficiency and purity (defined as ratio 



Strictly speaking this is not the transverse mass since there are two neutrinos in the decay chain of the charged Higgs boson 
we are considering, even so the characteristics of this mass are very similar to that of the true transverse mass. 



of the number of signal events divided by the sum of signal and background events) combinations. 
These combinations define the working point (number of expected background events for a given signal 
efficiency) and the latter can be optimized to maximize the discovery potential. 

In order to illustrate the effect of the hard cut on the missing transverse momentum (p™'*^*^ > 
100 GeV), which is imposed to suppress the QCD background, the final efficiency-purity plot of the IDA 
analysis is shown in Fig.|29]for mfj± = 80 GeV for two reference samples (red, long dashed: with spin 
effects in the r decay; red, dotted: without spin effects) without imposing the hard cut. As expected the 
achievable purity for a given efficiency decreases with the hard cut, therefore the spin effects become even 
more important to separate signal and background. In principle, by choosing the signal reduction rates 
in the previous IDA iterations, the signal and background rates in the final distributions can be varied 
appropriately. However, we have checked that a different number of IDA iterations and/or different 
efficiencies for the first IDA iteration have only a minor effect on the final result. 




Fig. 29: Efficiency as a function of purity for 
m^i = 80 GeV and ^ = 1.96 TeV. Tlie 
black lines are the results after applying the hard 
cut p™"""" > 100 GeV when not taking the spin 
effects in the r decay into account (dashed) and 
with spin effects in the r decay (solid). The red 
lines are the results without applying the hard cut 
on p™'°° when not taking the spin effects in the r 
decay into account (dotted) and with spin effects 
in the r decay (long dashed). 
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14.3 LHC Energy 

The simulation procedure and the emulation of the detector response are the same as those outlined in 
Sect. 2.1 for the Tevatron, as well as, for the preselection and IDA method, as described in Sects. 2.3 
and 2.4, respectively. Hence, only the expected LHC rates are discussed, followed by the description of 
changes in the distributions of the variables and the final IDA results. 

Unlike the case of the Tevatron, where only charged Higgs masses smaller than the top quark 
mass can be explored, and 2HDM/MSSM signatures practically rely on rur pairs only, at the LHC the 
phenomenology is more varied. Here, the search strategies depend strongly on the charged Higgs boson 
mass. If < mt — rrn, (later referred to as a light Higgs boson), the charged Higgs boson can 

be produced in top (anti-)quark decay. The main source of top (anti-)quarks at the LHC is again tt 
pair production (cr^f = 850 pb at NLO) [262]. For the whole (tan/?, m^/i) parameter space there is 
a competition between the bW^ and hH^ channels in top decay keeping the sum BR(t bW^) + 
BR(t — > bH^) at almost unity. The top quark decay to bW^ is however the dominant mode for most 
of the parameter space. Thus, the best way to search for a (light) charged Higgs boson is by requiring 
that the top quark produced in the tbH^ process decays to a W^. While in the case of decays r's 
will be tagged via their hadronic decay producing low-multiplicity narrow jets in the detector, there are 
two different decays that can be explored. The leptonic signature bbH^W^ bbrvlv provides 
a clean selection of the signal via the identification of the lepton / = e, In this case the charged 
Higgs transverse mass cannot be reconstructed because of the presence of two neutrinos with different 
origin. In this channel charged Higgs discovery will be determined by the observation of an excess of 
such events over SM expectations through a simple counting experiment. In the case of hadronic decays 
bbH^W^ — > bbrvji the transverse mass can instead be reconstructed since all neutrinos are arising 
from the charged Higgs boson decay. This allows for an efficient separation of the signal and the main 
tt — ^ bbW^W^ — > bbriyjj background (assuming 'mfj± > m\Y±). The absence of a lepton (e or /x) 



provides a less clean environment but the use of the transverse mass makes it possible to reach the same 
mass discovery region as in the previous case and also to extract the charged Higgs boson mass. Both 
these channels show that after an integrated luminosity of 30 fb^^ the discovery could be possible up to 
a mass of 150 GeV for all tan/3 values in both ATLAS and CMS [263-265]. 

If the charged Higgs is heavier than the top quark, the dominant decay channels are TV 
and tb depending on tan/3. They have both been studied by ATLAS and CMS [266-269]. The 

charged Higgs bosons are produced in the pp tbH^ channel. For the — > tb decay, a charged 
Higgs boson can be discovered up to high masses {mH± ~ 400 GeV) in the case of very large tan (3 
values and this reach cannot be much improved because of the large multi-jet environment. For the 
— > TV decay mode this reach is larger due to a cleaner signal despite a lower BR. In this case the 5a 
reach ranges from tan /3 = 20 for mfj± = 200 GeV to tan /? = 30 for mfj± = 400 GeV. 

For the LHC, signal and background events have been simulated in the same way as for the Teva- 
tron as described before, however, without implying any rescaling factor to match a measured ti cross 
section. Table [16] lists the resulting cross sections before (a^^) and after (a) applying the basic cuts 
p^'^'' > 20 GeV and the hard cut p™'^^ > 100 GeV. The LHC rates allow for the discovery to be less 
challenging than at the Tevatron in the region m/^± ~ m\Y±, yet the separation of signal events from 
background remains crucial for the measurement of the charged Higgs mass. 

Table 16: LHC cross sections of background qq, gg tt and signal qq, gg tbH^ for tan /? = 30 and mjj± — 80, 100, 130 
and 150 GeV into the final state 2b + 2j + rjct + pT''''' before (cr"') and after (a) the basic cuts (pt > 20 GeV for all jets) and 
the hard cut (p^'"" > 100 GeV). 
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The kinematic distributions for ^/s = 14 TeV are shown in Ref. [258]. The choice of variables 
is identical to the one for the Tevatron and allows for a one-to-one comparison, the differences being 
due to a change in CM energy (and, to a somewhat lesser extent, due to the leading partonic mode of 
the production procesQ). The main differences with respect to the Tevatron case are that the various 
transverse momenta and invariant masses have longer high energy tails. In particular, it should be noted 
that the effect of the spin differences between and events can be explored very effectively also 
at LHC energies, e.g. the ratio which is very sensitive to the spin effects. These observations 

lead to the conclusion that the same method using spin differences can be used to separate signal from 
background at both the Tevatron and the LHC. 

The distributions of the IDA output variables are shown in Ref. [258] for the study at ^/s = 
14 TeV for two steps with 90% efficiency in the first step. These distributions are qualitatively similar 
to those for the Tevatron The final achievable purity for a given efficiency is shown in Fig. [30l As 
for the Tevatron energy a good separation of signal and background events can be achieved with the 
spin sensitive variables and the IDA method even in case 'mfj± ~ m^r±. For heavier masses 
the separation of signal and background events increases due to the kinematic differences of the event 
topology. 

14.4 Conclusions 

The discovery of charged Higgs bosons would be a clear sign of physics beyond the SM. In this case 
study we have investigated charged Higgs boson topologies produced at the current Tevatron and LHC 
energies and compared them against the irreducible SM background due to top-antitop production and 

As the latter is dominated by qq annihilation at the Tevatron and gg fusion at the LHC. 



Fig. 30: Efficiency as a function of the purity 
when not taking the spin effects in the r decay 
into account for m^fi — 80 GeV (dashed) and 
with spin effects in the r decay for m^± — 
80, 100, 130, 150 GeV (solid, from left to right). 
Results are for the LHC. 
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decay. While sizable differences between signal and background are expected whenever mjj± ^ m^y±, 
near the current mass limit of about mii± « 80 GeV the kinematic spectra are very similar between 
SM decays and those involving charged Higgs bosons. In this case, spin information will significantly 
distinguish between signal and irreducible SM background. In fact, we have considered hadronic rv-j- 
decays of charged Higgs bosons, wherein the r polarization induced by a decaying (pseudo)scalar object 
is significantly different from those emerging in the vector (M^^) decays onsetting in the top-antitop 
case. For a realistic analysis which is not specific for a particular detector, a dedicated Monte Carlo 
event generation and a simplified multipurpose detector response approximation have been applied. The 
identification of a hadronic tau-lepton will be an experimental challenge in an environment with typically 
four jets being present. We have demonstrated how an IDA method can be an applied to separate signal 
and background when the differences between the signal and background distributions are small. Our 
results show that the IDA method will be equally effective at both the Tevatron and LHC. While only 
the dominant irreducible it background has been dealt with in detail, we have also specifically addressed 
the QCD background. A suitably hard missing transverse momentum cut has been applied to reject 
such jet activity and we have demonstrated that although the discriminative power is reduced by such 
a cut, the reduction is small compared to the gain from including the r polarization effects. Using the 
differences in r polarization between the signal and the dominant SM irreducible ti background is crucial 
for disentangling the former from the latter. 
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Part IV 

CP VIOLATING HIGGS BOSONS 

15. JET ASSIGNMENT STUDIES IN THE SEARCH FOR THE DECAY t bH+, H+ 
HfW+, Hi bb IN THE CPX MSSM SCENARIO^ 

15.1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) can have loop-induced CP-violation (CPX) if 
the Higgsino mass parameter, the gaugino masses and the trilinear couplings are complex. One of the key 
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features of the CPX scenario is the suppression of the couphngs of the neutral Higgs boson to both vector 
boson pairs and to ti pairs. The suppression of the HfVV coupUng effectively dilutes the limits set on 
the neutral Higgs using LEP data [270], allowing the existence of a light neutral Higgs boson (40 - 50 
GeV) and a relatively light charged Higgs boson {M{H^) < Mtop) at low tan (3. The suppression of the 
couplings also makes the usual search methods at hadron colliders unviable. However, the suppression 
of the H'IVV leads to the enhancement of the HiH^W^ coupling via a sum rule, making tt production 
events in which one of the top quarks decays via t bH^, — > HfW^, Hf — > bb one of the most 
promising search channels for the CPX scenario [271]. 

We present here a study of mass reconstruction and the impact of jet misassignment on this search 
using the CMS detector; A feasibility study for discovering the Higgs bosons in the CPX scenario also 
using the CMS detector is presented in Section [THH 

15.2 Event generation 

The signal event sample was generated using PYTHIA [183] and assuming the following parameters: 
M{H'^ = 51 GeV, M{H^) = 133 GeV, Mt = 175 GeV, tan /? = 5 and ^>cp = 90°. In each event, 
one of the top quarks was forced to decay in the usual way, ie. t — > bW, while the other was forced to 
decay via t — > bH^, H^W^, bb. All possible decays were allowed. The relevant 

branching fractions were calculated using CPSuperH [272] and were found to be: BR{t —>■ bH^) = 
0.01, BR(H+ HlW+) = 0.99 and BR(Hl bb) = 0.92. Taking the total ti production cross 
section to be 840 pb [273], this gives a cross section for this process of 8.68 pb. 

For the purposes of this study only the subset of signal events in which one decayed hadroni- 
cally and the other decayed leptonically (electron or muon) were considered, as this is the experimental 
signature that will be used to identify events in this analysis. 

15.3 Event selection and mass reconstruction 

This study was performed using only generator-level information. The iterative cone (IC) algorithm [274] 
with a radius of 0.5 was used for jet identification. The jets are formed out of stable generator-level 
particles, although neutrinos and muons are explcitly excluded from the process. Six or more jets must 
be found using the IC algorithm that satisfy the following requirements: p!^^ > 20 GeV and \r]jet\ < 2.4. 
Three of more must also satisfy p^^* > 30 GeV. In addition, an electron or muon that satisfies pip > 20 
GeV and \r]i\ < 2.4 must also be present and the missing Et reconstructed from generator-level particles 
must be greater than 20 GeV. 

Events that pass these selection requirements then undergo the mass reconstruction procedure. 
The events are searched for the two possible decay channels, namely, t bqq', t — > bbblv and t blv, 
t bbbqq' + (c.c). As the study presented here is performed using generator-level information, the 
true lepton and neutrino from the leptonically-decaying are used. This means that during the mass 
reconstruction procedure, the four-vector is calculated simply by summing the lepton and neutrino 
four- vectors. When the decays hadronically, the mass is reconstructed using jets and must lie within 
25 GeV of the nominal value. The corresponding mass constraint is placed on all reconstructed top 
masses. When reconstructing both of the top masses from a given jet combination, the jet associated 
with the b-quark {t — > bW^ or t — > bH^) must satisfy pT > 40 GeV. A number of jet combinations will 
pass these requirements in each event and therefore the best candidate for a given event is selected by the 
minimization of a based on the top masses and the mass of the hadronically-decaying candidate. 

It should be noted that this mass reconstruction procedure results in three possible jet combinations 
associated with the best candidate x^- This is because the jets associated to the three b-quarks produced 
in the t bbbW^ decay can be swapped around, but still give the same top mass value. However, the 
stricter jet px requirement applied to the jet associated to the b-quark from the t —>■ bH^ decay can cause 
one or possibly two of the three combinations to be rejected before the calculation is performed. All 
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Fig. 31: The top mass distribution from the decay t —* bqq' and the M{H1) distribution from the decay t hbhlv re- 
constructed from angular-matched jets. All distributions are made with angular-matched jets that satisfy Ai? < 0.5. The 
dashed histograms, in addition, have the top and mass constraints applied, while the solid lines have the Ai? requirements 
tightened on for the decay products of the qq and the Hi bb. 



the combinations corresponding to the best candidate that also satisfy the stricter jet pr requirement 
will be used when making the mass distributions. 

15.4 Mass reconstruction studies 

Before attempting to reconstruct masses at the detector level, it is important to understand whether good 
mass reconstruction is possible. This is done by identifying the jets associated to the quarks produced in 
the decay channel (these quarks are hereafter referred to collectively as "decay quarks") and reconstruct- 
ing the masses from these jets. 

The association of jets with the decay quarks is done using two possible matching procedures: 
Angular matching, in which the quantity Ai? = Ar^^ -\- A(/)2 is used to determine a unique set of jet- 
parton matches; or jet constituent matching, in which the particles assigned to a given jet are classified 
according to the top quark decay from which they originated. The fraction of the transverse momentum 
of a given jet, p!^^ , carried by the constituents originating from each decay quark can then be determined 
and used to create a unique set of jet-parton matches. 

Figure [3T] shows the top mass distribution from the t — > bqq' decay and the mass distribution 
from the t bhblv decay reconstructed using angular-matched jets. The points correspond to those 
made using only jets that satisfy Ai? < 0.5 and it can be seen that in both cases a clear peak is visible in 
the correct position, although the mass has a noticeable high mass tail. The dashed lines represent the 
distributions after some mass constraints have been applied: in the case of the top mass from the t bqq' 
decay, the light-quark jet pair must give mass within 25 GeV of the nominal value, while the mass 
from the corresponding t —>■ bbblv decay must lie within 25 GeV of the nominal value. In the case of the 
mass distribution both top masses and the hadronically-decaying must lie within 25 GeV of their 
nominal values. These mass constraints reduce slightly the high mass tail on the mass distribution. 
The solid lines do not have the mass constraints applied, but instead the Ai? requirement on the decay 
products of the and the hadronically-decaying boson have been tightened to Ai? < 0.1. This all 
but removes the high mass tail on the neutral Higgs mass distribution, suggesting that the tail is caused 
by problems in the jet-parton matching procedure. 

Given the large number of jets in these events, the most likely reason for having problems with 
jet-parton matching (and potentially more generally with mass reconstruction) is that the jets tend to 
overlap with each other This can be verified using the p^^* fractions used for jet constituent matching. 
These fractions are determined by tracing all the particles associated to a given jet back to the top quark 



jet p.^ fraction from ligfit quarks 



jet fraction from light quarks 



Fig. 32: (a) A comparison of the jet transverse momentum fractions for constituents coming from tlie b-quarlc and tire ligiit 
quarlcs in tiie decay t —> bqq' and (b) A comparison of the jet transverse momentum fractions for consituents coming from the 
decay products of the neutral Higgs and the hadronically-decaying in the decay channel t bbbqq . 

decay they came from. The transverse momenta of the particles associated to a given decay quark are 
then summed and the result divided by the jet transverse momentum, resulting in six pT fraction values 
per jet. 

Figure fy2\ a) compares the p^^* fractions for all jets with particles associated to the b-quark and 
to either of the light quarks in the decay t bqq'. Figure l32l b) shows the equivalent distribution for 
the decay channel t — > bbbqq', but compares the the fractions for all jets with particles associated to the 
Hf decay products and the light quarks coming from the decay. No jet angular matching has been 
applied. The two combinations are chosen because they represent the jets from decay quarks that are 
expected to be closest to each other. In the case of the SM top decay, the distribution shows that the jets 
are well separated, as the values are concentrated at very high or low values. In the case of the t bbbqq' 
decay, it is clear that the jets overlap significantly, as suspected. 

15.41 Jet assignment studies 

Although jet overlapping has been identified as a potential problem for mass reconstruction, the results 
in section ITSAl show that it is basically possible to reconstruct reasonable mass distributions. However, 
the impact of jet misassignment on the mass distributions must also be understood and ways found to 
minimize its effect. Jet misassignment arises from two different sources: the misassignment of jets as- 
sociated to the decay quarks and the misassignment of jets associated to other hard partons in the event, 
for example, gluons from initial state radiation or produced during parton showering. The contribution 
from these two sources can be studied by comparing the mass distributions from three different recon- 
struction procedures: those produced using jets matched to the decay quarks ("fully-matched"), those 
produced using the subset of jets matched to the decay quarks, but without using the knowledge about 
which jet belongs to which quark, ("partially-matched") and those produced using the standard mass 
reconstruction procedure ("unmatched"). Comparisons of "fully-matched" and "partially-matched" dis- 
tributions provide information about the misassignment of jets from decay quarks, while comparisons of 
"partially-matched" and "unmatched distributions" provide information about the misassignment of jets 
from other hard partons. 

Figure [33] shows the comparison of these three reconstruction methods for the top mass distribu- 
tion from the t bqq' decay and the mass distribution and the corresponding top mass distribution 
from the t — > bbblv decay. The three methods for the top mass from the SM top decay are broadly in 
agreement, indicating that the reconstruction procedure is working well. The mass distribution shows 
differences between all three methods, indicating that there are contributions from both sources of misas- 
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Fig. 33: A comparison of thie fully-matched, partially-matched and unmatched (with and without b-tagging) mass distributions: 
(a) The top mass distribution from the decay t bqq', (b) the Hi mass distribution from the decay t hhhlv and (c) the top 
mass distribution from the t — > bbblv decay. 



signed jets. However, in the case of the top mass distribution from the same decay the only difference is 
between the partially-matched and unmatched versions, ie. the contribution from the misassignment of 
jets associated to the decay quarks has disappeared. This indicates that the misassigned decay-quark jets 
observed in the mass distribution come from within the t — > bbblv decay chain. The high mass tail 
observed on the unmatched top mass distribution is therefore partially caused by the misassignment of 
jets associated to other hard partons. The remainder of the high mass tail, ie. the contribution that is also 
observed in the fully-matched distribution, is caused by overlapping jets, as discussed in section [TT4l 

One possible method of improving the jet assignment during the mass reconstruction procedure is 
to use b-tagging. To study what impact it may have at generator level, "perfect" b-tagging can be used. 
Perfect b-tagging means using only jets that have been matched to one of the b-quarks if a b-tagged jet 
is required, while only jets not matched to a b-quark are used when a light-quark jet is required. Perfect 
b-tagging has been applied to the unmatched distributions, as shown in Figl33l and it can be seen that the 
differences between the partially-matched and unmatched distributions are eliminated for all the mass 
distributions. This is consistent with the conclusion that this difference is a result of the misassignment 
of jets associated to other hard partons in the event, as the other hard partons are more likely to be gluon 
or light-quark jets than b-quark jets. 

15.5 Conclusions 

A study of jet reconstruction and assignment has been performed at generator level for the analysis 
of CP- violating Higgs production at LHC via the decay channel pp —>■ ttX, t — > bW^, t — > bH^, 
HfW^, Hi bb. It has been established that it is possible to reconstruct reasonable mass 
distributions for this decay channel, but studies of jet-parton matching show that overlapping jets are a 
significant problem for the supersymmetric top decay. This results in a high mass tail on the top mass 
distributions reconstructed from the t bbbW^ decay. 



Jet assignment has also been studied for this decay channel and it has been found to be good 
for the mass distributions reconstructed using the Standard Model top decay channels. However, in the 
case of the supersymmetric top decay, the Higgs mass distributions show that there are contributions 
from both the misassignment of jets associated to other decay quarks and from jets associated to other 
hard partons in the event. However, only the latter contribution is observed in the corresponding top 
mass distributions, indicating that it is jets associated to the supersymmetric top decay that are being 
misassigned, not those from the SM top decay. The misassignment of jets from other hard partons also 
results in a high mass tail on the top mass distributions. The use of perfect b-tagging (based on jet-parton 
matching) suppresses this effect. This is consistent with the assumption that the other hard paitons come 
from initial state gluon radiation or parton showering, as in this case the misassigned jets are much more 
likely to be gluon- or light quark-initiated jets. 

It may be possible to reduce the impact of overlapping jets on the mass distributions by using a 
smaller jet cone radius or by using another jet finder, such as the kt algorithm [155,275]. The impact of 
detector-level jet finding and lepton identification must also be assessed. 

16. SEARCH FOR THE t bH+, H+ HiW, Hi bb CHANNEL IN CPX MSSM SCE- 
NARIO IN CMS EZl 

16.1 Introduction 

CP violation (CPX) in the Higgs sector of the Minimal Supersymmetric Standard Model (MSSM), when 
the Higgsino mass parameter fi, the gaugino mass parameters Mj and the trilinear couplings Aj are 
complex, allows the existence of the light neutral Higgs boson (niHi ^ ^0 GeV) and relatively light 
charged Higgs boson {m^+ < nit) in low tan/3 region not excluded by the LEP data because of the 
reduction of HiZZ coupling [270]. In CPX scenario the usual search channels may not be useful, 
because of the simultaneous reduction in the couplings of the Higgs boson to the vector boson pair and 
to the top quark pair, as it affects the Higgs boson production and decays rates. The one of the promising 
search channels in the CPX scenario proposed in [271] is the ti production when one of the top quarks 
decays as t ^ bH^ , — > HiW, Hi bb. It is due to the suppression of the HiZZ coupling leads to 
the enhancement of the H~^W~ Hi coupling in order to satisfy the coupling sum-rule. We investigated 
a feasibility for the discovery of the Higgs bosons in this channel using the full CMS detector simulated 
data. The results shown are preliminary. 

16.2 Event generation 

The signal events were generated using PYTHIA [183] with m(=175 GeV, m,/f^=51 GeV and mu+=\33 
GeV, corresponding to tan/5=5 and CP mixing angle ^>(CP)=90° in the CPX MSSM. The following 
decays were forced in PYTHIA: ti bW , t2 bH+, H+ WHi, Hi bb and both W bosons 
from the top decays were allowed to decay into all possible modes. The decay branching fractions were 
calculated using CPsuperH program [272]. The total cross section was calculated taking the next-to- 
leading order cross section for an inclusive ti production 840 pb [273] and multiplying by the branching 
ratios, Br(t bH+)=0.0\, Br(iJ+ HiW)=0.567, Br{t bW)=0.99, Br(ifi bb)=0.92 which 
gives the cross section 8.68 pb. 

The major background processes for this channel are tt + jets and tibb. The + jets background 
was generated using ALPGEN [180] with the MLM prescription for jet-parton matching [181, 182] at 
the PYTHIA shower simulation. The ti +2 jets (exclusive), ti + 3 jets and ti + A jets(inclusive) with 
jetpr >20 GeV were generated. The cross sections for these processes are shown in Table [TTl The tibb 
background was not considered yet in this study. 
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16.3 Simulation and Reconstruction 

The CMS detector was simulated using full GEANT4 [276] simulation and the reconstruction was done 
using the CMS simulation and reconstruction software CMSSW. No pileup events were included. We 
summarize briefly the object reconstruction methods [277] used in this analysis. Muons are reconstructed 
from the muon chambers and the silicon tracker and electrons are reconstructed from the tracks in the 
silicon tracker and the clusters in the electromagnetic calorimeter. The loose electron identification 
criteria were applied. The lepton isolation was done using the tracker isolation such that leptons are 
selected if sum pT of the tracks in a cone around the lepton (inner radius 0.015 and outer radius 0.25) is 
less than 3 GeV. The jets were reconstructed from the calorimeter towers using an iterative cone algorithm 
with the cone size 0.5. The jet energy was corrected using the Monte Carlo jet energy corrections. The 
missing Et was reconstructed from the calorimeter towers and corrected for the jet energy scale. The 
missing Et was also corrected for the muons by adding the muon momenta to the calorimeter missing 
Et. 



16.4 Event selection 

16.41 Primary selections 

The final state considered in this analysis consists of two light quarks, four b quarks, one lepton and 
neutrino: £u + qq' + bbbb. Since the neutral Higgs boson Hi is very light (51 GeV), the b quarks from 
the Hi bb decay are very soft as seen in Figl34](a,b). Only ~ 36% of events have both b quarks from 
the Hi decay with >20 GeV. The final state quarks in the event fall very close to each other in (r/, (j)) 
space. Fig. [34] (c) shows the sepai^ation in (i], (p) space between two closest quarks. Because of these 
reasons it is difficult to reconstruct six jets in the event corresponding to the six final states quarks. The 
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Fig. 34: (a) pt distribution of b quarks from top quark and from Hi decays, (b) pr distribution of b quarks from Hi decay, (c) 
the AR separation in (77, (j)) space between two closest final state quarks in the event. 



events with one isolated lepton with pT >20 GeV and six or more jets with Et >20 GeV were selected. 
The number of leptons in the event (electrons with pT > 10 GeV and muons with pt >5 GeV) passing 
the identification and the isolation were counted and the events with more than one lepton were rejected. 
The jets were b tagged using the track counting b-tagging algorithm. The three dimensional impact 
parameter significance of the second highest significance track in the jet was used as the b-discriminator 
parameter. The four highest discriminator jets with the discriminator value greater than 2.95 were tagged 
as b jets. 



16.42 Top mass reconstruction 

One W boson in the event was reconstructed from the lepton and the missing Et- The z-component of 
the missing energy was calculated using the W mass constraint. This yields the real solutions in nearly 
66% events. The events with the imaginary solutions were rejected. There are two possible solutions for 
the z-component of the missing energy which gives two possible candidates for the leptonically decaying 
W boson. The W boson decaying hadronically was reconstructed from the jets not tagged as b jets. All 
jet pairs with the invariant mass within the mvi/± 20 GeV mass window were considered as possible 
candidates. The di-jets invariant mass for the jets matching to quarks from the W boson decay is shown 
in Figj35](a). The momenta of the two top quarks were reconstructed simultaneously from four b-tagged 
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Fig. 35: (a) the di-jet invariant mass of the jets matching to quarks from W decay, (b) the top-quark mass reconstructed from the 
bjj final state after the minimization of 5M. (c) the top-quark mass reconstructed from bbbjj final state after the minimization 

ofSM. 



jets, two W ^ iv candidates and candidates for the hadronically decaying W boson. The jets and the 
W boson candidates were assigned to the two top quarks by minimizing the 5M, where (5M is defined as 



5M = y' (nitopi - mtop)^ + (mtop2 - mtop)^ + (mw(hadronic) - mw)^, (31) 

there mtopi is reconstructed from one b-tagged jet and one W boson candidate, mtop2 is reconstructed 
from three b-tagged jets and one W boson candidate, mtop is the generated top-quark mass (175 GeV) 
and the mw is the W boson mass (PDG value). 

The top-quark mass distributions reconstructed from three jets {bjj) and five jets {bbbjj) after the 
minimization of 5M are shown in Figl35](b,c). One can see that the top-quark mass distribution from 
the bbbjj final state is very wide and has a big tail. It is because of wrong assignment of the jets or W 
candidates to the top while minimizing 5M. The events with the two top-quark reconstructed masses 
within the mtop± 30 GeV mass window were selected. Table [TT] shows the initial cross sections for the 
signal and background processes, the number of Monte-Carlo events remaining after each selection step 
and the cross sections after all selections. 

16.43 Reconstruction of the neutral Hi and charged Higgs bosons 

Since it is not known what pair of the b-tagged jets from the reconstructed top quark decay chain 
t — > bbbW is coming from the Hi bb decay, all three b-tagged jet pairs were considered as the 

^'Xhe W ^ £u and — + jj reconstruction step selects events with the positive solution for z-component of and 
with at least one jet pair having the di-jet mass within the m^yi 20 GeV mass window; the top-quark mass reconstruction step 
requires that the two top-quark reconstructed masses are within the mtopi 30 GeV mass window. 



Table 17: The initial cross sections for the signal and background processes, the number of Monte Carlo events remaining after 
each selection steps and the cross section after all selections. 





signal 


tt + 2 jets 


tt + 3 jets 


tt + 4 jets 






(exclusive) 


(exclusive) 


(inclusive) 


cross section, pb 


O ZTO 


100 


40 


61 


number of MC events analyzed 


193884 


241000 


71000 


94000 


(corresponding luminosity, fb^^) 


(22.35) 


(2.41) 


(1.775) 


(1.54) 


isolated lepton pr > 20 GeV 


41035 


57920 


16915 


22214 


> 6 jets Et > 20 GeV 


21389 


36315 


14479 


21866 


4 b-tagged jets with discr. >2.95 


881 


371 


248 


1069 


W ^ iu and ^ jj reco 


379 


158 


132 


602 


top-quark mass reconstruction 


83 


4 


1 


7 


cross section after all selections, fb 


3.71 


1.66 


0.56 


4.54 



possible candidates. The invariant mass of b-tagged jet pairs, mi,f) is shown in Fig. [36] (left plot) for the 
background and the signal plus background. The right plot in Fig. |36] shows, fitted by the Gaussian the 
nibb distribution of the signal plus background. The mean value of the fitted distribution is close to the 
generated mass of Hi. The charged Higgs boson was reconstructed from the two b-tagged jets and W 
boson, where the b-tagged jet pair was chosen as the jet pair with the invariant mass closest to the peak of 
the rribb mass distribution and within the window ±20 GeV around the fitted mean value. The invariant 
mass distribution of the charged Higgs boson reconstructed in this way, nibbw is shown in Figl37l 
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Fig. 36: The invariant mass of the b-tagged jet pairs from the reconstructed top quark decay chain t bbbW 



The available Monte Carlo statistics of tt + jets background events for this study was only order 
of ~ 2 fb^^, thus it can not be simply rescaled in order to produce the smooth shape of nibb and rubbw 
distributions expected for 30 fb^^ after all selections. We have obtained, however that the shape of rubb 
and nibbw distributions is almost the same after relaxing the cut on the b-discriminator value. Fig. [38] 
shows the nibb (left plot) and nibbw (right plot) distributions for four different b-discriminator cuts: 0, 
1.0, 1.5 and 2.0. 
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Fig. 37: The invariant mass of two b-tagged jets and W boson, where two b-tagged jets were chosen with the mass closest to 
the peak of the rribb mass distribution and within the window ±20 GeV around the fitted mean value 
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Fig. 38: The mtb (left) and rubbw (right) distributions after all selections for four different cuts on the b-discriminator value: 
0, 1.0, 1.5 and 2.0. 



16.5 Results 

The simple selection strategy described in the previous sections yields S=110 signal events and B=203 
ti + jets background events expected with 30 fb^^. The ttbb background still need to be taken into 
account. The uncertainty due to the Monte Carlo statistics on the ti + jets backgrounds is ~ 30%. The 
experimental systematic uncertainty was estimated by taking into account the systematic uncertainties on 
the lepton identification (2%), the b-jet tagging (5% per jet), the jet energy scale (5% per jet), the missing 
transeverse energy scale (10% on the raw calorimeter energy scale and 5% on the jet energy scale) and 
the luminosity uncertainty (5%). It leads to the total systematic uncertainty 22.5% (the uncertainty due 
to the jet and the missing Et scale only is 8.8%). The significance is calculated as S/^/B + AB'^, 
where AB is the experimental systematic uncertainty on the background. In order to get the pessimistic 
value for the significance, the Monte Carlo statistical uncertainty was added to the total background: 
B=203-i-60=263 events. The signal significance is then 110/V263 + 59^=1.8. The uncertainty on the 
theoretical leading-order cross section of the tt + n jets, (n>2) processes is >50%. 

One can see that the discovery potential is restricted by both the experimental and the theoretical 



uncertainties. The uncertainties can be partially reduced if the number of the background events and 
the 77166 and rrifefeVF mass shapes can be extracted from the data. The shapes can be evaluated from the 
data with the ratio S/B<<1 when the relaxed cut on the b-discriminator value is used (see Fig.[38l). The 
background normalization on the number of events with the relaxed b-jet tagging will eliminate the jet 
and the missing Et scale uncertainties, the luminosity uncertainty and partially reduce the b-tagging 
uncertainty which dominates the experimental uncertainty. It will also reduce the absolute background 
prediction uncertainty from the theory, since only the ratio of tt + jets and ttbb cross sections need to be 
used. The further, more detailed investigations of this channel is foreseen in CMS. 
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Part V 

NMSSM HIGGS BOSONS 

17. LES HOUCHES BENCHMARK SCENARIOS FOR THE NMSSM 
17.1 Introduction 

The next-to-minimal supersymmetric extension of the Standard Model (NMSSM) [278,279], in which 
the spectrum of the minimal supersymmetric extension (MSSM) is extended by one singlet superfield, 
is interesting in many respects. Compared to the MSSM, it solves in an elegant way the so-called /i 
problem, has less fine tuning and can induce a rather different phenomenology in the Higgs and neutralino 
sectors. Given the possibility of a quite different phenomenology, it is important to address the question 
whether such NMSSM specific scenarios will be probed at the LHC. In particular, it would be crucial 
to make sure that at least one Higgs particle should be observed at the LHC for the planned integrated 
luminosity or try to define regions of the NMSSM parameter space in which more Higgs states than 
those available within the MSSM are visible. However, a potential drawback of the NMSSM, at least 
in its non-constrained versions, is that it leads to a larger number of input parameters to deal with. In 
particular, it is clearly unfeasible to make multi-dimensional scans over the free inputs of the NMSSM 
when performing complete/realistic simulations to address the two points mentioned above. 

An alternative approach is to resort to a few benchmark scenarios which embodying the most 
peculiar/representative phenomenological features of the model's parameter space, which can be subject 
to full experimental investigation, without loss of substantial theoretical information. Building on the 
experience of Ref. [280], we define in this note benchmark points which fulfill the present collider and 
cosmological constraints using the most-up to date tools to calculate the particle spectra. We work in the 
framework of a semi-constrained NMSSM (cNMSSM) where the soft Supersymmetry (SUSY) breaking 
parameters are defined at some high scale, typically that of grand unification theories (GUTs). This 
approach leads to a much more plausible sparticle spectrum, allows to relate features of the Higgs sector 
to properties of the neutralino sector and, at the same time, still contains the distinctive phenomenological 
features of the NMSSM that are suitable for intensive phenomenological/experimental investigation. The 
emphasis is primarily on the different possible scenarios within the Higgs sector and the implication for 
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Higgs searches at the LHC. In particular, we propose five benchmark points which lead to Higgs-to-Higgs 
decays or a light Higgs spectrum but with reduced Higgs-gauge boson couplings, which are known to 
be rather difficult to probe at the LHC. 

17.2 The Model and Its Spectrum 

We confine ourselves to the NMSSM with a scale invariant superpotential given, in terms of (hatted) 
superfields with only the third generation (s)fermions included, by 

W = XSHuHa + f 53 + htQEj^R - hQHab'k - KLHaT^. (32) 

The first two terms substitute the jiHuHd term in the MSSM superpotential, while the three last terms are 
the usual generaUzation of the Yukawa interactions. The soft SUSY breaking terms consist of the scalar 
mass terms for the Higgs, sfermion and gaugino fields and the trilinear interactions between the sfermion 
and Higgs fields. In an unconstrained NMSSM with non-universal soft terms at the GUT scale, the three 
SUSY breaking masses squared for H^, and S are determined through the minimization conditions 
of the scalar potential. Thus, the Higgs sector of the NMSSM is described by the six parameters 

A , AC , Aa , ^K, tan/3 = {Hu)/{Hd) and //eff = X{S) . (33) 

As the number of input parameters is rather large, one can attempt to define a constrained (cNMSSM) 
model, similar to the minimal supergravity model or cMSSM, in which the soft SUSY breaking pa- 
rameters are fixed at the GUT scale, leading to only a handful of inputs. One can thus impose uni- 
fication of the gaugino, sfermion and Higgs mass parameters and the triUnear couphngs at Mqut'- 
Mi^2,3 = Mi/2,mp^ = niHi = rno,Ai = Aq. The fully constrained cNMSSM has two additional 
parameters, A and k, beyond the above and the correct Mz value imposes one constraint. Hence, a 
priori, the number of inputs in the cMSSM and the fully constrained cNMSSM is exactly the same. 

In practice, it is convenient to use the analytic form of the three minimization conditions of the 
NMSSM effective potential and, for given Mz, tan /3, A and all soft terms at the weak scale except 
for m|, these can be solved for |/iefr| (or | (S) |), k and m|; sign(^efr) can still be chosen at will. 
Here, we will relax the hypothesis of complete unification of the soft terms in the singlet sector, m| / 
rriQ and A,^ / Aq at Mgut, since the singlet could play a special role. In addition, for some of the 
benchmark points, we will also relax the unification hypothesis for rnjj^ and m|^^ and for one scenario, 
the hypothesis Ax = Aq. Such points in parameter space can have additional unconventional properties, 
whose phenomenology should also be investigated. 

Following the procedure employed by the routine NMSPEC within NMSSMTools [281], which 
calculates the spectra of the Higgs and SUSY particles in the NMSSM, a point in the parameter space of 
the cNMSSM is defined by the soft SUSY breaking terms at Mqut (except for the parameter m|), tan /? 
at the weak scale, A at the SUSY scale (defined as an average of the first generation squark masses) and 
the sign of the parameter fics- The parameters k, m| and |/ieff| are determined at the SUSY scale in terms 
of the other parameters through the minimization equations of the scalar potential. The renormalisation 
group equations (RGEs) for the gauge and Yukawa couplings and those for the soft terms are integrated 
between Mz and Mqut defined by gauge couplings unification. For the most relevant Standard Model 
parameters, we chose as{Mz) = 0.1172, rnh{mi,)^'^^ = 4.214 GeV and m^°p = 171.4 GeV. 

After RGE running is completed, the Higgs, gluino, chargino, neutralino and sfermion masses are 
computed including dominant one-loop corrections to their pole masses. All the Higgs decay branch- 
ing ratios (BRs) into SM and SUSY particles are determined including dominant radiative corrections. 

Subsequently, the following Tevatron and LEP constraints are applied: i) Direct searches for the LSP 
neutralino and invisible Z decay width, ii) direct bounds on the masses of the charged particles h^, x^, 
q, I and the gluino; in) constraints on the Higgs production rates from all channels studied at LEP. 



Light hi {i = 1,2) scalar states (with M/j- ^114 GeV) can still be allowed by LEP constraints, if 
the Z-Z-hi coupling is heavily suppressed or the lightest pseudoscalar a\ state has Ma^ < 10 GeV such 
that hi decays dominantly into a\a\ states but the U) decay of the a\ is impossible. Constraints from the 
decays hi — > oioi — > 4r allow for M/^^ down to ~ 86 GeV. Note that LEP constraints are implemented 
only for individual processes and that combinations of different processes could potentially rule out 
seemingly viable scenarios. Finally, experimental constraints from B physics are taken into account, and 
we require that the relic abundance of the NMSSM dark matter (DM) candidate, the lightest neutralino 
X]* which can be singlino-like, matches the WMAP constraint 0.094 < J^cdm^^ ^ 0.136 at the la level. 

17.3 The Benchmark Points 

In the Higgs sector of the NMSSM, two different types of difficult scenarios have been pointed out, 
depending on whether Higgs-to-Higgs decays are kinematically allowed or forbidden; see e.g. Ref. [280]. 

Within the first category, there are two possibilities, each associated with light scalar/pseudoscalar 
Higgs states: (i) The hghtest CP-odd a\ state is rather light, Ma^ < 40-50 GeV, and the hghtest CP- 
even hi particle has enough phase space for the decay into two ai particles, hi —>■ oioi, to be allowed 
and dominant. The ai state will mainly decay into r+r^ if M^^ < 10 GeV or to t^t~ (~ 10%) 
and bb (90%) states if Mq^ > 10 GeV. One would have then the possibilities hi — > aioi — > 4r and 
hi — > oiOi — > 4t, 46 and 2t26 for the hi state which can have a mass that is either close to its theoretical 
upper limit of 130 GeV or to the lower limit of 90 GeV. (ii) The lightest CP-even hi boson is relatively 
light, Mh-^ ^ 50 GeV, and decays into bb pairs (the situation where Mh^ ^ 10 GeV is very constrained 
by LEP data). In this case, the next-to-lightest CP-even /i2 state is SM-like with a mass below ~ 140 
GeV and can decay into two hi bosons leading to the final topologies /12 ^ /ii /ii — > 4r, 2r25 and 46. 

The second category of scenarios, where Higgs-to-Higgs decays are suppressed, includes regions 
of the parameter space where the five neutral Higgs particles are relatively light, with masses in the range 
90-180 GeV, which opens the possibility of producing all of the them at the LHC, but with couplings to 
gauge bosons that are reduced compared to the SM Higgs case. This scenario is similar to the so-called 
"intense coupling regime" of the MSSM [282] but with two more neutral Higgs particles. 

We propose five benchmark points of the NMSSM parameter space, PI to P5, in which the above 
mentioned scenarios are realized (see Ref. [283] for more details). Each point is representative of dis- 
tinctive NMSSM features. Points PI to P3 exemplify scenarios where hi decays into light pseudoscalar 
states decaying, in turn, into bb or r+r^ final states; these points can be realized within the cNMSSM 
with nearly universal soft terms at Mqut. the exception being the parameters m| and A^. P4 illustrates 
the NMSSM possibility of a very light hi and can be obtained once one relaxes the universality con- 
ditions on the soft SUSY breaking Higgs mass terms, Mh^ / 7^ mo. Point P5 corresponds to 
the case where all Higgs bosons are rather light and can be obtained if one allows additionally for the 
inequality Ax Aq. In all cases, the input parameters as well as the resulting Higgs masses and some 
decay information are given in Table [TSl the main characteristics of the Xi DM candidate are also given. 
Next, we summarize the most relevant phenomenological properties of the benchmark points. 

In the first two points PI and P2, the lightest hi CP-even state has a mass of Mhj ~ 120 GeV 
and is SM-like with couplings (relative to that of the SM Higgs) to gauge bosons Ri, top quarks ti 
and bottom quarks bi, which are almost equal to unity. The lightest CP-odd ai boson has a mass of, 
respectively, 40.5 GeV and 9.09 GeV. In both cases PI and P2, the decay channel hi —>■ aiOi is largely 
dominating with a BR very close to 90%, while the decays hi — > bb and r+r^ are suppressed by an order 
of magnitude when compared to the SM case. The most relevant difference between the two scenarios 
concerns the mass and decays of the lightest pseudoscalar state. In PI the ai boson decays into b quarks 
and T leptons with rates of ~ 90% and ~ 10%, respectively. In contrast, in P2 the pseudoscalar ai state 
with its mass M^^ ~9.09 GeV decays dominantly into r+r^ pairs, with a rate above 80%. 

For point P3, the same inputs of points PI and P2 are chosen except for the and A parameters. 



Table 18: Input and output parameters for the five benchmark NMSSM points. 



Point 


PI 


P2 


P3 


P4 


P5 


GUT/input parameters 




sign(/ieff) 


+ 


-1- 


-1- 




-1- 


tan/3 


10 


10 


10 


2.6 


6 


mo (GeV) 


174 


174 


174 


775 


1500 


Ml/2 (GeV) 


500 


500 


500 


760 


175 


Ao 


-1500 


-1500 


-1500 


-2300 


-2468 


Ax 


-1500 


-1500 


-1500 


-2300 


-800 


A^ 


-33.9 


-33.4 


-628.56 


-1170 


60 


NUHM: Mhj (GeV) 








880 


-311 


NUHM: Mh„ (GeV) 








2195 


1910 


Parameters at the SUSY scale 




A (input parameter) 


0.1 


0.1 


0.4 


0.53 


0.016 


K 


0.11 


0.11 


0.31 


0.12 


-0.0029 


Ax (GeV) 


-982 


-982 


-629 


-510 


45.8 


A^ (GeV) 


-1.63 


-1.14 


-11.4 


220 


60.2 


M2 (GeV) 


392 


392 


393 


603 


140 


p-eff v^evj 


968 


968 


936 


-193 


303 


CP-even Higgs bosons 




lUh^ (OeV) 


120.2 


120.2 


89.9 


32.3 


90.7 




0.072 


0.056 


7 X 10"* 


0.918 


0.895 


BR(/)i ^ r ' r ) 


0.008 


0.006 


7 X 10 " 


0.073 


0.088 


BR(/ii —> aioi) 


0.897 


0.921 


0.999 


0.0 


0.0 


nih^ (GeV) 


998 


998 


964 


123 


118 


ma, (GeV) 

"■3 \ ' / 


2142 


2142 


1434 


547 


174 


CP-odd Higgs bosons 




m„ (GeV) 


40.5 


9.09 


9.13 


185 


99.6 


BR(ai ^ 66) 


0.91 


0. 


0. 


0.62 


0.91 


BR(ai -» T+r") 


0.085 


0.88 


0.88 


0.070 


0.090 


(GeV) 


1003 


1003 


996 


546 


170 


Charged Higgs boson 




mh± (GeV) 


1005 


1005 


987 


541 


188 


LSP 




rn.o (GeV) 


208 


208 


208 


101 


70.4 


2 


0.099 


0.099 


0.130 


0.099 


0.105 



which are now varied as to have a lighter hi state. This again leads to a pseudoscalar ai boson which has 
approximately the same mass as in scenario P2, Ma-^ ~ 9.96 GeV, and which decays almost exclusively 
into r+r~ final states. The difference between P3 and P2 is the lightest CP-even Higgs boson hi, which 
has a mass M^^ ~ 90 GeV, lower than in scenarios PI and P2. In this case, and although hi is still 
SM-like, i.e. exhibiting couplings to gauge bosons, top and bottom quarks that are very close to those of 
the SM Higgs boson, it decays nevertheless almost exclusively into ai pairs, with a rate close to 100%. 
Another difference between P2 and P3 is that in the former case, the interesting decay mode hi aiZ 
is kinematically possible but the rate is rather small, BR(/ti aiZ) ~ 3%. 

Note that in all these first three points, the heaviest neutral Higgs particles /i2, ^3 and 02, as well 
as the charged Higgs states h^, all have masses close to, or above, 1 TeV The main decay modes are into 
hh and for the neutral and th for the charged states, as tan (5 is not too large and the ti-Higgs couplings 
are not very strongly suppressed, while the BRs for the neutral Higgs-to-Higgs decays, in particular the 
channels /t2 hihi and /12 a\ai, are very tiny, not exceeding the permille level. Regarding the 
properties of the DM candidate, PI, P2 and P3 exhibit a lightest neutralino which is bino-like, with mass 
is m o ~ 208 GeV. In all three cases, the correct cosmological density, ricDM^i^ — 0.1, is achieved 



through the co-annihilation with the fi slepton, which has a mass comparable to that of the LSP. 

Point P4 corresponds to a scenario in which the CP-even boson hi is very light, M^^ = 32.3 GeV 
and singlet-like and predominantly decays into bb pairs, with BR(/ii bb) = 92%, and to a smaller 
extent into r pairs with BR(/ii t^t~) 7%. The CP-even /i2 boson has a mass of Mh,-, ^ 123 GeV 
and is SM-like, with normalized couplings to W/Z and t/b states close to unity. However, it mostly 
decays into two hi bosons, BR(/i2 hihi) ~ 88% and the dominant SM-like bb decay mode occurs 
only at a rate less than 10%. The lightest CP-odd particle is not very heavy, Ma^ = 185 GeV, and decays 
mostly into fermion pairs, with BR(ai bb) ~ 61% and BR(ai t~^t~) ~ 7%; the other dominant 
decay is the interesting channel ai hiZ which has a rate of the order of 30%. Finally, the heaviest 
CP-even h^,, CP-odd 02 and the charged h^ particles have masses in the 500 GeV range and will mostly 
decay, as tan (3 is small, into tt/tb final states for the neutral/charged states. All these features make the 
phenomenology of point P4 rather different from that of points PI to P3 discussed above. To achieve 
a correct cosmological relic density, the common sfermion and gaugino mass parameters at Mqut are 
close to 1 TeV At the SUSY scale, one thus finds a higgsino-singlino-like neutralino LSP, whose mass 
is m^o ~ 100 GeV and LSP annihilation essentially leads to WW and Zhi final states. 

Al 

Finally, point P5 is characterized by having all Higgs particles relatively light with masses in the 
range 90 to 190 GeV. Here, both A and k are relatively small. The three CP-even Higgs bosons with 
masses of 91, 118 and 174 GeV, respectively, share the couphngs of the SM Higgs boson to the gauge 
bosons with the dominant component being taken by the /12 state. The pseudoscalar Higgs bosons have 
masses M^j ~ 100 GeV and Ma^ ~ 170 GeV, while the charged Higgs particle is the heaviest one with 
M^± ~ 188 GeV. Here, all the neutral Higgs-to-Higgs decays are kinematically disfavored; this is also 
the case of neutral Higgs decays into into lighter Higgs states with opposite parity and gauge bosons. The 
only non-fermionic two-body Higgs decays are thus /t^ Whi and /13 WW, but as the involved 
Higgs-gauge boson couplings are small, the BRs are tiny. Here, the LSP with a mass m^o ~ 70 GeV, is 
a bino-like neutralino but it has a small non-negligible higgsino component. The value i^cDM^^ — 0.1 
is achieved through the annihilation processes XiXi ~^ bb,T~^T~, with s-channel exchange of Higgs 
bosons. 

17.4 Expectations at the LHC 

In the cases discussed here, at least one CP-even Higgs particle hi has strong enough couplings to 

massive gauge bosons and top quarks, Ri,ti ~ 1, to allow for the production at the LHC in one of the 
main channels which are advocated for the search of the SM Higgs particle [284]: i) gluon-gluon fusion, 
gg — > hi, ii) vector boson fusion (VBF), qq qqW*W*,qqZ*Z* qqhi with two forward jets and a 
centrally decaying Higgs boson, in) Higgs-strahlung (HS), qq' Whi and qq Zhi, with the gauge 
boson decaying leptonically, iv) associated production with heavy top quark pairs qq/gg tthi. 

In scenarios PI to P3, this CP-even hi particle is the hi boson which has i?i ~ ti ~ 61 ~ 1, 
but which decays most of the time into a pair of light pseudoscalar Higgs particles, hi — > aioi, which 
subsequently decay into light fermion pairs, ai —>■ bb and T"^r~. In scenario P4, this particle is the h2 
boson which decays most of the time into a pair of hi particles, /12 h\h\, which again decay into light 
fermion pairs. In these four cases, the backgrounds in both gg ^ hi ^ 4: f and qq/gg tihi tt + Af, 
with f = b,T, processes will be extremely large and only the VBF (owing to the forward jet tagging) 
and eventually HS (due to the leptons coming from the decays of the gauge bosons) can be viable at 
the LHC. In P5, the particle that has couplings to gauge bosons and top quarks close to those of the SM 
Higgs boson is the /12 boson which decays into bb and r+r" final states with BRs close to 90% and 10%, 
respectively. Here again, the gg fusion and presumably associated production with top quarks cannot be 
used since the interesting decays such as /t2 WW*, ZZ* and 77 are suppressed compared to the SM 
case. Thus, in this case, only the channels qq qqh2 qqr'^r" and eventually qq' Wh2 IvlA) 
seem feasible. The state hi has still non-negligible couplings to gauge bosons and top quarks which 
lead to cross sections that are "only" one order of magnitude smaller than in the SM. Since here again. 



only the decays hi — > bb (90%) and r+r^ (10%) are relevant, the only channels which can be used are 
the VBF and HS processes discussed above, but one needs a luminosity 10 times larger to have the same 
event samples as in the SM. 

Several theoretical studies have been performed in the past to assess the potential of the LHC to 
observe NMSSM Higgs particles in some scenarios close to those presented here; see Ref. [283] for an 
account. Recently, the ATLAS and CMS collaborations started investigating some channels, the main 
focus being on the VBF production process pp qqhi and to a lesser extent HS via qq Whi £vhi, 
with the SM-like hi state decaying into hi — > aioi — > 4r, corresponding to scenario P2 and P3. The 
ATLAS collaboration is analyzing the Afi + Aur + channel from VBF, requiring three leptons to 
be observed and, for triggering, one or two high-pT leptons (pT > 20 or 10 GeV) [285]. CMS is 
currently considering the fi^fj-^T^^T^^ final state containing two same sign muons and two r jets [286]. 
Despite of the missing energy and the possibility of missing one lepton, the mass of the /i2 state could 
be reconstructed with the help of the coUinear approximation. The performance of the algorithms to 
observe the signals and the effects of the various backgrounds are under study. 

17.5 Conclusions 

The NMSSM is a very interesting supersymmetric extension of the SM as it solves the notorious 
problem of the MSSM and it has less fine tuning. It also leads to an interesting collider phenomenology 
in some cases, in particular in the Higgs sector, which is extended to contain an additional CP-even and 
a CP-odd state. Compared to the SM and MSSM, the searches for the NMSSM Higgs bosons will be 
rather challenging at the LHC in scenarios in which some neutral Higgs particles are very light, opening 
the possibility of dominant Higgs-to-Higgs decays, or when all Higgs bosons have reduced couplings 
to the electroweak gauge bosons and to the top quarks. These scenarios, for which we have provided 
benchmark points in a semi-unified NMSSM which involves a rather limited number of input parameters 
at the grand unification scale and which fulfills all present collider and cosmological constraints, require 
much more detailed phenomenological studies and experimental simulations to make sure that at least 
one Higgs particle of the NMSSM will be observed at the LHC. 
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18. PARAMETER SCANS IN TWO INTERESTING NMSSM SCENARIOS 
18.1 Introduction 

In the past, proposals for interesting points in the pai^ameter space of the Next-to-Minimal Supersymme- 
tric Standard Model (NMSSM) [287-292] have been made (see e.g. Refs. [293-295]). A new study 
proposes benchmark points for the constrained NMSSM [296]. To evaluate the discovery potential of 
NMSSM particles at collider experiments like the Large Hadron Collider (LHcjH it is furthermore de- 
sirable to define two-dimensional benchmark scans which include regions of typical and experimentally 
challenging NMSSM phenomenology. In the following, two such parameter scans over the Higgs sector 
of the NMSSM are proposed for this model. Both scans include a benchmark point from Ref. [295]. 

"Contributed by: I. Rottlander and M. Schumacher 

"""A proton-proton coUider with a design center-of-mass energy of 14 TeV. First physics runs are expected for 2008. 




Fig. 39: Hi mass [GeV] in the Reduced Couplings Scenario 




Fig. 41: ^1 mass [GeV] in the Reduced Couplings Scenario 




Fig. 40: H2 mass [GeV] in the Reduced Couplings Scenario 




Fig. 42: mass [GeV] in the Reduced Couplings Scenario 



18.2 The NMSSM 

In the Minimal Supersymmetric Standard Model (MSSM), the value of the Higgs-Higgsino mass param- 
eter /X is not confined by theory, but it is experimentally constrained to lie at the weak scale or else large 
fine-tuning is required (the so called ^-problem). In the NMSSM, an additional neutral singlet superfield 
S is added to the MSSM. After symmetry breaking, // is then given by the product of the vacuum expec- 
tation value of the bosonic component of S {<s>) and a new coupling constant A. Constraints from the 
Higgs potential minimization strongly prefer <s> to lie at the weak scale. The right value of /^i is thus 
obtained naturally. 

The resulting model contains the whole particle spectrum of the MSSM with an additional neutral scalar 
boson, a pseudoscalar boson and a neutral fermion ("singlino"). The two additional neutral scalar bosons 
contained in S mix with the MSSM Higgs bosons to form the five neutral Higgs bosons of the NMSSM: 
three CP-even bosons Hi, H2, and two CP-odd Higgs bosons Ai, A2. The neutral fermion mixes 
with the four neutralinos of the MSSM, thus, the model contains in total five neutral fermion states. Since 
no charged particles are added, the features of the other MSSM particles, including the charged Higgs 
boson H^, are only modified marginally. The maximally allowed mass of the lightest NMSSM scalar 
Hi is about 10 GeV higher than the bound for h in the MSSM [297]. 

In the NMSSM, the Higgs sector can at tree level be described by six parameters. Usually, these are 
chosen to be the coupling parameters of S (A, k, A\, A^), fj, and the ratio of the vacuum expectation 
values of the Higgs fields, tan (3. In the here defined two-dimensional parameter scans, A and k are 
varied. Variation of the other parameters also changes the features of the Higgs sector, however, a X-k 
variation was found to be sufficient to cover the most important phenomenology types in the two scans 
described here. 

To calculate the NMSSM particle spectra and exclusion constraints from theory and LEll^. the program 
NMHDECAY [298-300] was used. The mass parameters were chosen as Mi = 500 GeV, M2 = 1 TeV, 



'The Large Electron Positron Collider, which ran until 2000 at center-of mass energies up to 209 GeV. 



Fig. 43: H2^H\H\ branching ratio in tlie Reduced Couplings 
Scenario 




Fig. 45: H2 vector boson coupling relative to its SM-value in 
the Reduced Couplings Scenario 



Fig. 44: H\ vector boson coupling relative to its SM-value in 
the Reduced Couplings Scenario 




Fig. 46: Hz vector boson coupling relative to its SM-value in 
the Reduced Couplings Scenario 



M3 = 3 TeV and Mgusy = 1 TeV. The trilinear soft supersymmetry-breaking parameters were set to 
At = Ab = Ar = 1.5 TeV, the top quark mass to 172 GeV. 

18.3 The Reduced Couplings Scenario 

Due to the mixing with the gauge singlet states, the NMSSM Higgs bosons can have reduced gauge 
coupUngs and thus reduced production cross sections compared to the Standard Model (SM) or the 
MSSM case. A light scalar with reduced gauge couplings and a mass below 1 14 GeV is still unexcluded 
by LEP. 

The here proposed scenario is a \-k scan with parameters given in Table [T9l The point with A = 0.0163 
and K = —0.0034 is described as having the lowest statistical significance found in a region without 
Higgs-to-Higgs decays in Ref. [295]. 

The masses of all six Higgs bosons in this scenario are smaller than about 300 GeV The Hi is very 
light, down to values of about 20 GeV in an unexcluded region with small negative k (Figi39l). Since 
the H2 has a SM-like mass around 120 GeV in the entire plane (Fig|40l). there is a region where the 
decay H2^HiHi is allowed with a small branching ratio of at maximum 6% in the unexcluded region 
(Figj43]). The Ai mass ranges from about 55-100 GeV (FiglTTI) in the allowed parameter region, whereas 
the H^, A2, and are approximately degenerate in the entire plane, but with small differences in mass 
for large negative k. The mass of the ranges from about 150 to 300 GeV, the mass of the A2 from 
about 140 to 300 GeV and the charged Higgs boson mass from about 165 to 300 GeV in the unexcluded 
region (Figj42l). 

In Figures |44l|45]&|46l the vector boson couplings of the scalar bosons are given as an example gauge 
coupling. Higgs boson couplings to gluons and up-type fermions vary similarly. The Hi and H3 gauge 



Fig. 47: Hi mass [GeV] in the Light Ai Scenario 



Fig. 48: H2 mass [GeV] in the Light Ai Scenario 



coupling^ are highly suppressed in most of the parameter plane, reaching sizeable values only in the 
LEP excluded region at large negative k. The H2 has SM-like gauge couplings in large parts of the 
parameter plane. In the unexcluded region close to the benchmark point from Ref. [295], the vector 
boson couplings are reduced down to about 80% of their SM-value. Gauge couplings of the Ai and A2 
are highly suppressed for all considered parameter values. 

To summarize, this scenario is characterized by a region with a very light Hi close to the upper exclusion 
bound, where H2^HiHi decays are possible, a region with a SM-like H2 in the middle of the allowed 
parameter space, and a region with reduced couplings of the H2 at large negative k close to the lower 
exclusion bound. 



Table 19: Higgs sector parameters of the proposed scenarios 



Scenario 


A-range 


K-range 


Ax [GeV] 


A^ [GeV] 


^ [GeV] 


tan (3 


Reduced couplings 


- 0.025 


-0.005 - 


-70 


-54 


-284 


5.7 


Light Ai 


- 0.55 


-0.2 - 0.6 


-580 


-2.8 


-520 


5.0 




18.4 The Light Ai Scenario 

Unlike in the MSSM, the mass of the lightest pseudoscalar Ai is in the NMSSM not closely coupled to 
the masses of the scalar Higgs bosons and might thus lie well below the H1/H2 masses. In such a case, 

^*The term 'gauge couplings' here and in the following always excludes the Higgs boson coupling to down-type fermions 
which may be enhanced with respect to the SM-value, but are still too small to have an impact on the Higgs boson discovery 
potential with the here used tan /3 values around 5. 




the decay chain Hi/2^AiAi can be the dominant decay mode of the lightest scalars. 

The here described scenario is also a \-k scan with parameters given in Table [T9l The point with A = 

0.22 and n = -0.1 has been described in Ref. [295]. 

Here, the lightest scalar Hi has a mass around 120 GeV in the unexcluded region (FigflTl). The Ai is 
very light with masses up to about 60 GeV (Figj49l). so that the decay Hi^AiAi is possible in the entire 
parameter plane with exception of a small region at very small A and k (FigBTI). In the unexcluded region 
with large A and k, this decay reaches branching ratios above 90%. Areas with a smaller branching ratio 
exists for smaller A and k. The other Higgs bosons are rather heavy (Figs l48l50l) . with the H^, A2 and 

being approximately degenerate in large parts of the parameter plane. 
For Ai masses larger than 2mb, about 90% of the lightest pseudoscalar bosons decay to bottom quarks. 
In these regions, the decay chains Hi~^AiAi—^bbbb and Hi—>-AiAi—>-bbTT are important. In small 
regions at the borders of the unexcluded region, the Ai is so light that the decay chain Hi—>-AiAi—^tttt 
prevails (Figi52l). In the narrow unexcluded band around A « 0.25, the couplings of the Ai to fermions 
are heavily suppressed. Here, the decay chain Hi^AiAi — > 7777 is dominant. 

The gauge couplings of the Hi are SM-like in the entire allowed parameter region (Figj53l). The gauge 
couplings of the H2 are sizeable only in a small excluded region with k- values close to zero (Figl54l). 
All other Higgs bosons have highly suppressed gauge couplings in the entire parameter plane. 




Fig. 53: Hi vector boson coupling relative to its SM- value in Fig. 54: H2 vector boson coupling relative to its SM-value in 
the Light Ai Scenario the Light Ai Scenario 



18.5 Conclusions 

Two interesting two-dimensional NMSSM scans were described and proposed as possible benchmarks 
for NMSSM Higgs boson searches. These two scans cover the four main, for the NMSSM typical 
phenomenology types, for which a discovery of Higgs bosons at future experiments like the LHC might 
be difficult: 

• A region with very light scalar Hi. 



• A region with reduced gauge couplings of an otherwise SM-like scalar H2. 

• Regions with dominant Hi^AiAi^bbbb/bbTT decays of an otherwise SM-like scalar Hi. 

• Regions with dominant Hi^AiAi^tttt decay of an otherwise SM-like scalar Hi. 

Another example of an experimentally challenging phenomenology type not covered here is a dominant 
-ffi— >cc decay [294]. Also the region where the mass of the lightest scalar is maximal [297] could prove 
interesting for Higgs boson discovery . 

19. THE NMSSM NO-LOSE THEOREM AT THE LHC: THE SCOPE OF THE 4t CHANNEL 
IN HIGGS-STRAHLUNG AND VECTOR BOSON FUSION El 

19.1 Introduction 

As emphasised in Sect. 117.1 (see also [283]), the NMSSM has obvious advantages with respect to the 
MSSM. In constrast, it is not certain that at least one Higgs boson can be found at the LHC in such a 
scenario. In this respect, of particular relevance are hi aiai decays, as they have been claimed to 
be the only means to establish a no-lose theorem at the CERN machine for the NMSSM [185,293,294, 
301-307], at least over a region of parameter space where Supersymmetry (SUSY) partners of ordinary 
Standard Model (SM) objects are made suitable heavy. Here, oi states are rather light (of 10 GeV or 
less) while hi ones could well be below the LEP limit on the SM Higgs mass, of 114 GeV (albeit with 
weakened couplings to ordinary matter). The scope of hi — > aiai decays into jjr^T^ pairs (where j 
represents a jet of either heavy or light flavour and where the r's decay leptonically) has been found to be 
rather questionable [308]. Hence, in this contribution we investigate the scope of the 4r channel, wherein 
two r's are searched for in their muonic decays while the other two are selected via their hadronic ones. 
We will consider both HS and VBF as production channels. Finally, to enhance the yield of ai t^t^ 
decays, we limit ourselves to regions of NMSSM parameter space where Ma^ < 2m;, (light ai scenario). 

19.2 The Low-Energy NMSSM Parameter Space For The Light ai Scenario 

In this section we investigate the NMSSM parameter space setups which yield the M^^ < 2mb mass 
configuration, with particular interest to the cases where the aforementioned hi ^ aioi ^ 4t decays 
may be visible at the LHC if happening in conjunction with HS and/or VBF production processes of 
Higgs bosons. Notice that there are altogether fourteen parameters that uniquely define at the Electro- 
Weak (EW) scale the NMSSM Higgs sector for the purposes of our analysis. With reference to notation 
already defined elsewhere in this report, these are: tan A, k, Ax, A^, Mi, M2, M3, At, A},, Ar,Mf^ 
and Mf^, where M and Mj^ denote the soft SUSY breaking slepton and squark mass parameters. We 
will start by establishing the portion of NMSSM parameter space, defined in terms of the above inputs, 
that survives present theoretical and experimental constraints. 

1 9. 21 Full NMSSM Parameter Scan 

The numerical values over which the parameters introduced above have been scanned are: 

-1000 GeV <A^< 100 GeV, - 10 TeV < < 10 TeV, 100 GeV < /i < 1000 GeV, 

10"^ < A, K < 0.7, 1.5 < tan/3 < 50 (34) 
while the remaining parameters were fixed at 

M1/M2/M3 = 150/300/1000 GeV, At = A = = 2.5 TeV, M/^ = M/^ = 1 TeV. (35) 

We will call the scan performed over such intervals a 'wide' scan. This (as well as all those in the remain- 
der of this note) has been performed by using the NMSSMTools package [298-300], which calculates 

^'Contributed by: Alexander Belyaev, Stefan Hesselbach, Sami Lehti, Stefano Moretti, Alexander Nikitenko, Claire H. 
Shepherd-Themistocleous 
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Fig. 55: Result of the NMSSM 'wide' scan mapped onto thie planes: (a) [Ax, M^J, (b) [A^, MhJ, (c) [k, M^J, (d) [A, M^J, 
(e) [Mai 1 -^^hili (f) [Mai I "^x?]' Colour code: red - all constraints are satisfied but relic density (above WMAP constraint: 
Qh^ > 0.11); black - all constraints are satisfied, Ma > 10 GeV; green - all constraints are satisfied, Ma < 10 GeV. 

NMSSM spectra (masses, couplings and decay rates) and takes into account experimental inputs includ- 
ing LEP limits, i?-physics bounds as well as (cold) DM constraints. In Fig.|55]we present the results of 
this scan. Though only a few Af^^ < 10 GeV points survive, one can see from Fig.|55la) the preference 
for a large positive Ax while Fig. \55\ b) indicates that small | | 's are preferred. 

1 9. 22 Scan for Narrowed 

The results of Fig. [55] (specifically, the preference for small A^'s) motivated us to 'narrow' the range of 
the parameters, by scanning it over the intervals 

- 15 GeV <A^<20 GeV, - 2 TeV < < 4 TeV, 100 GeV < /i < 300 GeV, (36) 

and the rest of the parameters as in Eq. (l34l) . Fig. [56] makes the point that this is precisely the region 
where a large portion of NMSSM points with M„j < 10 GeV are found, consistent with all known 
constraints. Now we can clearly see certain correlations onsetting in the Mq^ < 10 GeV region: 1. 
values of > are preferred, see Fig. [53 a); 2. points with low M^^, ~ (Fig. [53b)) and small 
values of k (Fig. [561 c)) are preferred; 3. we can see interesting < 10 GeV points with also low, 
down to 20 GeV, Mh^ values (Fig.[56];e)). 

19.23 Final Scan For the Light ai Scenario 

We have then performed one 'final' scan over the NMSSM parameter space by requiring at the same 
time Ma^ < 10 GeV and Eq. ([36] ). Having already learnt the size of such portion of the entire NMSSM 
parameter space after experimental constraints, we now want to characterise it in terms of the quantities 
which enter the event rates for hi — > oioi ^ 4t decays produced via HS and VBF The results of this 
scan are shown in Fig. [57] Note that, here, the colours were chosen to indicate the measure of decoupling 
of the lightest CP-even Higgs boson, hi, from the SM limit (denoted simply by H). To this aim, we 
have defined the measure Rzzh = (fl'zz^f ^/•S'Izh)^' i-^-' ^^tio of the coupling strength (squared) 
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Fig. 56: Results of the NMSSM 'narrow' scan, i.e., analogous to Fig.[55]but for -15 GeV < A„ < 20 GeV, -2 TeV < 
A\ < 4 TeV, 100 GeV < fj, < 300 GeV. The individual plots and the colour code are the same as in Fig. 1551 

of the ZZh\ vertex in the NMSSM relative to the SM case (in fact, this is the same for couplings to 
gauge bosons). One should notice that both HS, Vhi, and VBF, pp jjV*V* jjhi, 

rates {V = Z, W^) are directly proportional to Rzzh and are suppressed in the non-decoupling regime 
whenever Rzzh is essentially smaller then unity. From Fig. [57] one can see the following important 
features of the Mq^ < 2mfe 10 GeV scenario: 1. the lighter the Higgs the more significant should be 
the NMSSM deviations from the SM case, e.g., for any M^^ < 50 GeV any Rzzh is limited to be < 0.5, 
as dictated by LEP constraints [4, 270] (this correlation is illustrated in a more clear way in Fig. [58la). 
presenting the Rzzh versus M^^ plane, which exhibits the typical pattern of the LEP Higgs exclusion 
curve [270]); 2. in the Mh^ < 40 GeV region Ax is always positive (Fig. [STta)), k < 0.1 (Fig.^c)) 
while A < 0.45 (Fig. [STl d)). In this case, one should notice the correlation between the singlet nature 
of the hi and the singlino component of the lightest neutralino, which is visually depicted in Fig. [58l b). 
Finally, it is also worth to point out the correlation between their masses in Fig. [58l c). From these plots, 
one can see a striking correlation between the lightest neutralino and Higgs boson whenever one has that 
Mh-^ < 50 GeV. In this connection, one should stress that the NMSSM model structure requires hi to be 
a singlet and Xi to be a singlino (for M^^ < 10 GeV and M^-^ < 10 GeV) in order to have a relic density 
consistent with current experimental constraints. In fact, over the NMSSM parameter space restricted 
to having Maj < 10 GeV and M^-^ < 10 GeV, the Xi-pair annihilation in the early Universe proceeds 
through the /ii -funnel region. So, in this region, 2m^o ~ Mh-^ as we observe from the lower-left part of 
Fig.[58tc). 

19.3 Phenomenology of the Light ai Scenario 

As final step of our analysis, we combined the production rates of HS and VBF with selection efficiencies 
determined by generating these processes within the PYTHIA Monte Carlo (including smearing effects). 
The latter have been estimated in presence of cuts, after parton shower, hadronisation plus heavy hadrons 
decays (and with underlying event turned on). For HS we enforced (assuming e, /x decays of the W^) 
r+ —>■ livv and hadrons and the selection cuts were 
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Fig. 57: Results of the NMSSM 'finar scan, i.e., with Ma^ < 10 GeV and with Eq. l l36t enforced. The black, blue and green 
colours indicate the cases Rzzh < 0.1, 0.1 < Rzzh < 0.5 and Rzzh > 0.5, respectively (where Rzzh is defined in the 
text). The individual plots are the same as in Fig. 1551 
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Fig. 58: Results of the NMSSM 'final' scan, i.e., with Ma^ < 10 GeV and with Eq. 
Mh, , (b) Rzzh ve 
same as in Fig. 1571 



• Trigger selection: isolated single muon or single electron found with thresholds 19 and 26 GeV, re- 
spectively, \ri\ < 2.5. • Muon px > V GeV, |r/| < 2.1. • Tau jet Et > 10 GeV, Ivl < 2-1- • Isolated 
1 -prong r's within Ai? < 0.6 from the muon using tracker isolation for tracks pT > 2 GeV. • Tau and 
muon oppositely charged. • Two tau -i- muon pairs found. 

For VBF the selection cuts were 

• Two same sign muons with px > 7 GeV, \r]\ < 2.1 and with one track of px > 2 GeV in cone 0.6 
around each muon. • Two r jets with Et > 10 GeV, \r]\ < 2.1. • Two jets with Et > 30 GeV, \rj\ < 4.5. 

The results in Fig. |59] show that, after our final scan, the population of parameter points is such that in 
both channels the highest cross sections are found for Mh-^ > 80 GeV, although in the case of VBF also 
lower hi masses can yield sizable rates, but never for values less than 40 GeV Independently of M/j^, 




Fig. 59: Cross sections (including all relevant branching ratios) for HS (a) and VBF (b) after the selection cuts described in the 
text. The population of points used correspond to that of the 'final scan' described previously. Black diamonds correspond to 
the benchmark P2 (right) and P3 (left) from Ref. [283]. 

the ai mass enables sizable event rates anywhere above 2mT-, but particularly just above the threshold. 
At high luminosity, 100 fb~^, the highest rates would correspond to 1000 events per year for HS and 
8000 for VBF. 

19.4 Conclusions 

We have shown that there is significant potential in establishing a no-lose theorem for the NMSSM at 
the LHC via (marginally) HS and (primarily) VBF production of the lightest CP-even Higgs boson hi 
decaying into aiai pairs in turn yielding four r leptons, searched for through their semi-leptonic/hadronic 
decays into muons and jets. To enhance the decay fraction into r's of the lightest CP-odd Higgs boson 
ai we have restricted ourselves to the case Ma^ < 2m{, (otherwise ai —>■ bb decays ai^e dominant). We 
have also found that the hi state can be very light, indeed at times lighter than the ai. However, this 
last configuration can only be achieved in a low-energy NMSSM setup, with no unification assumptions 
at the high scale. In fact, we are currently investigating whether such light hi masses can be found in 
a less constrained version of the cNMSSM discussed in [283]. Finally, with reference to the NMSSM 
benchmarks discussed in [283], we should like to point out here that those relevant to our 4r channels 
are P2 and P3. We have reported the cross section times efficiency rates for these two points in Fig. [59] 
(black diamond symbols, P2 to the right and P3 to the left). As it can be appreciated, they correspond to 
event rates that are mid range amongst all those explored, hence not particularly biased towards a far too 
favourable NMSSM setup, yet susceptible to experimental discovery. Our summary is preliminary, as 
only signal processes have been considered and only in presence of MC simulations, with no backgrounds 
and detector performance enabled. The latter clearly ought to be investigated before drawing any firm 
conclusions and this is currently being done. 
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20. INVESTIGATION OF THE LHC DISCOVERY POTENTIAL FOR HIGGS BOSONS IN 
THE NMSSM 

20.1 Introduction 

The Large Hadron Collider (LHC) will deliver proton-proton collisions at a center-of-mass energy 
of 14 TeV. First physics runs are expected for 2008. First, the LHC will operate at low luminosity 
(2 • 10'^^cm~^s~^). Later, the luminosity will be increased to its design value of 10'^^cm~^s^^. One of 
the main aims of the ATLAS [309] and CMS [98] experiments at the LHC is the search for the Higgs 
boson. In the Standard Model (SM) electroweak symmetry breaking is achieved via the introduction of 
one Higgs doublet. Only one neutral Higgs boson is predicted. Extended Higgs sectors, with additonal 
Higgs doublets and Higgs singlets give rise to several neutral and charged Higgs bosons, e.g. the two 
Higgs doublets of the Minimal Supersymmetric Extension of the SM (MSSM) yield three neutral and two 
charged Higgs bosons. Detailed studies have shown that the SM Higgs boson will be observable at AT- 
LAS and CMS [98, 108, 309]. The discovery of one or more Higgs bosons of the CP-conserving MSSM 
will be possible [310]. Previous studies claim that at least one Higgs boson of the Next-to-Minimal Su- 
persymmetric Standard Model (NMSSM) will most likely be observable at the LHC [293,302]. Here, 
we present an evaluation of the discovery potential for NMSSM Higgs bosons based on current ATLAS 
studies [108,263,309,311-319]. 

20.2 The NMSSM Higgs Sector 

In the framework of the NMSSM, the fi-problem of the MSSM is solved by the introduction of an 
additional neutral singlet superfield S [290]. The two additional neutral scalar bosons contained in S 
mix with the MSSM Higgs bosons to form the five neutral Higgs bosons of the NMSSM: three CP-even 
bosons Hi, H2, and two CP-odd Higgs bosons Ai, A2. The phenomenology of the charged Higgs 
boson is only modified marginally with respect to the MSSM. The Higgs sector of the NMSSM at 
Born level is determined by the four coupling parameters of the singlet superfield. A, k, Ax, A,^, and 
the two parameters fi and tan/3. For a more detailed description of the NMSSM Higgs sector see e.g. 
Refs. [290,320]. 

20.3 Evaluation of the Discovery Potential 

Two two-dimensional benchmark scenarios are investigated in this study: the Reduced Couplings Sce- 
nario and the Light Ai Scenario which were proposed during this workshop (for details see these pro- 
ceedings). The parameters A and k are varied in meaningful ranges whereas the other parameters are 
fixed as described previously in this report. The method of evaluation of the discovery potential is simi- 
lar to the study performed for the MSSM in Ref. [310]. 

20.31 Calculation of masses and events rates in the NMSSM 

NMHDECAY [298, 299] was used to calculate the masses, branching ratios and decay widths of the 
NMSSM Higgs bosons and the couplings of the neutral Higgs bosons to fermions and gauge bosons, 
relative to the respective SM couplings. Couplings to gluons relative to the SM couplings were calculated 
from the ratio of partial widths of H^gg in the NMSSM and the SM [188] as in Eqi37l 



For the neutral Higgs bosons, leading order SM cross sections [115] were scaled according to Eqj38l 



9Hgg,NMSSM _ ^{H ^ 99) NMSSM 
SHgg,SM ^ 99)sM 



(37) 



cnmssm = crsM ■ 



9nmssm 

9sM 



(38) 



36, 
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Table 20: Included search topologies with allowed mass ranges. 



Search Channel 


Mass Range [GeV] 


Refs. 


VBF, H^TT 


110-150 


[108] 


VBF, H^WW^lluv 


110-200 


[108] 


VBF, H^WW^lvh 


130-200 


[108] 


VBF, i7^77 


110-160 


[311] 


ttH, H^hh 


70-150 


[312] 


GGF, H^ZZ-^Al 


120-420 


[309] 


GGF, H^WW^llvv 


140-200 


[313] 


WH, H^WW^llvv, W^lz^ 


150-190 


[309] 


Inclusive H^jj 


70-160 


[309] 


Inclusive A—>-jj 


200-450 


[309] 


WH, ZH, ttH, iJ^77 


70-150 


[309] 


bbH, F/yl^rr^hh 


450-800 


[314] 


GGF, bbH, H/A^TT^luh 


150-800 


[315] 


GGF, bbH, H/A^fifi 


70-450 


[316,317] 


GGF, H^hh^-f-fhh 


230-270/70-100 


[309] 


GGF, H^ZA^Ubb 


200-250/70-100 


[309] 




175-600 


[263] 


gb^HH, H^^tb 


190-400 


[318] 


tt^HHWb^Tulubb 


90-165 


[309] 


ti^HHWb^Tuqqbb 


80-165 


[319] 



The charged Higgs boson gb—^tH^ cross sections in leading order were taken from Ref. [244] and were 
modified according to the i^^tb-couplings obtained with NMHDECAY. The branching ratio t^H^b 
was calculated with Feynhiggs [246]. For tt-production, a leading order cross section of 482 pb was 
assumed. The top quark mass was set to 172 GeV. Theoretical and LEll^ exclusion criteria (bounds 
from hZ and hA searches) were calculated by NMHDECAY. 

20.32 Significance Calculation 

The expected number of signal events is derived from the above discussed NMSSM cross sections. Signal 
efficiencies are taken from published ATLAS Monte-Carlo studies (Table l20l). The expected numbers of 
background events are also taken from published ATLAS MC studies. If MC studies at design luminosity 
exist, a data volume of 300 fb~^ is assumed; if only low luminosity studies are available, 30 fb"^ are 
used, and if both scenarios have been investigated, 30 fb^^ taken at low luminosity and 270 fb^^ taken at 
design luminosity are assumed. The current results only include SM background processes. Systematic 
uncertainties are neglected. For the significance calculation, the profile likelihood method [321] with 
asymptotic approximation [322] is used. To claim a discovery, a significance of at least 5a is required. 
The number of expected signal events is corrected for the effects of increased Higgs boson decay widths 
and the possibility of degenerate Higgs boson masses as described in the following. 

Corrections for large Higgs bosons widths 

In the NMSSM, the natural line width of the Higgs boson may be enhanced relative to the SM case. 
Thus, a larger fraction of signal events may lie outside a mass window cut than in the SM. To correct for 
this, the Higgs boson peak was described by a Voigt-function whose Breit-Wigner part is given by the 
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natural line width, the Gaussian part by the detector resolution. The ratio of the integral values over the 
mass window for the SM and the NMSSM case was used as a correction factor. 

Corrections for degenerate Higgs boson masses 

Higgs boson peaks were described by a Voigt function as previously. The peaks were assumed to be 
indistinguishable if their separation was smaller than 2.355-FWHM. In case of negligible Higgs boson 
width, this corresponds to a 2a separation of two Gaussians. Higgs bosons with overlapping mass win- 
dows were also considered indistinguishable to avoid double counting of events. In case of inseparable 
peaks, contributions from all Higgs bosons were added up for each boson's mass window. Only the high- 
est observed significance was kept and assigned to the Higgs boson with the largest fraction of signal 
events in that mass window. 

20.4 Search Topologies 

The combinations of production mechanisms and decay modes considered in the evaluation of the disco- 
very potential and the considered mass ranges are summarised in Table [200 Within the scenarios exam- 
ined here, only the VBF, H—^tt; ttH, H—^bb and i/— ^77 channels show significances greater 5o" at the 
given integrated luminosities in the theoretically allowed and yet unexcluded regions (see section |203] ). 

20.5 Results 

20.51 The Reduced Couplings Scenario 

In the Reduced Couplings Scenario, the H2 with a mass of about 120 GeV is SM-like in large parts of 
the parameter space. In an unexcluded region with large negative k close to the lower exclusion bound, 
the gauge couplings of H2 are reduced to about 80% of their SM-value. The Hi gets very light at 
the region close to the upper exclusion bound, so that the decay H2^HiHi is kinematically allowed. 
However, due to the small branching ratio for this decay mode of at maximum 6%, its effect on the 
discovery potential is negligible. The discovery potential for the H2 is shown in Figl60l The entire 
unexcluded region is covered by the ttH, H2^bb channel despite the coupling reduction. The inclusive 
ffg— >77 and the VBF, H2^tt channels also contribute. With 30 fb^^, the search for H2^tt will be 
the only sensitive channel. The region with reduced couplings will not be covered in that case. The 
gauge couplings of the Hi and H-^ are sizable only at large negative k. Here, the channels //s— ^77; 

Production modes are abbreviated GGF for gluonfusion, VBF for vector boson fusion and ttH, bbH, WH and ZH for 
associated production with top quarks, bottom quarks and vector bosons. 




Fig. 60: 5(T discovery contours of the H2 in the X-k plane for 
the Reduced Couplings Scenario 



Fig. 61: 5a discovery contours of the Hi and Hs in the X-k 
plane for the Reduced Couplings Scenario 




VBF, Hs^tt; ttH, Hi/s^bb and GGF, H^^ZZ^Al contribute in a region ruled out by LEP (FigllB- 
Since the charged Higgs boson is lighter than the top quark in the same region, it can be detected via 
the tt^H^hW^h^Tvlvbh and ti—>-H^bW^b—>-Ti>qqbb searches only in the LEP excluded region also 
(Figl62l). All other Higgs bosons have highly reduced gauge couplings and ai^e therefore unobservable. 

20.52 The Light Ai Scenario 

In this scenario, the Hi has a mass of about 120 GeV and SM-like gauge couplings. Since the Ai is light, 
Hi^AiAi decays are kinematically possible and often dominant. In the upper right unexcluded region, 
the branching ratio of Hi^AiAi is larger than 90%. Here, the Hi cannot be observed (see Figi63l). The 
branching ratio of Hi^AiAi drops for small A and k. Therefore, a discovery via the inclusive and asso- 
ciated Hi^^j; VBF, Hi^TT and ttH, Hi^bb modes is possible in that region (Figl64l). The outermost 
discovery contour of Hi^jj follows approximately the 60% line of the branching ratio of Hi^AiAi. 
The H2 has contributions from the channels i?2^77; VBF Hi^WW; GGF, H2^ZZ^4l and GGF 
H2^WW —^2l2u in the excluded region where it is light enough to be accessible (Figj65]). All other 
Higgs bosons have either highly reduced couplings or are too heavy to be observed in this scenario. 
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Fig. 64: 5(7 discovery contours of the Hi in the A-k plane for 
the Light Ai Scenario, restricted to low A and k values. 



Fig. 65: 5(7 discovery contours of H2 in the \-k plane for the 
Light Ai Scenario, restricted to low A and k values. 



20.6 Conclusions 

An evaluation of the ATLAS discovery potential for NMSSM Higgs bosons within two benchmark sce- 
narios was performed. At least one Higgs boson was found to be observable in regions without a light 
Ai or where the branching ratio of Hi^2~^AiAi is smaller than about 60%. In the other cases, searches 
for the decay chains Hi/2—>^AiAi—>-TTbb or Hii2^AiAi^At could be considered. 
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21. THE A^A^ 66r+T SIGNAL IN VECTOR BOSON FUSION PRODUCTION AT 

THE LHC 

21.1 Motivation 

In the Minimal Supersymmetric extension of the Standard Model (MSSM) at least one additional 
SU{2)l Higgs doublet is required compared to the SM in order to cancel gauge anomalies of the super- 
partners and to allow Yukawa couplings for all fermions. In order to address the fine-tuning "^u-problem" 
that appears in the MSSM, one can also add an extra complex singlet to these doublets. This last possi- 
bility, known in the literature as the Next-to-Minimal Supersymmetric Standard Model (NMSSM), has 
an interesting phenomenology (e.g. see [323]). 

In the NMSSM, one of the pseudoscalar states {A^) is the Goldstone boson of either a global 
[/(I) R-symmetry or a U{\) Peccei-Quinn symmetry in some limit of the model parameters. Since 
low-fine tuning scenarios predict a moderate breaking of these symmetries, the mass of A^ is expected 
to be relatively small compared to the mass of the lightest scalar (/i*^) such that the A^A^ decay 

is kinematically allowed. In [324], two different types of scenarios are considered, depending on if 
m^o > 2mij or m^o < 2mfe . Scenarios with m^o > 2mf, are disfavored when LEP data for Z2b and 
ZAb final states are taken into account. Indeed, the simultaneous analysis of both these channels excludes 
at better than 99% the possibility for to be lighter than ~ 108 GeV, and a heavier in turn requires a 
higher fine-tuning of model parameters. On the contrary, scenarios with m^o < 2mi, are favored by the 
same data and can even account for the 2a excess observed in the Z2b final state in the m^^o ~ 100 GeV 
vicinity. As a consequence, many NMSSM related analyses focus on the —>■ A^A^ — > t~^t~t^t~ 
decay which has the most favorable branching ratio if rrij^o < 2m;,. 

Nevertheless, besides the particular context of the NMSSM, many other possibilities remain open. 
If the MSSM scalar sector violates the CP symmetry, standard mass relations do not hold anymore and 
the decay of into two lighter Higgs bosons may be allowed [325]. In [326], a light A^ (i.e., between 
0. 1 and a few tens of GeV) decays predominantly into pairs of photons (or gluons) thanks to a vector-like 
quark loop. Another possibility arises in the context of the generic two-Higgs-doublet model (2HDM). 
As shown in [327], a moderately light A^ (i.e., between 10 and 100 GeV) can naturally satisfy the p 
parameter constraints thanks to a twisted realization of the custodial (or equivalently CP) symmetry. As 
emphasized in [328], a pseudoscalar in this mass range together with a moderate value of tan /? can also 
account (in type II 2HDMs) for the observed discrepancy between the experimental measurement of the 
muon anomalous magnetic moment and the SM predictions. 

Assuming m^o > 2mb, and that the coupling of A^ to fermions is proportional to the mass for 
down-type quarks and leptons (the up-type quark couplings being negligible for tan /? » 1), the main 
decay modes are A^ bb (BR ~ 0.92) and A^ t^t~ (BR ~ 0.08). Under the hypothesis that 
BR(/i° A°^°)~ 1 (which is a reasonable approximation in many models), this gives a total branching 
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ratio of ~ 0.85 for /i° 46, ~ 0.15 for /i° A°A° 262r and less than one percent for 

/i*^ —>■ A^A^ —>■ At. Since the four r final state signal is suppressed at least by a factor of a hundred 
compared to the m,y^o < 2mh scenario studied in Sectior fT9l the LHC discovery of hP and AP in this 
channel is probably difficult. On the other hand, the four h final state has a large BR, but suffers from 
important QCD backgrounds. This final state has been investigated in direct production mode at the 
Tevatron (where it is overwhelmed by the backgrounds [329]) and in W/Z associated production [307]. 
At the LHC, a discovery significance may still be reached in this last mode [307, 330]. 

In the current work, we focus on the intermediate 262r final state, which has a smaller but still 
sizable BR than the 46 final state, together with a much lower background. This final state has been 
considered in the framework of the associated production of hP with a W/Z boson at the Tevatron 
in [100, 307]. However, in this case, only a few events could be observed after a few fb~^ due to 
the cuts and 6/r tagging necessary to remove the large reducible background. Similar difficulties with 
the reducible background are also expected at LHC [307]. In the present study, we concentrate on 
the Vector Boson Fusion (VBF) production mode for hP, which has been shown to be a promising 
channel at the LHC for the SM decay hP t^t^ both in parton-level analysis [113, 331] and after full 
detector simulation [108, 332, 333]. After the end of the redaction of this work, it has been brought to 
our knowledge that a study on similar lines in the context of the NMSSM, using parton shower based 
simulations, can be found in [293, 304]. 

21.2 Signal and Background 

The signal and background Monte-Carlo simulation has been carried out at tree level using Mad- 
Graph/MadEvent v4 [153] for the parton-level event generation. 

In the framework of this preliminary analysis, some simplifying assumptions are made. The SM- 
like Higgs, hP , shares all SM Higgs boson couplings plus an additional coupling to the pseudoscalar AP 
large enough to ensure BRC/i*^ — > A^A^) ~ 1. Its mass is fixed at 120 GeV, i.e. this is above the best 
LEP limit to avoid de facto all possible direct constraints, but is still light enough to ensure a sizable 
production cross-section. The light pseudoscalar mass is fixed at 50 GeV in order to lie below the m^o /2 
threshold, while still being large enough to guarantee a good angular separation of decay products. 

As mentioned in the previous section, the coupling of AP to fermions is assumed to be proportional 
to their mass for down-type quarks and charged leptons, giving a total branching ratio for hP —>■ A'^A'^ —>■ 
262r of about 0.15, which may be compared with the SM expectation BRC/i*^ —>■ r+r^) of ~ 0.08. This 
is only true if the coupling to up-type quarks is strongly suppressed, for example, due to an additional 
tan P factor in a type II 2HDM. If this is not the case, the considered total branching ratio can be reduced 
by up to a factor two. 

In order to improve signal to background separation, few kinematical cuts such as minimum px 
of 10 GeV for all 6-jets, 20 GeV for all non-6 jets, and 10 GeV for all leptons have been applied. In 
addition, acceptance cuts such as the maximum pseudorapidity of 5 for jets, and of 2.5 for 6-jets and 
leptons, and a minimal separation cut, i.e. AR > 0.3, on all objects pairs have been applied at the 
parton level. Furthermore, to narrow ourselves to the particular kinematic configuration of signal events, 
standard VBF cuts are applied, i.e. | Ary| > 4 and mjj > 700 GeV for the two forward jets. Finally, 
a maximum invariant mass cut, < 80 GeV, is imposed on all leptons pairs to avoid the Z peak in 
some backgrounds. 

The signal is characterized by a final state populated with two central 6 jets, two central r's and 
two forward jets. To avoid triggering issues, we focus on final states in which both r's decay leptonically. 
The associated tree level cross-section (after r decays and cuts) is 9.5 fb. The irreducible background, 
where the r pair is coming from an off-shell photon or Z, and the 6 pair from a gluon splitting, is rather 
low, with a 1 fb cross-section. The same process with an e or /i pair replacing the r pair has a more 
sizable cross section of around 8.7 fb, due to the absence of the r branching ratio suppression. The most 




Fig. 66: Invariant mass of the final-state (oppositely charged) di-lepton pairs and di-bottom quarks before the final kinematical 
cuts[39l Each signal or background contribution is normalized by cross-section. 
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Fig. 67: Invariant mass Mbui of the four-body final state. The signal and background histograms are stacked and normalized 
by their corresponding cross-sections. 



dangerous background is ti + 2jets with fully leptonic top decays (through an intermediate r or not). 
However, even though the total cross-section is almost three orders of magnitude larger than the signal 
(3.2 pb), the associated kinematics, and in particular the invariant mass distribution of 6's and leptons, 
are very different. 

21.3 Results 

Figure[66lshows the invariant mass combination of any oppositely charged di-lepton pairs and any bottom 
quark pairs. Only the cuts described in the previous section have been applied. 

Looking at the kinematic distributions of the signal and background samples (described in the 
previous section), it is evident that a cut based technique can be defined to achieve separation. The 
chosen selection criteria are: 

Mil < 30, 40 < Mbb < 60, ARii < 2, and ARbb < 2. (39) 

Figure [67] shows the invariant mass Mbui, of the four body final state after these simple cuts. The 
signal and the background considered are stacked and normalized by cross-section. A simple estimate 
of the significance around Mbui, in the region 50 < M^bii < 110, yields S/^/B = 4 for an integrated 
luminosity of 100 fb^^, with approximatively 100 signal events. B-tagging efficiency will impact the 




number of both signal and background events, and reduces this significance by a factor of ~ 2 if a b- 
tagging efficiency of 50% is assumed. Of course this simple generator-level estimate is merely a crude 
check on the feasibility of studying A^A^ ^ bbr^r^ in VBF. 

Since after reconstruction we expect the significance to decrease even further, this parton-level re- 
sult motivates the use of a statistics-based multivariate approach in order to further discriminate between 
signal and background. Preliminary results demonstrating the discriminating power of the technique be- 
tween the signal, the irreducible background and part of the ti + 2j background are shown in Figure l68l 

Figure |6l] (a) shows the relative contribution of the various input variables to signal and back- 
ground separation. A framework for parameter space optimization, PARADIGM, is utilized for the above 
task [334]. The two most effective variables for signal/background separation in this decay mode are 
the invariant mass of the 6-jets and leptons M^fi^u and the invariant mass of the 6-jets A/{,f„ as was also 
observed in the cut based study. Although PARADIGM allows the reduction of parameter space, we do 
not eliminate any of the variables since the dimensionality of the initial feature space considered is lower 
than the degrees of freedom of the model. Therefore, it is likely that the classifier performance can be 
further enhanced by the addition of more variables. 

The decision tree classifier output is shown in Figure[68](b). The measure of discrepancy between 
the background-only hypothesis and the background plus signal hypothesis (assuming a normal error 
distribution and using the classifier itself as the test statistic in a two-tailed test [335]) is found to be 
0.0086 ± 0.0058 at 95% CL. This is a statistically significant result. 

21.4 Conclusions 

We showed that the — > A^A^ — > bbr^T^ signal in VBF production at the LHC is potentially feasible 
with an integrated luminosity of 100 fb^^. Using a simple cut based technique, we found approximately 
25 signal events and a significance of ~2 for this luminosity (taking into account a 50% ^-tagging effi- 
ciency) . This result motivates the use other techniques, such as a multivariate analysis, to further enhance 
the feasibility of this search at the LHC. A more robust multivariate analysis that includes different mass 
hypotheses, a full set of reducible backgrounds as well as fast detector simulation and evaluation of 
systematic uncertainties is envisaged by the authors. 
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